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ABSTRACT 


In  this  thesis,  we  develop  the  forward  channel  model  for  a  DS-CDMA  cellular 
system  operating  in  a  slow-flat  Rayleigh  fading  and  lognormal  shadowing  environment, 
which  incorporates  the  extended  Hata  model  to  predict  median  path  loss.  Forward  error 
correction  is  integrated  into  the  model  by  applying  convolution  encoding  with  soft- 
decision  decoding.  The  worst-case  probability  of  bit  error  for  a  mobile  user  at  the  edge 
of  the  center  cell  of  a  seven-cell  cluster  is  developed  using  Gaussian  approximation.  In 
estimating  the  probability  of  bit  error,  we  develop  a  statistical  model,  which  approximates 
the  sum  of  d  multiplicative  chi-square(two  degrees  of  freedom)-lognormal  random 
variables  as  a  multiplicative  chi-square(with  2 d  degrees  of  freedom)-lognormal  random 
variable.  Using  this  approximation,  we  examine  the  performance  of  the  cellular  system 
under  a  range  of  shadowing  conditions,  for  various  user  capacities  and  with  antenna 
sectoring  as  they  compare  with  Monte  Carlo  simulated  results.  Next,  we  modify  our 
worst-case  performance  analysis  to  accommodate  users  that  are  distributed  in  the  cell 
according  to  a  specified  distribution  and  compare  results  with  the  worst-case 
performance.  Finally,  we  introduce  fast  power  control  into  the  forward  channel  and 
explore  system  performance  with  power  control  under  a  range  of  operating  conditions  as 
it  compares  with  the  fixed-power  performance. 
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EXECUTIVE  SUMMARY 


In  the  near  future,  the  field  of  wireless  mobile  cellular  communications  will  enter 
the  domain  of  high-speed  Internet  access  with  the  deployment  of  the  third-generation 
cellular  system  (3G).  While  the  second-generation  systems  currently  in  use  are  digital 
and  support  data-type  traffic  at  low  data  rates,  they  were  built  to  carry  principally  voice- 
type  traffic.  The  third-generation  system  is  designed  to  support  data-type  traffic  at  much 
higher  data  rates  than  the  second-generation  systems.  Specifically,  manufacturers  are 
promising  download  rates  on  the  order  of  2.4  Mbps:  roughly  twice  as  fast  as  wired 
broadband.  For  fast-traveling  users,  the  expected  data  rate  is  reduced  to  144  kbps:  better 
than  a  modem  over  a  conventional  phone  line.  The  actual  performance  of  third- 
generation  systems,  however,  is  yet  to'be  determined.  Accordingly,  this  body  of  work  is 
aimed  at  building  models  and  tools  for  predicting  and  enhancing  the  performance  of 
cellular  systems  such  as  3G. 

The  third-generation  cellular  system  will  employ  a  multiple  access  technique 
known  as  Direct  Sequence  Code  Division  Multiple  Access  (DS-CDMA).  Multiple  user 
access  is  a  key  feature  in  cellular  system  operations  since  it  provides  for  the  efficient  use 
of  assigned  bandwidth  to  accommodate  many  users  simultaneously.  DS-CDMA  systems 
assign  mobile  users  orthogonal  Walsh  codes  and  spread  the  frequency  spectrum  of  the 
user  traffic  using  the  designated  Walsh  code.  All  user  traffic  is  spread  across  the  same 
frequency  band  and  transmitted  simultaneously  in  a  composite  forward  signal  from  the 
base  station.  Mobile  users  extract  their  traffic  from  the  composite  signal  using  their 
assigned  code.  In  this  thesis,  we  explore  the  algebraic  group  properties  of  a  set  of  Walsh 
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codes,  or  Walsh  functions.  We  introduce  the  idea  of  extending  the  orthogonality  property 
of  Walsh  functions  to  ensure  orthogonality  between  any  set  of  three  Walsh  functions. 
Finally,  we  examine  the  autocorrelation  functions  for  a  set  of  Walsh  functions  and  we 
develop  an  expression  for  the  average  autocorrelation  for  a  set  of  Walsh  functions. 

This  thesis  specifically  addresses  the  forward  channel  of  a  DS-CDMA  cellular 

system.  The  forward  channel  carries  traffic  from  the  base  station  to  a  wireless  mobile 

user  in  the  cell.  The  performance  of  the  forward  channel  is  critical  in  third-generation 

systems  since  it  will  be  used  to  download  data  from  the  Internet.  Modeling  the  forward 

channel  is  complicated  by  the  fact  that  the  receiver  is  mobile.  In  general,  the  signal 

transmitted  by  the  base  station  to  the  mobile  user  will  lose  a  portion  of  its  power  based  on 

the  distance  it  must  travel  and  the  terrain  over  which  it  travels.  Accordingly,  we  predict 

the  average  power  received  by  a  mobile  user  at  a  specific  distance  using  large-scale 

propagation  models.  However,  mobile  users  at  the  same  distance  could  have  terrain  in 

their  immediate  vicinity  that  is  vastly  different  and  dynamically  changing  due  to  their 

mobility.  This  phenomenon,  known  as  shadowing,  results  in  actual  path  loss  being 

different  for  users  that  are  at  the  same  distance  from  the  base  station.  We  reflect  this 

variability  in  path  loss  using  a  lognormal  random  variable  in  our  channel  model.  In 

addition  to  the  large-scale  propagation  effects,  a  mobile  user  typically  receives  multiple 

copies  of  the  forward  signal  that  have  been  reflected  off  obstacles  along  various  paths. 

The  relative  motion  between  the  mobile  and  the  base  station  induces  a  doppler  shift  on 

each  of  the  multipath  components.  Furthermore,  the  obstacles  that  reflect  the  multipath 

components  may  be  moving,  which  also  causes  a  time-varying  doppler  shift  of  the 

received  signal.  The  combined  effects  of  multipath  propagation  and  the  doppler  shifts 
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result  in  rapid  fluctuations  in  the  strength  of  the  received  signal,  which  is  called  small- 
scale  fading.  We  represent  the  effects  of  small-scale  fading  in  our  channel  model  as  a 
Rayleigh  random  variable.  Accordingly,  we  combine  the  effects  of  large-scale  lognormal 
shadowing  and  small-scale  Rayleigh  fading  on  the  forward  signal  into  a  single  unified 
channel  model,  which  we  use  to  analyze  the  performance  of  the  system  under  various 
operating  conditions. 

In  DS-CDMA  cellular  systems,  the  total  signal  that  a  mobile  user  will  receive 
includes  the  forward  signal  from  his  base  station  and  from  the  base  stations  of  the  six 
adjacent  cells.  In  our  analysis,  we  initially  examine  the  worst-case  scenario  by  placing 
the  receiving  mobile  user  at  the  edge  of  the  hexagonal  cell  in  the  center  of  a  seven-cell 
cluster.  By  doing  so,  our  mobile  user  receives  the  forward  signal  from  his  base  station  at 
roughly  the  same  power-level  as  that  of  the  forward  signal  from  the  two  nearest  adjacent 
cells.  The  forward  signals  from  the  six  adjacent  cells  constitute  intercell  interference. 
The  mobile  user’s  receiver  must  overcome  the  propagation  effects,  the  intercell 
interference,  and  additive  noise,  in  order  to  properly  decode  the  intended  traffic.  To 
improve  the  chances  of  accurately  decoding  the  signal,  we  incorporate  Forward  Error 
Correction  (FEC)  into  our  model  in  the  form  of  convolutional  encoding  with  soft- 
decision  decoding.  Furthermore,  we  include  sectoring  in  our  analysis,  which  reduces 
intercell  interference  by  splitting  the  cell  into  sectors  and  utilizing  highly  directional 
antennas.  Finally,  we  develop  the  probability  of  bit  error,  a  measure. of  system 
performance,  for  the  forward  channel  of  a  DS-CDMA  cellular  system  operating  in  a 
Rayleigh  fading  and  lognormal  shadowing  environment  that  employs  FEC  and  sectoring. 
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We  examine  the  performance  of  the  forward  channel  under  a  range  of  operating 
conditions. 

In  order  to  obtain  a  less  pessimistic  view  of  performance  of  the  forward  channel, 
we  modify  the  model  by  distributing  the  mobile  users  within  the  cell  according  to  a 
probability  distribution,  rather  than  fixing  the  mobile  at  the  worst-case  position.  We 
compare  the  distributed  user  results  with  that  of  the  worst-case  scenario. 

In  the  analysis  described  above,  we  assume  that  the  power  in  the  forward  signal  is 
fixed  for  every  channel  and  for  all  base  stations.  We  modify  our  analysis  to  allow  for  fast 
power  control  on  the  forward  channel  in  order  to  enhance  performance  of  the  system. 
Fast  power  control  ensures  that  the  power  received  by  all  mobile  users  is  at  a  target 
power  level,  which  is  the  same  for  all  users.  Based  on  feedback  from  the  mobile  users, 
the  base  station  adjusts  the  power  in  each  user  channel  to  that  which  is  necessary  to 
achieve  the  target  level.  We  develop  the  probability  of  bit  error  for  the  DS-CDMA 
system  with  fast  power  control  operating  in  a  Rayleigh-lognormal  channel  with  EEC  and 
sectoring.  We  compare  the  performance  results  of  the  system  with  power  control  to  those 
of  the  fixed  power  system. 

The  forward  channel  model  described  above  incorporates  many  of  the  system 
features  that  are  expected  to  be  included  in  the  third-generation  cellular  systems,  such  as 
EEC,  sectoring  and  fast  power  control.  More  importantly,  it  accounts  for  the  combined 
propagation  effects  of  Rayleigh  fading  and  lognormal  shadowing,  which  will  dominate 
the  operating  environments  of  the  third-generation  cellular  systems.  Accordingly,  we 
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expect  that  this  research  will  aid  system  designers  in  the  deployment  of  and  future 
upgrades  to  the  3-G  cellular  system. 
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I.  INTRODUCTION 


A.  BACKGROUND 

At  the  dawn  of  the  twenty-first  century,  wired  internet  users  have  come  to  expect 
instantaneous  access  to  e-mail,  voice-mail,  news,  weather,  music  and  video 
teleconferencing  in  the  form  of  text,  web  pages,  streaming  audio  and  video.  Wireless 
manufacturers  are  working  diligently  to  bring  a  subset  of  these  wired  Internet  features  to 
mobile  wireless  users  through  the  cellular  telephone  systems.  For  the  moment,  second- 
generation  cellular  systems  predominate  around  the  world  and  employ  various  types  of 
multiple  access  techniques.  While  the  second-generation  systems  are  digital  and  can 
accommodate  data  services,  they  were  principally  designed  for  voice-type  traffic. 
Accordingly,  second-generation  systems  typically  offer  a  meager  10  kbps  in  data  transfer 
capability  [1].  This  data  transfer  rate  frustrates  mobile  web  users  since  it  is  significantly 
slower  than  a  modem  with  a  conventional  phone  line:  the  slow-lane  on  the  wired 
information  highway.  Accordingly,  wireless  manufacturers  are  working  on  ways  to 
improve  the  data-transfer  capability  of  the  second-generation  systems.  The  transition  of 
second-generation  cellular  systems  to  accommodate  faster  data  rates  is  a  stopgap  measure 
pending  the  worldwide  deployment  of  the  third-generation  cellular  system  (3G.) 

The  third-generation  system  is  designed  to  support  circuit  and  packet  data  at 
higher  data  rates  than  the  second-generation  systems  by  employing  Wideband  Code 
Division  Multiple  Access  (W-CDMA)  technology  [2],  Manufacturers  are  promising 
download  rates  on  the  order  of  2.4  Mbps.  For  highly  mobile  (vehicular)  users,  however. 
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the  expected  data  rate  is  reduced  to  144  kbps  [3].  The  actual  performance  of  third- 
generation  systems  remains  to  be  seen.  The  deployment  of  third-generation  cellular 
systems  in  Japan  is  underway.  In  Europe  third-generation  cellular  systems  are  expected 
to  be  on-line  in  2002,  and  in  the  United  States  (possibly)  by  2003.  These  deployment 
dates  are  targets  for  metropolitan  centers,  while  worldwide  deployment  would  come 
much  later  (2010),  if  ever. 

B.  OBJECTIVE 

As  system  designers  gear  up  for  the  deployment  of  the  third-generation  systems, 
they  will  require  design  tools  and  models  that  reflect  the  mobile  radio  channel 
environments  that  they  will  encounter.  Furthermore,  as  we  search  for  ways  to  improve 
3G  system  performance,  tools  that  model  the  environment  precisely  will  become 
increasingly  more  important.  This  thesis  provides  such  tools  for  Direct  Sequence  CDMA 
(DS-CDMA)  cellular  systems  operating  in  fading  and  shadowing  environments  that  are 
associated  with  mobile  wireless  communications.  Specifically,  we  develop  a  tight  upper 
bound  on  the  probability  of  bit  error  for  the  forward  channel  of  a  CDMA  cellular  system 
in  a  slow  flat  Rayleigh  fading  environment,  which  is  also  affected  by  lognormal 
shadowing.  In  the  performance  analysis,  we  account  for  forward  error  correction  in  the 
form  of  convolutional  encoding  with  soft  decision  decoding.  Furthermore,  we 
incorporate  other  performance  enhancing  techniques  such  as  antenna  sectoring  and  power 
control  of  the  forward  channel. 
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C.  RELATED  WORK 


Most  of  the  related  work  in  DS-CDMA  cellular  performance  has  been  focused  on 
the  reverse  channel  as  in  [4]  through  [6].  In  some  instances  the  results  of  the  reverse 
channel  analysis  can  be  extended  to  the  forward  channel.  While  these  investigations 
consider  both  shadowing  and  fading  effects,  they  do  not  incorporate  forward  error 
correction  in  the  analysis.  In  publications  where  forward  error  correction  with  soft 
decision  decoding  has  been  considered  as  in  [7]  through  [10],  the  combined  effects  of 
shadowing  and  fading  are  not  included  in  the  analysis. 

The  performance  of  the  forward  DS-CDMA  channel  with  power  control  in  a 
Rayleigh  fading  channel  was  investigated  in  [11]  without  considering  the  effects  of 
shadowing  or  forward  error  correction.  Furthermore,  [12]  considers  the  effects  of 
lognormal  shadowing  and  Rayleigh  fading  separately  and  includes  forward  error 
correction  in  the  form  of  Golay  codes  with  hard  decision  decoding.  In  [13],  Nakagami 
fading  channels  were  considered  using  a  weak  Golay  Code  with  hard  decision  decoding 
and  the  effects  of  shadowing  were  not  included. 

In  summary,  previous  performance  analyses  of  the  DS-CDMA  cellular  system 
typically  have  considered  only  a  subset  of  the  channel  effects  when  incorporating  the 
benefits  of  forward  error  correction  with  soft  decision  decoding.  Alternatively,  the 
analyses  have  considered  the  combined  effects  of  fading  and  shadowing,  but  have  not 
incorporated  forward  error  correction. 
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D.  THESIS  OUTLINE 


In  Chapter  II,  we  examine  the  properties  and  generation  of  a  set  of  Walsh 
functions,  which  are  at  the  heart  of  DS-CDMA  techniques.  We  prove  the  closure 
property  for  a  set  of  Walsh  functions.  We  introduce  the  idea  of  extending  the 
orthogonality  property  for  a  reduced  set  Walsh  functions  so  that  product  of  any  three 
Walsh  functions  from  the  reduced  set  is  identically  zero  when  integrated  over  the  bit 
period.  We  conclude  the  chapter  by  developing  a  set  of  autocorrelation  functions  for  the 
set  of  Walsh  functions  of  order  N,  and  by  describing  the  average  autocorrelation  function 
for  the  entire  set. 

In  Chapter  DI,  we  develop  the  forward  channel  model  for  the  DS-CDMA  cellular 
system  using  a  seven-cell  cluster.  Our  channel  model  takes  into  account  both  large-scale 
and  small-scale  propagation  effects.  We  use  the  extended  Hata  model  to  predict  path  loss 
and  incorporate  Rayleigh  fading  and  lognormal  shadowing  into  the  channel  model.  We 
develop  the  composite  signal  received  by  a  mobile  user  in  the  center  cell  that  includes  the 
signal  intended  for  the  mobile,  signals  intended  for  other  mobile  users  in  the  center  cell 
(intracell  interference),  co-channel  interference  from  the  adjacent  six  cells,  and  additive 
noise.  Finally,  we  place  our  mobile  user  at  the  comer  of  the  hexagonal  cell,  in  order  to 
develop  the  performance  for  the  worst-case. 

In  Chapter  IV,  we  develop  a  Signal-to-Noise  plus  Interference  Ratio  (SNIR)  and 
the  bit-emor  probability  for  the  DS-CDMA  system  in  a  Rayleigh-lognormal  channel  by 
using  Gaussian  approximation.  We  incorporate  Forward  Error  Correction  (FEC)  in  the 
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form  of  convolutional  encoding  with  soft-decision  decoding  into  our  Rayleigh-lognormal 
forward  channel  model  and  develop  a  revised  upper  bound  on  the  probability  of  bit  error. 
In  calculating  the  probability  of  bit  error  for  the  system  with  EEC  analytically,  we  require 
the  probability  distribution  function  for  a  sum  of  d  multiplicative  chi-square(two 
degrees  of  freedom)-lognormal  random  variables.  To  simplify  the  performance  analysis, 
we  approximate  this  sum  of  d  multiplicative  random  variables  as  a  multiplicative  chi- 
square(with  2 d  degrees  of  freedom)-lognormal  random  variable.  We  compare  our 
modeled  performance  results  with  results  that  are  simulated  using  Monte  Carlo  methods. 
We  conclude  the  chapter  by  incorporating  antenna  sectoring  into  our  model.  Lastly,  we 
analyze  the  performance  of  our  Rayleigh-lognormal  forward  channel  model  under  a  wide 
range  of  shadowing  conditions  and  for  various  numbers  of  users  per  cell. 

In  Chapter  V,  we  modify  the  performance  analysis  of  the  previous  chapter  to 
incorporate  user  distribution  within  the  cell.  That  is,  rather  than  assuming  our  mobile 
user  is  located  at  the  worst-case  position  in  the  comer  of  the  cell,  we  assume  that  the 
position  is  random  according  to  a  specified  distribution.  We  revise  the  upper  bound  on 
the  probability  of  bit  error  for  distributed  users  in  order  to  gain  a  more  realistic  view  of 
typical  performance  for  the  Rayleigh-lognormal  forward  channel  model  with  FEC.  We 
conclude  the  chapter  by  introducing  forward  channel  power  control  as  a  means  of 
enhancing  the  overall  performance  of  the  system.  We  analyze  the  performance  of  the 
system  with  power  control  under  a  wide  range  of  shadowing  conditions  and  for  various 
numbers  of  users  per  cell. 

In  Chapter  VI,  we  summarize  our  conclusions  and  provide  areas  of  further 


research. 
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II.  WALSH  FUNCTIONS 


Walsh  functions  are  at  the  heart  of  DS-CDMA  cellular  communications 
technology.  They  provide  orthogonal  cover  on  the  forward  traffic  channel  within  each 
cell  to  overcome  intracell  interference.  Each  mobile  user  in  the  cell  has  a  unique  Walsh 
function  assigned  which  encodes  his  traffic  coming  from  the  base  station  to  the  mobile 
handset.  By  applying  his  unique  Walsh  function  to  the  traffic  coming  from  the  base 
station,  the  mobile  user  despreads  and,  in  effect,  decodes  only  the  traffic  that  is  intended 
for  him.  When  his  unique  code  is  applied  to  the  remaining  intracell  traffic  in  the  channel, 
the  orthogonality  between  his  code  and  other  user  codes  zeros  out  their  interference.  A 
simplified  representation  of  this  process  is  depicted  in  Figure  2.1. 
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Figure  2.1  Direct  Sequence  Spread  Spectrum  With  Walsh  Orthogonal  Cover. 
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If  non-orthogonal  functions,  such  as  pseudorandom  noise  sequences,  were  used  to 
spread  the  signal,  the  intracell  traffic  intended  for  other  users  would  remain  spread  across 
the  channel’s  frequency  band.  This  intracell  interference  would  contribute  to  the  noise 
(and  now  interference)  floor  as  depicted  in  Figure  2.2  instead  of  being  eliminated.  In 
Chapter  3,  we  develop  the  cellular  model  and  look  in  depth  at  the  construction  of  the 
forward  channel  using  Walsh  functions  combined  with  pseudorandom  noise  sequences. 
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Figure  2.2  Direct  Sequence  Spread  Spectrum  Without  Orthogonal  Cover. 


Walsh  functions  have  many  desirable  properties,  which  we  will  examine  in  this 
chapter.  We  will  look  at  ways  that  Walsh  functions  can  be  generated,  how  the  set  of 
Walsh  functions  form  an  algebraic  group  under  multiplication,  their  autocorrelation 
functions,  and  an  extension  to  the  orthogonality  property  which  may  prove  useful  to  the 
communications  engineer. 
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A.  THE  SET  OF  WALSH  FUNCTIONS  AS  AN  ALGEBRAIC  GROUP 


Walsh  functions  were  first  conceived  by  J.  L.  Walsh  [14]  as  a  closed  set  of 
normalized  orthogonal  functions  over  the  interval  (0,1).  Since  that  time,  much  work  has 
been  done  to  determine  ways  to  generate  these  powerful  functions.  Today,  there  are  at 
least  four  methods  of  efficiently  generating  a  set  of  Walsh  functions,  and  each  method 
comes  with  its  own  enumeration  or  indexing  scheme,  as  follows: 

♦>  Using  the  Walsh  binary  index  and  symmetry  to  instantly  generate  any 
member  function,  called  Instant  Walsh  functions. 

❖  Using  a  product  of  Rademacher  functions,  which  we  describe  in  this 
section. 

❖  Using  Hadamard  matrices. 

❖  Using  generator  matrices. 

The  above  methods  are  well  described  in  [15],  and  we  will  only  outline  those  that  are 
necessary  to  make  our  points.  Using  these  innovations  to  construct  and  describe  the  set 
of  Walsh  functions  more  concisely,  we  have  gained  better  insight  into  the  algebraic  group 
properties  of  the  Walsh  functions. 

1.  Walsh  Functions  Defined 

We  define  the  set  of  orthogonal  Walsh  functions  of  order  N  as 
WN  =  {w .  (t) ;  te  (0 ,T) ,  j  =  0,1,...,  N  - 1} ,  consisting  of  N  =  2k  elements  that  are 

functions  of  time  and  that  have  the  following  properties  as  in  [15] ; 
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\Nj  (t)  takes  on  the  values  {+1,-1}  except  at  a  finite  number  of  points  of 
discontinuity;  where  it  is  defined  to  be  zero. 


❖  \Nj  (0)  =  1  for  all  j. 

❖  \Nj  (r)  has  precisely  j  sign  changes  in  the  interval  (0,7). 


i 

*  Jw;  (0W*(0*  = 


[o,  if  j-*k 
\T,  if  j  =  k 


❖  Each  function  w;.  (r)  is  either  odd  or  even  with  respect  to  the  midpoint  of 
the  interval  (7/2). 

Figure  2.3  depicts  the  four  member  functions  of  Walsh-4, 
WA  ={w0(t),w1(r),w2(r),w3(r)}.  In  this  case,  the  interval  (0,7)  is  broken  into  N=4 

sections,  with  each  section  being  TC  =  T/N  =  7/4  time  units  long.  We  use  Walsh  binary 


w0(t) 


Wj(t) 


w2(t) 


w3(t) 


w0-(0  000) 


Wi=(0  Oil) 


w2-(0  110) 


w3=(0  101 ) 


Figure  2.3  Set  of  Walsh  Functions  of  Order  Four. 
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indexing,  which  implies  that  the  functions  are  ordered  and  numbered  according  to  the 
number  of  sign  changes  (or  visually  zero  crossings)  on  the  interval.  Other  indexing 
schemes  are  related  to  the  way  in  which  the  Walsh  functions  are  generated. 

The  set  of  Walsh  functions  of  order  N  forms  an  algebraic  group  under 
multiplication  [15];  we  will  prove  this  fact  later  in  this  section.  We  define  a  group 
(G,*) ,  as  in  [16],  to  consist  of  a  set  of  elements  G,  together  with  an  operation  *  that  is 

defined  on  G  and  that  satisfies  the  following: 

❖  Closure  property.  For  all  a,b£  G,  a*b  e  G. 

❖  Associative  property.  For  all  a,b,ceG,  (a*b)*c  =  a*(b*c). 

♦>  Identity  property.  There  exists  an  element  e  e  G  such  that 
e*  x  =  x*e  =  x  for  all  x  E  G . 

❖  Inverse  property.  For  each  a  e  G ,  there  exists  an  element  a  e  G  that 
satisfies  a* a  =  a  *  a  =  e . 

In  the  case  of  Walsh  functions  we  have  a  set  of  N  elements,  or  member  functions, 
together  with  the  algebraic  operation  of  multiplication  (•)  forming  a  abelian 

(commutative)  group  (WN,*)  •  The  last  three  group  properties  are  easy  to  see  for  Walsh 

functions  of  any  order  and  for  t  e  (0,T)  as  follows: 

♦>  The  associative  property  is  satisfied,  since  we  are  dealing  with  simple 
multiplication  of  time  functions. 

❖  The  identity  element  is  clearly  w0  (r)  =  1 ,  for  Walsh  functions  of  any 
order. 
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***  Each  Walsh  function  serves  as  its  own  multiplicative  inverse,  since 
\Nj  (?)*W;.  (?)  =  1  for  ally. 

❖  The  order  of  the  group  is  N  =  2k .  Furthermore,  the  order  of  any  non¬ 
identity  element  of  the  group  is  two,  which  implies  that  any  element  of  the 
group  can  generate  a  subgroup  of  (at  most)  order  two. 

The  closure  property  of  Walsh  functions  remains  to  be  proven.  We  will  continue 
to  lay  the  groundwork  toward  that  end,  and  present  the  proof  in  Section  2.A.4.  In  fact, 
we  will  go  a  step  further  to  show  that  the  product  of  any  two  Walsh  functions  is  not  only 
a  member  of  the  set,  but  that  it  is  a  member  that  can  be  predicted  by  its  binary  index 
without  ever  looking  at  the  waveforms.  We  will  show  that  the  product  of  any  two  Walsh 
functions,  w,  (?)*w;.  (?) ,  is  a  Walsh  function  wk  (?)  that  can  be  predicted  by  converting 

the  Walsh  indices  i,  j,  k  to  their  corresponding  Walsh  binary  indices  Qn  Qj,  Qk ,  and 
adding  them  modulo  2,  i.e., 

Qk=Qi®Qj,  (2.1) 

where  ©  denotes  addition  modulo  2. 

The  Walsh  binary  indices  are  ^-tuples  where  k  =  log2  N .  For  example,  using  the 
set  of  Walsh  functions  W4  shown  in  Figure  2.3,  it  is  clear  that  Wj  (?)*w2  (?)  =  w3  (?) .  We 
could  have  predicted  this  result  without  looking  at  the  waveforms  by  applying  (2.1)  to  the 
binary  Walsh  indices  as  follows:  (2j=(0l);  <22  =  (l0);  Qx  ®  Q2  =  (ll)  =  <23,  which 

implies  that  w1  (?)«w2  (?)  =  w3  (?) .  This  is  a  powerful  property,  which  we  will  use  in  the 
following  sections  before  we  have  actually  proven  it. 
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2. 


Walsh  Sequences  Defined 


Walsh  functions  may  also  be  viewed  and  treated  as  binary  sequences  called  Walsh 
sequences.  The  set  of  Walsh  sequences  of  order  N  is  denoted  by 

W*N  ={«c;.;j  =  0, 1, N-l}  =  {i£>0JM’1,...,*rw_1}.  Each  element  w;.(f)eW„  has  a 
corresponding  binary  Walsh  sequence  w-s  6  represented  as  an  IV-tuple  that  is  a 

member  of  the  set  of  Walsh  sequences.  Since  each  of  the  N  sections  of  any  Walsh 
function  is  constant  and  continuous,  each  section  can  be  represented  separately  by  a 
binary  number  in  the  Walsh  sequence.  Typically,  each  section  of  the  Walsh  function  that 
has  the  value  +1  is  represented  by  a  binary  zero  in  the  Walsh  sequence.  The  sections  of 
the  Walsh  function  that  are  -1  are  represented  by  a  binary  1  in  the  Walsh  sequence. 
Figure  2.3  shows  the  Walsh  sequences  that  correspond  to  the  Walsh  functions  of  order 
four. 

Each  Walsh  function  is  mapped  as  described  above  to  its  corresponding  Walsh 
sequence,  such  that  w;.  (r)  Wj ,  or  =  (w;.(?))  M .  This  notation  means  that  the 

map  or  function  JM  acting  on  w;.  (r)  produces  the  corresponding  Walsh  sequence  wr 
Accordingly,  the  mapping,  JM :  WN  — >  W*N  takes  the  entire  set  of  Walsh  functions  to 

their  corresponding  Walsh  sequences  in  a  one-to-one  and  onto  fashion. 

The  set  of  Walsh  sequences  of  order  N  forms  a  group  under  modulo  2  addition; 
we  denote  this  group  by  (W^,,©).  Once  again,  the  last  three  group  properties  are 

obvious,  while  the  closure  property  is  still  elusive.  Just  as  the  product  of  two  Walsh 
functions  can  be  predicted  using  (2.1)  with  the  Walsh  binary  index,  we  can  predict  the 
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modulo  2  addition  of  two  Walsh  sequences,  using  the  same  binary  index  and  (2.1).  This 
attribute  guarantees  that  the  operation  is  preserved  under  the  mapping  as  follows: 

(w,(r)»wy(f))*flf  =  (w,{f))jf  ®  (2.2) 

To  prove  that  the  function  is  preserved  as  defined  in  (2.2),  we  use  (2.1),  as  follows: 

=  wk 

=  Wi  ©  Wj 

=  (w,(0)^f®(w,(0)4f. 

Accordingly,  we  see  that  our  mapping,  JM ,  is  one-to-one  and  onto  and  preserves 
the  operation.  Therefore,  the  group  of  Walsh  functions  under  multiplication  is 

isomorphic  to  the  group  (|fF*A,,©^  of  Walsh  sequences  under  modulo  two  addition  as 

defined  in  [16].  We  will  use  Walsh  sequences  to  represent  the  Walsh  functions  for  clarity 
in  some  instances,  without  any  loss  of  generality. 

3.  Generating  Walsh  Functions  from  Rademacher  Functions 

One  of  the  ways  that  Walsh  functions  can  be  generated  is  as  a  product  of 
Rademacher  functions.  For  any  set  WN  of  Walsh  functions  of  order  N  =  2k ,  there  is  a 
corresponding  set  of  Rademacher  functions  of  order  1  +  log2  N  =  \  +  k,  whose  unique 
products  completely  form  the  set  WN  of  Walsh  functions.  We  will  define  a  set 

(?);  t e  (0,r),  n  =  0,1, .../£■}  of  orthogonal  Rademacher  functions  of  order  l+k , 
as  adapted  from  [16]  and  [17],  where 
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^(r)  =  ^(2n_1r)  n  =  2,3,...,k. 

Each  Rademacher  function  ( t ) ,  where  n  >  1 ,  is  made  up  of  a  single-cycle  sine¬ 
like  square  wave  with  a  period  of  that  is  repeated  2" 1  times  over  the  interval 

(0,7).  From  the  preceding  definition,  we  see  that 

^M-^2"')-  (2'3) 

since 

^™(0=^(2"""‘') 

and  if  we  let  r=2mt ,  then 

*.«(0='3(2"“T) 

-*.w 

=  ^.(2"0- 

Every  set  of  Rademacher  functions  under  multiplication  has  an  identity 
element,  /?0  (r) .  Each  element  serves  as  its  own  multiplicative  inverse,  although  a  set  of 

Rademacher  functions  does  not  form  an  algebraic  group  under  multiplication,  since  it  is 
not  closed.  Specifically,  the  product  of  any  two  Rademacher  functions  is  not  likely  to  be 
contained  in  the  set  ^ .  Furthermore,  the  product  of  any  number  of  non-identity 
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Rademacher .  functions  will  not  be  an  element  of  ,  unless  they  are  all  the  same 
Rademacher  function. 

Figure  2.4  shows  the  set  of  Rademacher  functions,  which  are  used  to 
construct  Walsh  functions  of  order  eight.  In  order  to  build  the  Walsh  functions,  we  first 


Figure  2.4  Rademacher  Functions  for  k  =3. 


convert  the  Walsh  binary  index,  Qt ,  to  a  Gray  code  index  [15],  denoted  G, .  Both  types 
of  indices  are  ^-tuples,  where  k  =  log2  N ,  and  are  expanded  as  <2,  = 

G;  =  (gii’£;2> •••£;*)  •  The  conversion  from  <2,  to  Gi  is  accomplished  using  the 
following: 


8 a  ~  <h> 

8  ij  ~  tfi  j- 1  ®  tfjj  • 


(2.4) 
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For  example,  if  we  wanted  to  find  the  Gray  Code  index  for  w6  (?)  €  Wg ,  we  look 
at  the  Walsh  binary  index,  Q6  =  (q6iq62q6^)  =  (1 10) .  Converting  bit-by-bit,  we  find  that 


&6,i  ~  q&i — i* 

§6,2  =  4,1  ®  4,2  =  1  ©  1  —  0, 

$6,3  =4,2  0  4,3  =  100  =  L 

Accordingly,  G6  =  (101).  Table  2.1  gives  the  conversions  from  Qt  to  G,  for  all 
\Nj  (?)  e  Wg ,  where  k  =  3. 


Qi  = 

% 

4,2 

4,3 

<?,= 

8i, i 

8i, 2 

8 1,3 

Qo  = 

0 

0 

l 

G0  = 

0 

0 

0 

Qi  = 

0 

0 

i 

G:  = 

0 

0 

1 

Qi  = 

0 

1 

0 

g2  = 

0 

1 

1 

03  = 

0 

1 

1 

g3  = 

0 

1 

0 

04  = 

1 

0 

o 

= 

1 

1 

0 

e5  = 

1 

0 

1 

G5  = 

1 

1 

1 

^6  = 

1 

1 

0 

G6  = 

1 

0 

1 

Qi  = 

1 

1 

1 

G7  = 

1 

0 

0 

Table  2.1  Gray  Code  Index  Conversion  for  k  -  3. 


Using  the  Gray  code  index,  we  can  now  construct  a  set  Walsh  functions  of  any 
order  N  =  2*,  using  a  set  of  Rademacher  functions  of  order  k+\,  by  applying  the 
following: 


r  k 


i  =  0. 


(2.5) 
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In  forming  the  product  in  (2.5)  for  the  i  *  0  case,  we  only  select  the  Rademacher 
functions  that  correspond  to  the  non-zero  elements  of  the  Gray  code  index.  The 
construction  process  defined  in  (2.5)  is  accomplished  by  associating  each  element  of  the 
Gray  code  index  with  a  non-identity  element  of  as  follows,  gtJ  — >  /?M_.  (?) .  If  you 

consider  the  Gray  code  index  Gt  to  be  in  ascending  order  according  toy  =  1,  2, 

then  the  associated  Rademacher  functions  are  paired  in  descending  order  as 
k  - 1,  ...,  2, 1, 0 .  Table  2.2  shows  the  association  for  the  case  of  k  =  3 . 


Q  -  Six  Si ,2  Si ,3 

Si,2 

Si, 3 

^k+\-  j  e  *^3 

*»(') 

-*?(») 

Table  2.2  Association  of  Gray  Code  Index  with  Rademacher  Function  for  k=  3. 

Using  the  previous  example  of  w6(f)e  Ws,  we  know  that  G6=(101). 
Accordingly, 

w6  (0  =  ^3  ('MM 

=/^vP3. 

Figure  2.5  shows  the  two  Rademacher  waveforms  and  the  resulting  Walsh  function 
waveform.  It  is  clear  that  the  product  of  the  Rademacher  functions,  when 

multiplied  section  by  section,  is  the  resulting  waveform  labeled  as  W6  (?) .  We  see  that 

this  waveform  also  has  six  zero  crossings,  which  defines  w6  (?) . 
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We  have  described  how  to  go  about  generating  Walsh  functions  from  a  set  of 
Rademacher  functions,  but  a  more  powerful  fact  is  that  any  finite  product  of  integer 
powers  of  Rademacher  functions  results  in  a  Walsh  function  [17].  We  can  select  any 
combination  (with  replacement)  of  Rademacher  functions  from  the  set  ,  and  their 
product  will  still  be  a  Walsh  function  from  the  set  WN  of  order  N  =  2k .  This  is  clear 
when  you  consider  that  even  powers  of  Rademacher  functions  are  unity,  and  odd  powers 
of  Rademacher  functions  are  the  equal  to  the  original  function: 

/P„  fori  odd  _  (26) 

/?0  =  1,  for  i  even 
since  for  i  =  2k +  1  (odd  for  k=  0, 1,  2  ....)  , 
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(*.)' = (*.)"*' = = M**)'  • 

but 

(^,(t))2 = (m  '"<))2 = km)2 = /?(<) = -?„ =1, 

SO 

K)' = M*’)1 = w = *.. 

and  for  i  =  2k  (even  for  k=  0, 1, 2  . . ..) , 

(^)‘=(^)2*=(^2)‘=K)*=^=i- 


Furthermore,  when  considering  a  set  of  Rademacher  functions,  J^,  there  are 

exactly  N  =  2*  unique  products  as  defined  in  (2.5),  which  correspond  directly  to  the  N 
Walsh  functions  that  we  seek  to  build.  Table  2.3  shows  the  correspondence  between  the 


Walsh  function 

Walsh  index 

Gray  Code  Index 

Products 

W,M 

Q  —{8iv8i2*“m*8ik} 

k 

f  -A -  ^ 

k 

W»M 

(00... 00) 

(00...  ooj 

w,(») 

(00...01) 

(00...01) 

w  i(>) 

(00...10) 

(00...  11) 

W  ,(<) 

(00...11) 

(00...10) 

... 

... 

... 

... 

W  N-i(t) 

(11-11) 

(10.. .00) 

*> 

Table  2.3  Walsh  Functions  from  Rademacher  Functions. 
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Walsh  functions  of  order  N  and  the  products  of  Rademacher  functions.  As  shown  above, 
other  possible  combinations  of  Rademacher  functions  that  repeat  elements  are  not  unique 
(e.g.,  since>  =  !•) 

The  last  three  sections  have  been  concerned  with  defining,  describing  and 
motivating  the  Walsh  functions  and  their  component  Rademacher  functions.  With  this 
foundation,  we  can  now  address  closure  as  a  final  step  in  proving  the  algebraic  group 
properties.  We  will  prove  that  the  product  of  any  two  Walsh  functions  is  not  only 
contained  in  the  set,  but  is  a  predictable  member  function  as  described  by  (2.1). 

4.  The  Closure  Property  of  Walsh  Functions 

The  set  of  Walsh  functions  of  order  N  forms  an  algebraic  group  under 
multiplication  provided  that  the  properties  of  closure,  associativity,  identity,  and  inverse, 
as  defined  in  Section  II.A.1,  are  satisfied.  By  defining  the  individual  Walsh  functions 
w;.  (?)e  WN  in  terms  of  their  component  Rademacher  functions  by  (2.5),  we  can  prove 

that  the  set  is  closed.  That  is,  for  any  w,  (?) ,  w;.  (?)  6  WN ,  w,.  (?)»w;  (t)eWN. 

Clearly,  if  either  of  the  Walsh  functions  w ;(?)  or  w;  (?)  are  the  identity,  w 0(?), 
then  closure  holds.  Accordingly,  we  let 

w.(r)=  tl^-M 

m- 1 
£«>*  0 

where  w,.  (r)  *  w0  (?)  has  a  corresponding  Gray  Code  index  of  Gf  =  (&pgi2>— >&*) » and 

p= i 
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where  wy  (r)^  w0(r)  has  a  corresponding  Gray  Code  index  of  Gj  =  (gn,gj2 . gJk )  ■ 


Then  their  product  is 


w<  (<)•»;  w  -  n  (o-  ii  (0- 

m=l  /?=1 

£(>*0 


When  we  combine  the  factors,  any  Rademacher  function,  /^+1_n,  that  is  in  one  but  not 


both  of  the  two  products  will  survive:  (gin=l  a  ghn  =  0)  v  ( gin  =  0  a  gjn  =  1  ] 


implies  that  /^+1_n  survives  the  multiplication.  If  a  Rademacher  function  is  in  both 


product  terms,  it  will  not  be  in  the  final  product  term  since  =  /?0  =  1 .  Accordingly, 


(glVt  =  1  a  ghn  =  l)  implies  that  /^+1_w  will  not  be  in  the  final  product.  Clearly,  if  a 


Rademacher  function  is  in  neither  product  term,  =0  a  gj/t  =  0),it  will  not  be  in 


the  final  product. 


Accordingly,  the  combination  of  the  two  products  can  be  rewritten  as  one  product 


using  a  new  Gray  code  index,  call  it  Gr  =  (grVgr2,-.,gri) ,  where  the  gr/l  =1  if  /^+1_„ 


survives  the  multiplication  of  products,  and  grfl  =  0  if  k?M_n  does  not  survive,  or  was 


never  a  part  of  the  one  of  the  products.  From  the  description  above,  we  know  that 


Sr,n  =  gi„  ®  gj/,  f°r  n=l>  2, ...,  k.  We  continue  with  w,.  (r)*wy  (r)  as  follows: 


(')•«;  (»)  =  ri  (<)•  n  m 

m= 1  p- 1 

-ft 

n=\ 

Sr..*  0 
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where  grn  =  gin  ©  gjn .  We  recognize  this  final  form  to  be  that  of  the  Walsh  function 
wr  (?) ,  whose  Gray  Code  index  is  determined  by 

Gr  =  (grv8r2’-’8ri)  =  (&1  @  8 jvSil  ©  8 . 8ik  ®  8 jk)  =  Gi  ®  Gj~ 

Since  w,.  (?)*w;  (?)  =  wr  (?)  e  Wn  ,  we  know  that  the  set  of  Walsh  functions,  WN , 

is  closed  under  multiplication.  Furthermore,  since  Gray  coding  as  defined  by  (2.4)  is  a 
one-to-one  and  onto  mapping,  we  see  that 

Gr  =  G,  ©  Gj  <r^Qr=  Qi  ©  Qj, 

which  proves  (2.1). 

B.  WALSH  ORTHOGONALITY  EXTENDED 


The  most  significant  property  of  the  Walsh  functions  as  applied  to  the  DS-CDMA 
cellular  problem  is  that  of  orthogonality  between  two  different  Walsh  functions.  That  is 
for  any  w;  (?) ,  w*  (?)  e  WN  , 


if;  *  k 
if;  =  k 


In  some  applications,  the  integral  of  a  product  of  three  Walsh  functions  is  of  interest.  Just 

T 

what  is  the  value  of  the  integral,  J  w;  (?)  w;  (?)  w*  (t)dt?  The  answer  depends  on  the 

0 

relationship  between  the  three  Walsh  functions  in  the  integral.  For  example,  if 

w.  M  w;  (0  =  w*  (0  ’ that  is  2,  ®  Qj  =  Qk » then  we  know that 
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T  T 

J  w«  W  w;  (0  w*  (0  dt  =f  w*  (r)  Wt  (f )  dt  =  r. 

0  0 

On  the  other  hand,  by  closure  we  know  that  if  wf  (t)  w;.  (/)  *  w*  (?) ,  it  is  equal  to  some 
Walsh  function,  call  it  w;  (?) ,  where  l  #  k .  Accordingly, 

T  T 

J  w<  M  w;  M  w*  W  =J  w;  (0  w*  (7) dt  =0-  (2-7) 

0  0 

The  fact  that  the  integral  value  depends  on  the  relationship  of  the  Walsh  functions  is  a 
troublesome  result  in  the  analysis  of  cellular  BER  performance  because  it  requires  a  high 
degree  of  specificity,  when  the  result  needs  to  be  applicable  to  the  general  case. 
Furthermore,  when  applying  Walsh  functions  to  eliminate  intracell  interference,  the 
product  of  three  Walsh  functions  integrated  over  the  bit  period  needs  to  be  (predictably) 
zero. 

One  way  around  this  problem  is  to  eliminate  some  of  the  Walsh  functions  from 
the  complete  set  WN ,  and  force  (2.7),  which  we  use  to  define  extended  orthogonality,  to 

be  true  for  every  Walsh  function  in  the  reduced  set  WN .  The  downside  to  this  solution  is 

that  by  eliminating  Walsh  functions  from  our  useable  set,  we  reduce  the  number  of 
mobile  users  that  can  be  accommodated  in  the  cell.  This  is  because  each  user  requires  a 
unique  Walsh  function;  when  we  have  assigned  all  Walsh  functions  in  the  useable  set,  no 
new  users  may  join.  Accordingly,  we  need  to  ensure  that  (2.7)  is  enforced  while 
eliminating  the  minimum  number  of  Walsh  functions  and  thereby  maximizing  the  size  of 
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our  new  set,  WN .  Furthermore,  once  we  determine  the  size  of  the  maximum  set,  we  will 
propose  a  method  of  generating  the  optimal  reduced  set,  WN . 

1.  Walsh  Subgroups  of  Order  Four 

We  went  to  great  lengths  in  Section  RA  to  prove  the  closure  property  for  the  set 
of  Walsh  functions  WN .  In  order  to  extend  the  orthogonality  of  the  Walsh  functions  as 

described  by  (2.7),  we  need  the  opposite  of  closure  for  our  reduced  set  WN .  That  is,  for 
all  the  (non-identity)  Walsh  functions  in  our  reduced  set  w; (r),w; (r) 6  WN ,  we  need  to 
ensure  that  their  product  is  not  in  the  reduced  set, 

wi(0*w;(r)g  (2-8) 

for  i  j  and  i,  j  *  0 . 

To  accomplish  (2.8)  we  need  only  ensure  that  there  are  no  Walsh  subgroups  of 
order  four  contained  in  the  reduced  set.  A  subgroup,  S  of  a  group  G,  is  a  subset  of  G 
which  is  a  group  itself  (as  defined  in  Section  n.A.l).  The  subgroup  is  under  the  same 
operation  as  G,  but  restricted  to  the  elements  of  S.  Any  Walsh  subgroup,  Qt?Jk,  of  order 

four  is  of  the  form 

^HWoM-W;(0>W;M’W*  (*)}• 

The  subgroup  is  required  to  have  the  multiplicative  identity  element  w0(r)  and  the 

elements  of  the  subgroup  must  satisfy  the  inverse  and  closure  properties.  Since  each 
element  serves  as  its  own  inverse,  it  is  the  closure  property  that  distinguishes  subgroups 
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of  order  four  from  simple  four-member  subsets  of  Ww .  In  order  to  satisfy  the  closure 
property  for  subgroups,  the  three  non-identity  elements  w,.  (?),w;.  (?),w*  (?),e  Q%Jk  must 
satisfy  the  following: 

w;M-w,W  =  w*W-  (2-9) 

which  implies  that  Qt  ®  <2;  =  Qk .  Accordingly,  we  know  that  the  product  of  any  two  of 
the  non-identity  elements  results  in  the  third  non-identity  element,  since  by  modulo  2 
addition,  Qi  ©  Qj  =  Qk  means  <2,  ®Qk  =  Qj  and  <2,  ®Qk  =  Qn  which  imply 

w,.  (?)w*  (?)  =  w;.  (?)  and  w;  (?)  \Nk  (?)  =  w;  (?) ,  respectively.  Furthermore,  we  know  the 
product  of  any  element,  ws  (?)  e  Q%jk  in  the  subgroup,  with  the  identity  is  itself,  namely 
W1  (?)  w0  (?)  =  Wf  (?) .  Finally,  the  product  of  any  W5  (?)  6  Ci^.  with  itself  is  the  identity 
element,  ws(?)w5(?)  =  w0(?).  Therefore,  (2.9)  guarantees  closure  for  a  subgroup 
G&?k  c  WN  of  order  four. 

Consequently,  we  see  that  any  two  non-identity  Walsh  functions  uniquely 
determine  the  final  non-identity  member  of  the  subgroup.  For  an  example,  let  us  look  at 
a  subgroup  of  order  four  from  W& .  For  ease  of  notation,  we  will  use  Walsh  sequences, 

w0,wl,tc2,tcvwi,u>5,w6,w1] .  We  know  that  the  identity  is  always  a  member  of 
the  subgroup.  If  we  select  wv  and  w6  as  the  first  two  non-identity  members,  the  third  is 
uniquely  determined  by  wx  ©  w6 ,  or  by  Q1®Q6=Q1,  as  shown  in  Table  2.4. 
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2; 

wi 

wx 

(001) 

(00001111) 

U>6 

(110) 

(01011010) 

wn 

(111) 

(01010101) 

Table  2.4  Building  the  subgroup  {w0 ,  , w6 ,  w1 } . 


If  we  had  started  with  any  of  the  two  non-identity  elements  in  the  above  subgroup 
(e.g.,  wvw1  or  w6,w1 ),  we  would  have  produced  the  same  subgroup,  since  any  sum  of 

the  two  non-identity  elements  results  in  the  third.  Continuing  this  process  for  all  the 
pairs  of  non-identity  elements  results  in  a  complete  listing  of  all  the  subgroups  of  order 

four  from  : 

{w0tt\w2w3 }  ,\wQwxwAw3  ]  ,\w^wxwew-, ) , 

{tt>0fr2*r4*r6  },{w0w2w5w2 }, 

We  will  now  consider  the  elimination  of  Walsh  functions  in  order  to  break  all 
subgroups  of  order  four  from  the  full  set,  Wv . 

2.  A  Reduced  Set  of  Walsh  Functions 

When  forming  our  reduced  set  WN  of  Walsh  functions  from  the  full  set  WN ,  we 
want  to  break  all  subgroups  of  order  four  contained  in  WN  in  order  to  satisfy  (2.8).  Our 
resulting  set  will  no  longer  be  an  algebraic  group  but  will  extend  the  orthogonality 
property  to  satisfy  (2.7).  So  exactly  how  many  subgroups  of  order  four  are  there  in  WN  ? 
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We  can  see  from  Section  n.B.l  that  for  W8,  there  are  seven  subgroups  of  order 

four.  We  built  those  subgroups  using  all  possible  combinations  of  two  non-identity 
elements,  and  eliminating  redundant  subgroups.  The  number  of  different  combinations 


of  size  k  from  a  set  n  is  the  binomial  coefficient  defined  by 


(  n ^ 


KkJ 


n\ 


k\{n-k)\ 

Since  there  are  N-l  non-identity  elements  for  WN ,  and  we  are  generating  combinations 


of  two  elements,  the  number  of  possible  subgroups  is  based  the  binomial  coefficient, 


(jV-lWtf-l)(iV-2) 

l  2  J=  2 


(2.10) 


The  binomial  coefficient  (2.10)  determines  how  many  different  pairs  of  non-identity 
elements  are  contained  in  the  set  WN ,  or  ~W*N .  For  any  single  subgroup  of  order  four, 

say  Ofc  ={w^Vjtck},  we  know  there  are  three  pairs  of  non-identity  elements,  namely 
WjWj,  u>jWk,  and  w,wk .  In  effect,  each  subgroup  is  counted  three  times  by  (2.10).  We 
scale  (2.10)  accordingly  and  determine  that  for  any  set  of  Walsh  functions,  WN ,  order  N, 
there  are  SN  subgroups  of  order  four,  which  is  defined  by 


Sn~3 


f  N  -1^ 


(2.11) 


V  “  ) 

Applying  (2.11)  to  ,  we  calculate  the  number  of  subgroups  of  order  four  to  be 
1*6 

Sz  =  =  7 ,  which  is  how  many  we  found. 
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Now  that  we  know  how  many  subgroups  of  order  four  are  associated  with  WN , 
namely  SN ,  we  can  start  eliminating  Walsh  functions  from  WN  to  form  the  reduced  set 
WN .  To  demonstrate,  we  fall  back  to  the  subgroups  of  Tf^.  In  order  to  simplify  the 

notation  of  the  subgroups,  we  will  use  a  tri-graph  to  represent  each  subgroup  of  order 
four.  The  tri-graph  will  reflect  the  Walsh  index  of  the  three  non-identity  members  as 
shown  for  W\ : 

jk-wd  {wow\waw5\  K w\w6wi) 

1-2  —  3  1  —  4--5  1-6-7 


{w0w2w4w6}  {w0w2w5u,2}  {w0w3«,4 K.?}  {k>0«»3«>5k.6} 
2-4-6 


"2-T-7 


W3WAW1. 
3-4-7 


3-5-6 


We  will  now  begin  to  eliminate  Walsh  functions  from  Ws,  in  order  to  form  one 


possible  reduced  set  Ws.  As  we  will  see,  the  WN  are  not  unique.  Further,  the  number  of 
elements  that  must  be  eliminated  to  break  all  SN  subgroups  varies  depending  upon  how 
the  elements  eliminated  are  chosen,  and  how  that  choice  relates  to  past  choices.  For 
example,  the  elimination  of  wx  from  W\  results  in  the  following  subgroups  being 

broken : 


and  leaves  unbroken: 


1-2-3, 1-4-5,  and  1-6-7 


2-4-6, 2-5-7,  3-4-7,  and  3-5-6 

The  elimination  of  w2  results  in  the  following  additional  subgroups  being  broken: 


2-4-6  and  2-5-7 
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and  leaves  unbroken: 


3-4-7  and  3-5-6 

The  elimination  of  wz  results  in  the  final  subgroups  being  broken: 

3-4-7  and  3-5-6 

By  eliminating  wv  w2,  and  w3  from  ,  we  effectively  created  a  reduced  set  of  Walsh 

functions  Ws  ,  which  has  none  of  the  subgroups  of  order  four 

remaining  from  W\ ,  and  accordingly  satisfies  (2.7).  In  this  instance  we  only  had  to 

eliminate  three  elements  to  break  all  subgroups,  leaving  five  useable  Walsh  functions 
with  extended  orthogonality.  However,  eliminating  a  different  set  of  three  elements  may 
or  may  not  break  all  the  subgroups. 

Let  us  begin  again  with  W\ ,  and  this  time  eliminate  w5,  w6,  and  tc7 .  The 
elimination  of  w5  results  in  the  following  subgroups  being  broken  : 


1-4-5,  2-5-7,  and  3-5-6 


and  leaves  unbroken: 


1-2-3, 1-6-7,  2-4-6,  and  3-4-7 

The  elimination  of  w6  results  in  the  following  additional  subgroups  being  broken: 

1-6-7  and  2-4-6 

and  leaves  unbroken: 


1-2-3  and  3-4-7 

The  elimination  of  wn  results  in  the  subgroup  3-4-7  being  broken,  but  leaves  the 

subgroup  1-2-3  unbroken.  In  order  to  break  all  subgroups  we  must  select  one  more 

element  to  eliminate  from  the  yet-unbroken  subgroup.  If  we  choose  to  eliminate  wz ,  we 
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end  up  with  a  different  reduced  set,  Ws  =  {w0,wvw2,wA)  which  satisfies  (2.7),  but  has 
one  fewer  elements.  It  is  clear  that  we  need  to  determine  the  minimum  number  of 
eliminations  necessary  to  break  all  subgroups,  thereby  generating  a  maximal  set  WN ,  and 
develop  a  method  to  generate  that  maximal  set  consistently. 

3.  Maximizing  the  Reduced  Set  WN 

We  saw  that  three  channel  eliminations,  {w{w2w2) ,  could  produce  the  desired 
reduced  set  in  the  case  of  W\.  On  the  other  hand,  any  three  channel  eliminations, 
((*05*r6*r7),  for  example)  did  not  produce  the  desired  reduced  set,  since  they  did  not 

break  all  the  subgroups  of  order  four  in  W\ .  The  difference  between  these  two  scenarios 

is  that  the  choices  in  the  first  case  maximize  the  number  of  subgroups  eliminated  with 
each  choice,  while  the  second  set  of  choices  did  not.  Accordingly,  we  need  to  determine 
how  many  subgroups  are  broken  with  each  elimination,  if  we  assume  we  can  maximize 
the  number  broken  with  each  channel  elimination. 

If  we  consider  the  elements  of  WN ,  we  can  see  that  each  non-identity  element 

\Nj  (t)eWN  is  contained  in  exactly  N/ 2-1  subgroups.  This  is  because  any  single 

element  w;  (t)  can  combine  with  the  remaining  N-2  elements  in  only  N-2  combinations 

w,  (t)vjj  (r)  wherein  j  and  i,j*0  as  follows: 
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w,(t)w,(t)  ' 
*,(>)*.(') 


w,W 


(') 


\N-2. 


w<W 


w*W 


Moreover,  any  one  subgroup  is  counted  twice  as  w,.  (r)w;.(r)  and  wi(?)wjt(r). 

Accordingly,  we  see  that  w(.  (?)  is  contained  in  exactly  l/2(N -2)  =  A/2-1  subgroups. 

For  at  least  the  first  channel  elimination,  N/ 2-1  subgroups  are  broken.  For  subsequent 
channel  eliminations,  the  maximum  number  of  subgroups  broken  depends  on  how  many 
channel  eliminations  have  already  taken  place.  In  order  to  better  describe  the  channel 
elimination  process,  let 

<?,  =  the  Walsh  channel  removed  on  the  zth  elimination, 

%  =  the  set  of  subgroups  which  contain  the  eliminated  channel  e( , 

4  =  the  number  additional  subgroups  broken  by  the  zth  channel 
elimination, 

the  total  number  of  unique  subgroups  broken  by  the  first  z 

i 

channel  eliminations,  i.e., 

j=l 

%  =  the  set  of  subgroups  yet  unbroken  after  the  zth  elimination,  and 
[S]  =  the  number  of  elements  in  the  set  S. 


We  know  that  when  the  number  of  broken  subgroups  equals  the  total  number  of 
subgroups,  i.e.,  ^  =  SN ,  all  subgroups  have  been  broken,  and  we  have  a  valid  reduced 

set,  WN .  Let  us  use  once  again  to  illustrate  the  above  notation,  and  demonstrate 
how  to  maximize  the  number  of  subgroups  broken  with  each  channel  elimination. 


Choosing  e1  =  w5  results  in  the  following  : 
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§[  =  {1-4-5, 2-5-7,  3-5-6}  (all  the  subgroups  containing  ws), 

*i  =3, 

=  *!  =  3,  and 

^  ={1-2-3, 1-6-7, 2-4-6,  3-4-7}. 

For  the  first  elimination  we  know  that  *x  =  [if],  since  all  of  the  subgroups  in  ^  are  newly 
broken.  We  can  see  that  the  set  of  unbroken  subgroups,  S' ,  has  four  elements.  We  want 
to  choose  the  next  channel  elimination  such  that  the  maximum  number  of  subgroups  are 
broken.  We  see  that  for  any  choice  e2  =  wj  e  {wlw2wzwAw(w1}  exactly  two  subgroups 

will  be  broken. 

Choosing  e2  =  w6  results  in 

%  =  {l -6-7, 2-4-6,  >5^f}  previously  eliminated , 

*2  =  2, 

=  *i  +  *2  =  5,  and 
=  {1-2-3,  3-4-7}. 

Since  we  already  broke  the  subgroup  3-5-6,  by  eliminating  w5 ,  it  is  not  counted  in  *2 , 
the  number  of  subgroups  broken  by  e2 .  Examination  of  ^  reveals  that,  in  order  to  break 
the  maximum  number  of  remaining  subgroups,  we  must  next  choose  e~  =  w3 .  This 
choice  breaks  the  remaining  two  subgroups.  Any  other  selection  for  e3  breaks  only  one 
subgroup.  Accordingly,  choose  e3  =  w2 ,  which  results  in 

5=  {1-2-3,  3-4-7,  >5^}, 

*3  =  *2  =  2, 

^  =  ^  +  *2  +  *3  =  *x  +  2*2  =  7,  and 
%  =  0. 
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We  see  that  for  the  case  of  W\ ,  eliminating  a  minimum  of  N/2-1  =  3  channels 
broke  all  Ss  =  7  subgroups.  If  we  had  not  optimized  our  choices,  it  could  have  taken 
more  than  N/2-1  eliminations,  but  we  can  find  no  way  to  eliminate  all  subgroups  with 
less  than  N/2-1  eliminations.  We  have  found  that  this  result  holds  true  in  general  for 
WN ,  as  we  will  sketch  below. 

For  the  general  case,  WN ,  we  assume  that  with  each  channel  elimination  we  break 
the  maximum  possible  number  of  subgroups  as  shown  for  W\ .  Furthermore,  we 
observe  the  following  properties: 

[if]  =  N/ 2-1  =  4  for  all  i,  and 

| if  n  |  =  1,  since  et  and  e;.  are  both  contained  in  at  most  one  subgroup. 

For  i=l,  eliminate  any  e1.  Since  ^contains  N/2-1  elements  which  represent  distinct 
subgroups,  then  we  know  that  =  |^J  =  N/2-1  subgroups  are  broken  by  eliminating 

e\  • 

For  i= 2,  eliminate  any  e2.  The  set  ^  also  contains  N/2-1  elements.  However, 
there  is  one  element  in  ^  that  was  contained  in  If ,  since  one  of  the  subgroups  must  have 
both  ex  and <?2 .  Accordingly,  we  know  that  ^2=|^]-[^n^]  =  (iV/2-l)-l  =  lV/2-2 

additional  subgroups  have  been  broken.  A  total  of  ^  =  &x  +  &2  subgroups  have  been 
broken  after  the  second  elimination. 

For  i= 3,  eliminate  e3  such  that  |^n^n^]  =  l  in  order  to  maximize  the  number 
of  broken  subgroups.  Accordingly,  we  know  that 
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and 


^=^+^+^3  =  ^4-  2^2 . 

For  i=4,  eliminate  e4  such  that  n 3  n $]  =  n g  n 5]  =  [$  n  g  n $]  =  0 

in  order  to  maximize  the  number  of  broken  subgroups.  Accordingly,  we  know  that 

and 

^=4+2^  +  ^. 

If  we  carry  on  in  this  manner,  maximizing  the  number  of  broken  subgroups  with 
each  channel  elimination,  we  find  that  after  i  =  N/ 2-1  channel  eliminations,  every 
subgroup  in  WN  is  broken.  The  final  ^ ,  which  we  will  call  £3S ,  at  channel  elimination 

/  =  i  =  N/2  - 1 ,  is  of  the  form 

—  6- j  +2^  +  4^4  +  8#g  + 16^16  +  ...  +  -^-^i,  (2.12) 

4  4 

=  2^,  j  =  2°,2\22,...,2k-2=i, 

j 

where  £j=N/2-j. 

The  pattern  represented  in  (2.12)  can  be  observed  in  Table  2.5  for  W16 . 
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1-2-3 

1-4-5 

1-6-7 

1-8-9 

1-10-11 

1-12-13 

1-14-15 

2-4-6 

2-5-7 

2-8-10 

2-9-11 

2-12-14 

2-13-15 

3-4-7 

3-5-6 

3-8-11 

3-9-10 

3-12-15 

3-13-14 

4-8-12 

4-9-13 

4-10-14 

4-11-15 

5-8-13 

5-9-12 

5-10-15 

5-11-14 

6-8-14 

6-10-12 

6-11-13 

7-9-14 

7-10-13 

7-11-12 

Table  2.5  The  Set  of  Subgroups  of  Order  Four  for  W16 . 


In  the  table,  we  have  organized  the  entire  set  of  subgroups  for  Wl6  according  to 
the  smallest  non-identity  Walsh  index  in  the  subgroup.  There  are  exactly  J[6  =  =  35 

subgroups  of  order  four  in  W16 .  If  we  let  el  =wi  for  z  =  l,2,...,Af/2-l,  we  successively 

eliminate  each  row  in  the  Table  2.5,  which  means  we  have  broken  all  the  subgroups  in 
the  zth  row  with  the  zth  elimination.  The  number  of  subgroups  broken  with  the  z'th 
elimination,  ,  is  the  number  of  elements  in  the  zth  row.  For  W16 , 

4  =  7  =  N/  _  i 

42=<Z3=6  =  %-2, 

^4=^=^=^7=4  =  i%-4, 

2&2  +  4^  =  35  =  S16. 

Accordingly,  after  z'  =  iV/2-l  =  7  channel  eliminations,  all  S16  =  35  subgroups  are 
broken. 

We  now  show  that  for  a  set  of  Walsh  functions,  WN ,  of  any  order  N, 
i  =  N/ 2-1  channel  eliminations  break  all  SN  subgroups.  Using  (2.12),  we  know  that 
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&f=il  +  2d2  +  4^4  +  8^g  + 16^16  + ... +-f-4 

4  4 


• _  /^O  M  _  JL 

J  —  £  ,A  ,Z,  4? 


where  £j=N/2-  j .  If  we  let 


7  =  2”,  and 


then 


•»/=  E^=i2'«.-S2’(J%-2")- 

j  n=0  n= 0 

,  *-2  *-2 
=  N/2jJ2n-jJ2ln  , 


(geometric  series) 


=  AI/ 


/t=0  n=0 


/ 1  o*-i 


1-2' 


V  1-2  2 


'2 

N2  -  3N  -2 


^  1  _  4*_1  ^ 


=  %(2--l)- 


^(2*-1)2-l^ 


“X 


(A-l)(A-2) 


2)Z 

-  V 

(N-l'] 

J 

~/3 

to 

- 


We  have  found  that  N/ 2-1  channel  eliminations  is  the  minimum  number 
required  to  break  all  subgroups  of  order  four  contained  in  WN .  Accordingly,  we  know 
that  the  remaining  N/ 2  +  1  elements  from  Ww  will  make  up  the  maximum  set  WN.  That 
is,  any  maximum  set  WN  will  have  exactly  A/2  +  1  member  functions.  In  the  next 
section,  we  will  propose  a  method  for  generating  this  maximum  set  WN . 
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4.  Generating  a  Maximum  Set  WN 

We  saw  in  Section  HB.3  that,  by  eliminating  the  first  N/ 2-1  channels  from  WN , 

we  can  break  all  the  subgroups  of  order  four  from  WN  and  generate  a  maximum  set,  WN , 
which  has  extended  orthogonality  as  defined  by  (2.7)  for  all  elements  in  the  set.  By 
eliminating  the  first  N/2  - 1  channels  our  reduced  set  looks  like 

WN  =  {w0  (t),wf  (t),  (t) . w^j  (f)}. 

For  example, 

^6={WoW-Ws(0>W9W’WloW’WllW>W12W’W13W’W14W-W15W]- 
While  this  method  of  generating  WN  certainly  works,  it  is  rather  rigid.  We  need  a 
method  that  is  more  adaptable  to  various  uses.  For  example,  we  may  want  to  devise  a 
valid  reduced  set  WN  which  includes  specific  Walsh  functions  that  are  already  in  use,  or 

permanently  assigned  for  various  control  functions.  We  may  want  to  use  the  reduced  set 
concept  when  assigning  Walsh  functions  from  the  full  set,  as  a  best  effort  design.  That  is, 
to  the  maximum  extent  possible,  the  Walsh  functions  in  use,  or  assigned,  are  from  a  valid 
reduced  set,  which  can  be  adapted  as  channels  (functions)  are  added  or  dropped,  until 
user  capacity  exceeds  that  which  can  be  supported  by  the  reduced  set. 

There  are  many  possible  combinations  of  N/2  +  1  Walsh  functions  which  form  a 

valid  reduced  set  WN .  We  submit  the  following  more  flexible  method  of  generating  the 

set  WN: 
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❖  Start  by  selecting  non-identity  channels  from  the  pool  of  channels 
from  WN  to  include  in  WN.  The  identity,  w0(?) ,  will  always  a 

member  of  the  maximum  set,  and  therefore  will  not  be  considered 
during  our  selection  process. 

❖  Each  new  channel  selected  dictates  the  channels  that  must  be 
discarded.  At  each  step  the  newly  selected  channel  is  multiplied 

by  each  of  the  previously  selected  members  of  WN  to  determine 

which  elements  from  the  pool  of  undesignated  channels  must  be 
discarded. 

❖  After  exactly  N/ 2-1  distinct  channels  have  been  discarded 
(eliminated),  all  undesignated  channels  from  the  pool  are  included 
in  WN. 

♦>  If  the  number  of  channels  eliminated  exceeds  the  target  number 
N/ 2-1  during  step  i,  go  back  and  re-select  another  channel  to 
keep  for  step  i  using  the  repeat  final  choice  procedure  detailed 
below. 

❖  To  repeat  the  final  choice,  select  a  channel  that  is  the  product  of  a 
previously  kept  channel  with  a  previous  discarded  channel, 
ws  (?)  =  w*.  (?)  wd  (?)  g  WN  which  is  not  already  a  channel  chosen 
to  be  kept. 

For  example,  we  will  build  a  valid  reduced  set  W*l6  from  as  depicted  in  Table  2.6. 
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Step 

Pool  of  undesignated  channels 

Keep  in  W]6 

Discard 

[M] 

1 

J  WXW2W3W4W^W2W%  ] 
\w9w10wnwnw13wuwl5j 

wx 

none 

0 

2 

(w2w3w4w5w6w2ws  I 
\w9w10wuw12w13w14w15  J 

wx  ©  w3  =  w2 

1 

3 

[w4w3w^w1w%  I 

\w9w10wnwl2wx3wuw15  J 

wn 

wl@w1=  w6 
w3  ©  w2  =  w4 

3 

■ 

{  w5  W%  W9Wl0WuWl2Wl3  wx  4  IL\  5  j 

wn 

ww  w%  wl2 

6 

X 

{w5w9w13wl4w15} 

X 

X 

5 

{tv5w9wl3wuwl5} 

w5 

w4  w6  w2  wu 

7 

6 

{w9wl3wl5}-+ 

w9wnwl5 

none 

7 

Table  2.6  Building  a  Valid  Reduced  Set  W\6 . 

For  Step  1,  all  of  TT^S  is  in  the  pool  of  undesignated  channels  from  which  we  can 
choose.  We  select  wx  to  keep  in  W\6 ,  and  need  not  discard  any  channels.  In  Step  2,  we 
take  wx  out  of  the  pool  of  undesignated  channels,  and  choose  to  keep  w3 .  The  choice  of 
w3  dictates  that  w2  =  w1®  w3  must  be  discarded  or  eliminated.  At  this  point  the  total 
number  of  channels  eliminated,  |{ue,.}J,  is  1.  We  continue  in  this  manner  until  Step  5. 
In  our  first  attempt  at  Step  5,  we  selected  a  channel  to  keep  which  resulted  in  too  many 
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channel  eliminations.  We  require  [{uej]  =  N/2-l  =  7  channel  eliminations  for  W]6 

in  order  to  build  a  maximum  set  W\6 .  Our  choice  in  the  initial  Step  5  results  in  a  total  of 
10  channel  eliminations.  Accordingly,  we  ignore  that  result  and  select  w5=w1®  w4 
using  the  repeat-final-choice  procedure.  This  time  we  achieve  the  desired  total  number 
of  channel  eliminations,  and  accordingly  build  a  maximum  set 

W\6  =  { w9wxw2w$w9w9wx ! 3*015 } .  It  turns  out  that  in  Step  5  any  other  selection  from 
w5,  w9,  wn,  or  wl5 ,  other  than  wu  would  have  produced  the  desired  result  of  one  new 
channel  elimination,  namely  tvu .  That  is,  any  one  of  the  following  versions  of  Step  5 
could  be  used. 


5v3 

{w5w9w13w14w15} 

w9 

wl0  wu  w2 

1 

5v4 

{wsw9wnwuwl5} 

1 

5v5 

{w5w9wl3wl4w15} 

w15 

wu  wi2  IL\  WA 

1 

Let  us  consider  the  set-building  process  from  the  point  of  view  of  maximizing  the 
subgroups  broken  with  each  elimination  as  we  did  in  Section  n.B.3.  After  Step  4  is 
complete,  the  only  subgroups  that  had  not  yet  been  broken  were,  1-14-15,  3-13-14,  5-11- 
14,  7-9-14.  It  is  obvious  that  channel  14  needs  to  be  eliminated,  and  choosing  any  other 
element  to  keep  besides  wu  gives  us  the  desired  result. 

To  show  that  we  have  broken  all  35  subgroups  of  order  four  from  W\6 ,  we  have 
reorganized  the  subgroups  (from  Table  2.5)  around  the  channels  eliminated,  ei ,  and  in 
the  order  in  which  they  were  eliminated,  as  shown  in  Table  2.7. 
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sr 

II 

2-9-11 

2-12-14 

2-13-15 

6-8-14 

6-9-15 

6-10-12 

6-11-13 

1-4-5 

3-4-7 

4-8-12 

4-9-13 

7-10-13 

sr 

ii 

1-8-9 

3-8-11 

5-8-13 

7-8-15 

1-12-13 

3-12-15 

5-9-12 

7-11-12 

H 

oT 

1-14-15 

3-13-14 

5-11-14 

7-9-14 

Table  2.7  The  Set  of  Subgroups  of  Order  Four  for  W16 . 

The  fact  that  we  eliminated  ex  through  e1  proves  that  we  have  broken  all  subgroups  in 
Table  2.5,  one  row  at  a  time.  Accordingly  we  have  generated  a  valid  maximum  set, 
W^6,  whose  elements  have  extended  orthogonality  as  defined  by  (2.7). 

In  Section  HB  we  have  introduced  the  notion  of  extended  orthogonality 
for  Walsh  functions.  We  have  detailed  the  Walsh  function  subgroups  of  order  four,  as 
they  play  a  key  role  in  extending  the  orthogonality  property.  We  have  determined  how 
many  subgroups  of  order  four  exist,  and  how  many  channel  eliminations  are  required  to 
break  those  subgroups  for  a  set  of  Walsh  functions  of  order  N  to  generate  a  reduced  set  of 
Walsh  functions.  Finally,  we  have  proposed  an  adaptable  method  of  generating  a  valid 
reduced  set  of  Walsh  functions  that  guarantees  the  maximum  number  of  channels  from 
the  original  set. 
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C.  AUTOCORRELATION  OF  WALSH  FUNCTIONS 


The  autocorrelation  function  of  a  deterministic  signal,  such  as  a  signal  built  using 
Walsh  functions,  quantifies  the  degree  of  similarity  between  the  signal  and  a  time  shifted 
version  of  itself.  In  contrast  to  a  non-deterministic  signal,  the  autocorrelation  function  of 
a  deterministic  signal  can  be  evaluated  precisely.  In  this  section  we  will  examine  the 
autocorrelation  functions  for  the  Walsh  functions,  and  consider  the  average 
autocorrelation  for  a  set  of  Walsh  functions. 

1.  Autocorrelation  of  Walsh  Functions 

In  order  to  explore  the  autocorrelation  function  of  Walsh  functions,  we  will 
extend  our  definition  and  consider  each  Walsh  function  to  be  periodic  with  a  period  of  T, 
which  is  consistent  with  their  use  in  practice.  On  the  forward  channel,  each  mobile  user 
is  assigned  a  specific  Walsh  function.  The  data  destined  for  that  mobile  user  is  encoded 
by  successively  applying  the  user’s  entire  Walsh  function  to  each  data  bit.  Accordingly, 
the  Walsh  function  is  periodic,  with  period  T  =  bit  duration  =  l/bit  rate ,  and  we  will  see 
that  the  corresponding  autocorrelation  functions  are  also  periodic  with  period  T. 

We  define  the  normalized  autocorrelation  function,  ai  (w),  for  any  periodic 


Walsh  function  w;  (f)  e  WN  by 


Although  (2.13)  does  not  shed  much  light  on  the  Walsh  autocorrelation  functions, 
we  can  generate  them  numerically  based  on  this  definition.  Let  us  look  at  the  set  of 
Walsh  functions  Ws,  depicted  in  Figure  2.6,  and  the  corresponding  autocorrelation 

functions,  at  (w),  shown  in  Figure  2.7. 

As  shown  in  Figure  2.7,  there  are  only  N/ 2  +  1  =  5  unique  autocorrelation 
functions,  since  (^(u)  =  a2(u) ,  a3(u)  =  a4(u) ,  a5(u)  =  a6(u) .  In  spread  spectrum 
systems  we  define  chip  duration,  Tc ,  to  be  the  bit  duration  divided  by  the  number  of 
chips  per  bit,  Tc  =  t/n.  If  we  consider  the  Walsh  functions  of  order  eight,  to  be  periodic 
with  period  T,  and  chip  duration,  Tc  =T/8 ,  we  see  in  Figure  2.7  that, 
w2 (?)  =  w,  (t -  2TC)  and  w4 (?)  =  W3  (? -  Tc)  and  w5  (?)  =  w6 (?  -  27;) .  Since  a  time 

shift  in  a  deterministic  signal  has  no  effect  on  its  autocorrelation  function,  the  Walsh 
functions  that  are  time-shifted  versions  of  each  other  have  identical  autocorrelation 
functions.  We  found  that,  for  this  reason,  in  any  set  WN  of  Walsh  functions  there  are 

exactly  N/2  +  1  unique  autocorrelation  functions.  They  are 

«o («).  aN-M)  and  cc2i{u)  =  a2i_x{u),  for  i  =  l,2,...,N/2 - 1. 

Furthermore,  the  autocorrelation  functions  of  Walsh  functions  are  piecewise  continuous 
and  linear  within  the  intervals  ?  6  [nTc,[n  +  1)7;)  for  n  =  0,±1,±2,... 

While  these  observations  regarding  the  autocorrelation  functions  for  Walsh 
functions  are  interesting,  they  are  not  helpful  in  generating  a  representative 
autocorrelation  function  for  a  set  of  Walsh  functions.  Accordingly,  we  will  develop  an 
average  autocorrelation  function  for  a  set  of  Walsh  functions,  WN . 
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2. 


Average  Autocorrelation  of  a  Set  of  Walsh  Functions  WN 


We  define  the  average  autocorrelation  function,  ( « ) .  for  a  set  WN  of  Walsh 

functions  by 


™  i=0 


where  at(u)  is  defined  by  (2.13).  If  we  generate  r/g  (u)  from  the  Walsh  functions  in  the 
set  Wz,  we  get  the  average  autocorrelation  function  depicted  in  Figure  2.8. 


Figure  2.8  Average  Autocorrelation  Function,  (u)  for  the  Set  Ws . 

We  find  that  the  form  of  ^  (u)  is  similar  to  the  form  of  the  autocorrelation  function  for 

a  random  binary  signal.  A  random  binary  signal  is  generated  from  an  infinitely  long 
binary  random  sequence  in  which  the  bits  are  independent  and  identically  distributed 
random  variables.  The  resulting  normalized  autocorrelation  function  for  the  random 
binary  signal  is  defined  after  [15]  by 
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(2.14) 


P(u)  = 


M 


l-H  M<rc 

71  11 


0, 


otherwise 


This  autocorrelation  function  is  an  aperiodic  function,  since  it  is  generated  from  an 
aperiodic  signal.  On  the  other  hand,  we  know  that  y?z(u)  is  periodic  with  period  T. 


Consequently,  (w)  can  be  described  analytically  over  the  period  (~% > %)  by 


0, 


\u\<T 


rc<|w|<-  =  -^  =  AT 

C  II  2  2 


We  found  experimentally  that  this  result  can  be  generalized  to  a  set  of  Walsh  functions 
of  any  order  N.  That  is,  the  average  normalized  autocorrelation  function,  y?N(u),  can  be 
written  analytically  over  one  period  as  follows: 

^N(u)  =  < 

where  (u)  is  periodic  with  period  T  =  NTC ,  as  depicted  in  Figure  2.9. 

Accordingly,  we  now  have  an  analytical  expression  representing  the  average 
autocorrelation  function,  y?N  (n) ,  of  a  set  WN  of  Walsh  functions.  We  will  also  consider 
another  average  autocorrelation  function  which  excludes  the  identity  Walsh  function, 
w0(r),  from  the  average.  In  a  DS-CDMA  system,  the  identity  Walsh  function  is 
typically  reserved  for  the  pilot  tone,  and  as  such  it  cannot  be  assigned  to  a  mobile  user. 


1-t1,  \u\<Tc 
Tc  '  1 

o,  rc<|«|<-  =  ^ 
11  2  2 


(2.15) 
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Figure  2.9  Average  Autocorrelation  Function,  r4N  (w)  for  the  Set  WN . 

Therefore,  to  determine  the  average  autocorrelation  function  of  Walsh  functions  that  may 
be  assigned  to  mobile  users,  we  should  exclude  w0(f)  and  consequently  a0(u)  from  the 
set,  and  develop  a  different  average  autocorrelation  function.  Let  us  define  the  revised 
function,  y?  n  (u)  ,  which  represents  the  average  autocorrelation  of  the  non-identity 
Walsh  functions,  by 

j  AM 

^ n  {u)  =  ^—^aXu). 

We  return  to  W&  for  an  example,  and  find  the  revised  average  autocorrelation 
function,  r/s  (w) ,  to  be  as  shown  in  Figure  2.10.  The  form  of  y? 8  (n)  is  similar  to  the 

autocorrelation  function  for  pseudorandom  or  pseudo-noise  (PN)  signal.  PN  signals  are 

generated  by  PN  sequences,  which  are  maximal  length  deterministic  sequences.  The  PN 

. — .  ^ — . 

signals  are  periodic  with  period  NTC ,  where  N  is  the  length  of  the  PN  sequence  that 
generated  the  PN  signal.  When  PN  signals  are  used  in  conjunction  with  Walsh  functions 
for  DS-CDMA,  N  »  N .  We  will  discuss  PN  sequences  in  more  detail  in  Chapter  3. 


49 


Figure  2.10  Revised  Average  Autocorrelation  Function,  s  (u)  for  the  Set  Ws . 


The  autocorrelation  function  for  a  PN  signal,  SPN  (u) ,  as  adapted  from  [15],  is 
defined  by 


^ PN  ( u )  ~ 


1-— 

Tc 

-1 

N 


i  +  ij,  0  <  |u|  <  Tc 

T<\u\<^ 


(2.16) 


which  is  periodic  with  period  NTC  as  shown  in  Figure  2.11. 

^ — . 

When  the  length  of  the  generating  PN  sequence  gets  large,  N  — >  °° ,  the 
autocorrelation  function  of  the  PN  signal  approaches  the  autocorrelation  function  of  a 
binary  random  signal,  #PN  ( u )  — >  f3(u) .  Accordingly,  for  simplicity,  the  autocorrelation 

function  of  the  PN  signal,  BPN  (m)  ,  is  often  modeled  as  the  autocorrelation  function  of 
the  random  binary  signal,  /?  (u) ,  as  in  [18]. 
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SPN  (w) 


-t,  o  r 


T  =  NT 


Figure  2.11  Autocorrelation  Function  of  a  PN  Signal. 


We  found  that  for  our  revised  average  autocorrelation  function,  ^n(u)  ,  the  form 


of  (u)  held  true  for  the  general  case  as  shown  in  Figure  2.12. 


-Tc  0  Tc 


T  =  NT 


Figure  2.12  Revised  Average  Autocorrelation  Function,  Ws  («)  for  a  Set  W 


Accordingly,  our  revised  average  autocorrelation  function  for  the  set  of  non¬ 


identity  Walsh  functions  can  be  described  analytically  by 


W N  (w)  = 


u 

1  1 

1 

1 

1  + 

,  0  <  u 

Tc 

N-l) 

-1 

T.<  \u\ 

.N-l 

NT 


(2.17) 


We  now  have  two  possible  analytical  autocorrelation  functions, 
WN  («)  and  Wn  (m)  ,  which  can  represent  the  average  autocorrelation  functions  of  a  set 
of  Walsh  functions,  WN . 


D.  SUMMARY 

In  Chapter  2,  we  investigated  the  definition  and  properties  of  Walsh  functions  as 
they  apply  to  DS-CDMA  systems.  We  explored  the  algebraic  group  properties  of  a  set 
WN  of  Walsh  functions  over  multiplication,  and  we  proved  the  closure  property  using 

Rademacher  functions  as  building  blocks  for  the  Walsh  functions.  We  introduced  the 
notion  of  extending  the  orthogonality  property  for  Walsh  functions  to  ensure  that  the 
integral  of  any  three  non-identity  Walsh  functions  over  the  bit  interval  is  identically  zero. 

We  developed  a  method  to  generate  a  reduced  set  Wn  of  Walsh  functions,  which  extends 
orthogonality  while  maximizing  the  number  of  Walsh  functions  in  the  reduced  set.  We 
explored  the  autocorrelation  functions  for  a  set  of  Walsh  functions,  and  developed  two 

average  autocorrelation  functions,  WN  (u)  and  Wn  (u)  ,  for  the  set  of  Walsh  functions 
and  the  set  of  non-identity  Walsh  functions,  respectively. 
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III.  FORWARD  CHANNEL  MODEL 


The  forward  channel  in  a  DS-CDMA  cellular  system  is  the  channel  that  carries 
traffic  from  the  base  station  to  mobile  users,  as  depicted  in  Figure  3.1.  As  we  transition 
from  cellular  systems  that  carry  principally  voice  traffic  to  systems  that  must  support 
data-type  traffic  at  high  data  rates,  the  forward  channel  becomes  more  critical.  For 
example,  when  mobile  users  attempt  to  download  large  amounts  of  data,  it  will  be  the 
forward  channel  that  carries  this  data.  In  Chapter  HI,  we  will  set  up  a  channel  model  for 
a  typical  DS-CDMA  cellular  system,  complete  with  channel  fading  and  shadowing 
effects  built  into  the  model.  We  will  ultimately  characterize  the  input  to  a  mobile  user’s 
receiver,  which  includes  the  traffic  intended  for  the  user,  together  with  the  noise  and 
interference  terms. 


Forward  Channel 


Figure  3.1  CDMA  Forward  and  Reverse  Channels. 
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Analysis  of  the  forward  channel  is  complicated  by  two  facts;  the  channel  is 
wireless  and  our  users  are  mobile.  These  conditions  give  rise  to  reflection,  diffraction, 
and  scattering  of  the  forward  signal.  The  average  signal  power  received  by  the  mobile 
user  on  the  forward  channel  varies  as  a  function  of  distance  from  the  base  station  and  the 
terrain  between  them.  Additionally,  as  the  forward  signal  emanates  (in  many  directions) 
from  the  base  station,  it  reflects  off  buildings,  trees,  mountains,  etc...  It  arrives  at  the 
receiver  of  the  mobile  user  via  many  paths  of  different  lengths,  and  at  various  signal 
strengths.  The  sum  of  the  multipath  signals  changes  as  the  user  moves.  This  results  in 
rapid  fluctuations  in  the  strength  of  the  received  signal. 

Large-scale  propagation  models,  which  are  used  to  determine  the  average  signal 
power  received,  are  discussed  in  Section  III.B.  Small-scale  fading  models,  which 
characterize  the  rapid  fluctuations  of  the  received  signal  due  to  multipath  effects,  are  also 
discussed  in  Section  IH.B.  We  combine  these  effects  into  our  model  in  Section  m.C. 

A.  BUILDING  THE  DS-CDMA  FORWARD  SIGNAL 

We  model  the  cells  in  our  CDMA  system  as  hexagons,  using  a  basic  seven-cell 
cluster  as  depicted  in  Figure  3.2.  The  hexagon  in  the  center  represents  the  cell  of 
primary  interest  in  our  analysis.  We  assume  that  a  base  station  is  located  in  the  center  of 
each  cell.  We  let  K  represent  the  total  number  of  mobile  users  (or  active  channels)  in  the 
center  cell,  while  Ks  is  the  total  number  of  mobile  users  (or  active  channels)  in  the 
adjacent  cells  i  =  1,2,..., 6 . 
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Figure  3.2  Seven-Cell  Cluster. 

In  this  section,  we  will  describe  the  forward  signal  S  (?)  that  is  transmitted  by  the 

base  station  in  each  cell.  It  includes  the  traffic  for  all  active  channels  in  the  cell.  The 
power  spectrum  of  this  traffic  has  been  spread  by  a  factor  of  N  and  then  BPSK 
modulated.  From  the  perspective  of  an  individual  mobile  user  in  the  center  cell,  the 
transmitted  signal,  50  (?) ,  contains  traffic  intended  for  him  plus  intracell  interference. 
Traffic  that  is  received  by  a  mobile  user  and  that  is  intended  for  other  users  in  his  cell 
constitutes  intracell  interference.  In  addition,  the  signals  Si  (?)  (l  <i<  6)  transmitted 

by  the  base  stations  in  the  six  adjacent  cells  constitute  intercell  interference  to  our  mobile 
user.  Intercell  interference,  or  co-channel  interference,  is  the  traffic  intended  for  mobile 
users  in  the  adjacent  cells. 


55 


1. 


Direct  Sequence  Spread  Spectrum  CDMA 


Spectrum  Spreading  is  a  process  in  which  the  power  spectrum  of  a  base  band 
signal  is  spread  over  a  specified  transmission  bandwidth.  Traditionally,  spectrum 
spreading  has  been  done  for  a  variety  of  reasons,  including  anti-jam  and  low  probability 
of  intercept.  In  DS-CDMA  cellular  systems,  we  use  spread  spectrum  as  a  multiple  access 
technique.  That  is,  given  a  finite  bandwidth  for  many  users  in  the  cellular  system,  we 
optimize  the  use  of  the  bandwidth  by  spreading  all  the  user  traffic  across  the  band  using 
PN  spreading  signals  and  Walsh  functions  for  orthogonal  cover.  Each  user  despreads 
only  that  traffic  which  is  destined  for  him,  using  the  cell’s  PN  signal  and  his  unique 
Walsh  function. 

In  Direct  Sequence  Spread  Spectrum  (DSSS)  systems  the  spreading  of  the 
information  signal  is  accomplished  by  multiplying  the  information  signal  by  a  spreading 
signal.  We  will  represent  the  information  signal  intended  for  mobile  user  k  as  bk  (r) , 

which  is  a  stream  of  binary  data  or  bits,  where  bk  (t)e  |±l]  within  each  bit  interval  of 

duration  T.  We  spread  each  data  bit  by  a  factor  of  N  using  a  spreading  signal  whose  chip 
duration  is  TC=T/N . 

To  illustrate  spectrum  spreading,  we  can  spread  a  random  binary  signal, 
x(r)  — >  xs  (r) ,  by  a  factor  of  N.  As  shown  in  Figure  3.3,  the  main  lobe  of  the  power 

spectral  density  is  spread  by  a  factor  of  N,  while  the  magnitude  at  &  (0)  is  reduced  by  a 

factor  of  l/ N2 .  We  used  a  spreading  factor  of  N=  2  to  develop  Figure  3.3,  while  our  DS- 
CDMA  system  will  use  a  spreading  factor  of  128. 
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Figure  33  Spectrum  Spreading  by  a  factor  of  N. 

As  an  example  of  DSSS,  let  us  look  at  the  first  two  bits  of  data  intended  for 
mobile  user  3.  We  will  use  the  two-bit  binary  sequence  (01)  whose  waveform  as  shown 
in  Figure  3.4  is  defined  by 

f  J.  te(0,T) 

K  J  [-1,  te(T,2T) 

To  demonstrate  direct  sequence  spectrum  spreading  we  can  use  a  Walsh  function 
w,  (f)e  to  spread  b~(t)  by  a  factor  of  N= 4.  By  multiplying  each  bit  in  b3(r)  by 

w3  (?),  we  get  a  spread  signal  that  represents  the  information  in  b3(t)  as  shown  in 
Figure  3.4.  The  new  signal  has  a  power  spectrum  that  is  spread  by  a  factor  of  four. 
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Figure  3.4  Direct  Sequence  Spreading  of  Two  Data  Bits  Using  w3  (t ) . 

We  can  also  view  the  spreading  process  in  terms  of  sequences,  if  we  expand  our 
information  sequence  so  that  each  information  bit  is  represented  by  N=  4  elements  in  the 
information  sequence.  That  is,  b3  =  (OOOOllll) .  The  spreading  process  of  the  sequence 

b3  is  a  simple  modulo  two  addition  with  the  Walsh  sequence,  w3 ,  as  depicted  in  Figure 
3.5. 


Figure  3.5  Spectrum  Spreading  of  Two  Data  Bits  Using  Sequences. 
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In  DS-CDMA  systems,  PN  signals  are  used  in  conjunction  with  Walsh  functions 
to  spread  the  data  signal.  In  the  same  manner  that  Walsh  functions  combat  the  intracell 
interference,  the  PN  signals,  to  a  lesser  degree,  combat  the  intercell  interference.  That  is, 
each  base  station  is  assigned  a  unique  PN  signal.  These  PN  signals  can  be  made  to  be 
orthogonal  to  one  another,  but  are  not  necessarily  so.  A  truly  random  PN  signal  for  each 
base  station  would  be  ideal  but  impractical,  since  it  must  be  known  to  both  the  base 
station  and  all  the  potential  users.  Accordingly,  we  generate  PN  sequences 
deterministically  from  linear  feedback  shift  registers  (LFSR). 

We  represent  the  PN  signal  in  the  center  cell  as  c(t) ,  and  adjacent  cell  PN  signals 

as  ct  ( t )  for  i=  1,2,... 6.  Since  the  mobile  users  in  the  center  cell  are  not  synchronized 

with  base  stations  in  adjacent  cells,  the  use  of  even  orthogonal  PN  signals  does  not 
completely  eliminate  the  intercell  interference  in  the  same  manner  that  Walsh  functions 
eliminate  intracell  interference.  However,  the  use  of  PN  signals  does  ensure  that  the 
intercell  traffic,  which  results  in  intercell  interference,  remains  spread  over  the  band  and 
minimizes  it's  effect  on  SNR. 

Application  of  the  PN  signal  also  ensures  that  all  of  the  user  channels  are  spread 
the  same  chip  duration,  Tc  =  T/N .  The  application  of  the  Walsh  function  w,(r)e  WN 
alone  does  not  guarantee  a  spreading  factor  of  N.  For  example,  if  we  had  applied 
w,(t)  vice  w3(r)  to  spread  data  for  mobile  user  1,  where  by  (?)  =  £>3  (?) ,  the  spreading 
would  have  only  been  by  a  factor  of  two  rather  than  four,  as  shown  in  Figure  3.6. 
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bx  =  (00001111) 


wx=  (00110011) 


Y 

©“► 


bx  ®wx  =(00111100) 


Figure  3.6  Spectrum  Spreading  of  Two  Data  Bits  Using  \nx  (t ) . 


Accordingly,  the  use  of  Walsh  functions  provides  for  orthogonal  cover  for 
intracell  interference,  while  the  PN  signal  accomplishes  the  true  spreading  and  minimizes 
intercell  interference. 

The  baseband  spreading  and  despreading  process  for  one  data  bit  intended  for 
user  one  in  our  DS-CDMA  systems  is  depicted  in  Figure  3.7.  We  assume 
synchronization  with  the  base  station. 


Spread  Traffic 


Figure  3.7  DS-CDMA  Spreading  and  Despreading. 


We  can  recover  6,  (t)  from  the  spread  traffic  bx  (r)c(r)  because  Wj  (?)  w,  (r)  =  1 
and  c(t)c(t)  =  1. 
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If  we  look  at  the  process  in  terms  of  sequences,  we  can  build  a  table  which  shows 
the  spreading  of  one  bit  of  information  bx  =  (0000) ,  using  both  the  Walsh  function, 

wx  =  (001 1) ,  and  a  PN  sequence  (for  this  particular  bit  interval),  c  =  (101 1) . 


(0000) 

wx 

(0011) 

bj  ©  wx 

(0011) 

c 

(1011) 

bJ®wl®c 

(1000) 

Table  3.1  Spreading  of  One  Information  Bit  using  wl  and  c . 

At  the  receiver,  in  order  to  despread  the  sequence  and  recover  bx ,  we  apply  the  PN  and 
Walsh  sequences  again  as  shown  in  Table  3.2. 


bj  ©  wx  ©  c 

(1000) 

C 

(1011) 

bj  ©  wx®  c  ©  c  =  bj  ©  wx 

(0011) 

wx 

(0011) 

bj  ®w1®wl=  bj 

(0000) 

Table  3.2  Despreading  of  One  Information  Bit  using  wx  and  c . 
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When  another  user  with  assigned  Walsh  function  w2(?),  for  example,  in  the 
center  cell,  applies  the  center  cell’s  PN  signal  and  his  Walsh  function,  c(?)w2(?) ,  to  the 
spread  traffic,  the  output  is  still  spread,  as  shown  in  Table  3.3. 


bj  ®  wx  ©  c 

(1000) 

c 

(1011) 

bJ®w1®c®c  =  bJ®wl 

(0011) 

w2 

(0110) 

b1@wl®  w2 

(0101) 

Table  3.3  Application  of  w2  to  intracell  traffic  covered  with  wx . 


Alternatively,  if  a  user  with  assigned  Walsh  function  Wj  (?)  in  an  adjacent  cell 
i-6,  for  example,  received  the  spread  traffic  and  applied  his  PN  signal  and  Walsh 
function,  c6(t +  r)w1(t +  t)  ,  the  output  would  also  remain  spread, 

bx  (?)  c(?)  Wj  (?)  c6(?  +  t)  wx  (?  +  t)  .  Here,  t  represents  a  timing  discrepancy  due  the  lack 

of  synchronization  between  the  center  cell  and  adjacent  cell  i=6.  In  fact,  this  output 
would  likely  be  spread  even  further  since  the  product  of  nonsynchronous  PN  signals 
and/or  Walsh  functions  would  contain  some  pulses  with  duration  less  that  Tc  as  shown  in 

Figure  3.8.  Therefore  the  spreading  factor  would  be  greater  than  N,  since  the  spread  of 
the  power  spectrum  is  based  upon  the  duration  of  the  smallest  pulse. 
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Figure  3.8  Product  of  Non-Synchronous  Walsh  Functions. 

2.  The  Transmitted  Signal  S  (t ) 


The  base  station  in  each  cell  must  convert  the  information  signal  for  user  k  into  a 
DS-CDMA  signal,  amplify  it  to  the  desired  power  level,  and  BPSK  modulate  it  for 
transmission,  as  shown  for  the  center  cell  in  Figure  3.9. 

W* 

K(t)  - yx> - *x) — ► - KX) - >\U) 

\Nk(t)  c(t)  COS  (2  7lfct) 

Figure  3.9  Spreading  and  Modulating  the  Information  Signal. 

Accordingly,  the  spread  spectrum  orthogonally-covered  BPSK-modulated  signal  ok  (t) 
that  will  be  transmitted  for  user  k  is  defined  by 
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where 


k  =  the  mobile  user  or  channel  k  in  the  center  cell, 

Puk  =  the  average  transmit  power  in  the  kth  channel, 

bk  (?)  =  the  information  signal  consisting  of  binary  data  stream 
for  the  M  user  channel  in  center  cell, 
wt  (f)  =  the  Walsh  function  for  the  &th  user  channel  in  center  cell, 
c(r)  =  PN  spreading  signal  for  the  center  cell,  and 
fc  =  the  carrier  frequency  of  the  signal. 


The  base  station  in  the  center  cell  generates  a  signal,  as  defined  in  (3.1),  for  every 
user  or  active  channel  in  the  cell,  to  form  a  composite  signal  S0  (r)  that  is  transmitted  into 
the  channel  as  follows: 

*0  (o (o = xV2/^*  W' w*  (0c(0c°s(2?r/c0’  (3-2) 

*=0  k= 0 

where  K  is  the  number  of  active  channels  in  the  center  cell. 

Once  the  signal  ^  (t)  enters  the  channel,  it  is  attenuated  due  to  path  loss  and  it  is 

subject  to  fading  and  shadowing.  In  the  next  two  sections  we  will  discuss  large-scale 
path  loss  and  small-scale  fading  of  the  forward  channel  as  they  are  used  in  our  model. 


B.  PROPAGATION  IN  THE  MOBILE  RADIO  CHANNEL 

The  signal  transmitted  from  the  base  station  to  the  mobile  users  loses  a  portion  of 
its  power  along  the  way  as  a  function  of  the  distance  it  must  travel  and  the  terrain  type 
across  which  it  travels.  This  loss  in  signal  power  is  called  path  loss,  which  we  will 
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generally  represent  as  L  or  £  in  dB .  We  use  large-scale  propagation  models  to 
determine  an  average  and/or  median  power  received  based  on  the  separation  distance 
between  the  transmitter  and  receiver  and  a  general  notion  of  the  surrounding  terrain/area, 
as  discussed  in  Sections  m.B.l  and  IH.B.2.  These  models,  however,  do  not  account  for 
that  fact  that  mobile  users  at  a  specific  distance  with  the  same  general  type  of  terrain  may 
have  local  (dynamically  changing)  environments  that  are  vastly  different  from  one 
another.  These  differences  result  in  variability  of  the  path  loss  predicted  by  our  large- 
scale  propagation  models.  We  model  this  variability  as  a  lognormal  random  variable  in 
Section  DI.B.3.  Small-scale  propagation  models  are  used  to  model  the  effects  of 
multipath  propagation  and  are  presented  in  Section  III.B.4. 


1.  Large-Scale  Path  Loss 

Large-scale  propagation  models  are  based  on  the  type  of  environment  through 
which  the  signal  must  travel.  For  example,  in  communications  systems  where  there  is  a 
direct  line-of-sight  between  transmitter  and  receiver  we  use  the  free-space  propagation 
model  as  in  [19].  We  can  predict  the  power  received  using  the  Friis  free-space  equation 
defined  by 


(4x)2d%' 


(3.3) 


where 
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Pb  =  the  free  space  power  received, 

Pt  =  the  transmitted  power, 

Gt  =  the  transmitter  antenna  gain, 

Gr  =  the  receiver  antenna  gain, 

Xc  =  the  wavelength  of  the  carrier,  Xc=  — , 

fe 

d  -  the  separation  distance  between  transmitter  and  receiver, 

Ls  =  the  system  (hardware)  loss  factor,  and 
fe  =  the  carrier  frequency. 

Since  our  aim  is  to  model  the  channel  without  regard  to  antenna  gains  and  hardware 
losses,  we  will  take  Gt  =  Gr  =  Ls  =  1 .  The  free-space  path  loss  <£fr  can  be  modeled  as 

4(d)  =  10log^.  =  -20logf-i-l  (dB).  (3.4) 

Pb  1^4  nd) 

The  free-space  power  received  PfT,  as  described  by  (3.3),  and  consequently  the 
free-space  path  loss  calculated  by  (3.4),  are  limited  to  distances  d  that  are  in  the  far  field 
of  the  transmitting  antenna.  Accordingly,  large-scale  models  use  a  power  reference  point 
calculated  at  a  reference  distance  d0 ,  which  can  then  be  related  to  the  free-space  power 
received  by 

=/»(d0)(AT  a ,<d. 

\d  ) 

The  free-space  reference  power,  Pb  (d0) ,  can  either  be  based  upon  close-in 
measurements  or  calculated  using  (3.3). 

The  mobile  communication  channel  cannot  be  accurately  modeled  using  only  the 
free-space  propagation  model.  Typically  a  line-of-sight  between  the  base  station  and  the 
mobile  user  does  not  exist,  due  to  buildings  and  other  obstacles  that  lie  in  the  path. 
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Propagations  models  for  other  than  free  space  have  been  developed  both  analytically  and 
empirically.  Both  approaches  indicate  that  average  received  power  decreases 
logarithmically  as  a  function  of  distance  [19].  The  average  path  loss  can  be  described  as 
a  function  of  a  path-loss  exponent  n  using 


M4 


fdX 

■dn 


(3.5) 


The  path-loss  exponent,  n,  characterizes  the  rate  of  path  loss  and  is  dependent 


upon  the  propagation  environment  such  as  terrain,  and  obstructions  to  the  path.  Table  3.1 
from  [19]  shows  typical  path-loss  exponents  for  various  conditions. 


Environment 

Path-loss  exponent,  n 

Free  Space 

2 

Urban  area  cellular  radio 

i 

2.7  to  3.5 

Shadowed  urban  cellular  radio 

i 

3  to  5 

In  building  line-of-sight 

1.6  to  1.8 

Obstructed  in  building 

4  to  6 

Obstructed  in  factories 

2  to  3 

Table  3.4  Path-loss  exponents  in  Various  Environments  from  [19]. 


By  applying  the  appropriate  path-loss  exponent  to  our  propagation  model,  we  can 
determine  an  average  path  loss,  L„ ,  which  represents  the  ensemble  average  of  all 
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possible  path  losses  for  the  given  distance  d.  This  average  path  loss,  Ln,  or  £„  in  dB, 
can  be  described  by 


J£„(rf)  =  je(d0)  +  10#*log 


(dB), 


(3.6) 


V  V 


where  <£(d0)  is  average  loss  at  the  power  reference  distance  which  can  be  specifically 
measured  or,  based  on  a  free-space  assumption,  calculated  using  (3.4)  as 
•£(d0)  =  <£fr  (d0) .  The  key  to  using  (3.6)  in  analysis  is  determining  what  path-loss 

exponent  is  appropriate  for  a  given  environment.  While  Table  3.4  is  a  good  start,  we  can 
be  more  specific  about  our  operating  environment  and  consequently  the  loss  associated 
with  it.  Accordingly,  in  the  next  section  we  will  explore  the  Hata  model,  which  allows 
for  more  specificity. 


2.  Hata  Model 

Propagation  models  are  typically  developed  using  both  analytical  and  empirical 
methods.  Empirical  methods  use  measured  data  to  produce  a  graphical  or  analytical 
expression,  which  well  represents  the  actual  data.  As  pointed  out  in  [19],  one  of  the 
simplest  and  most  accurate  empirical  models  for  path  loss  in  cluttered  cellular  systems  is 
Okumura’s  model.  Okumura’s  model  is  based  entirely  on  measured  data,  which  is 
represented  graphically  in  [20]  as  a  set  of  curves  from  which  median  path  loss  can  be 
predicted. 

The  Hata  model  takes  Okumura’s  graphical  path  loss  data  and  represents  it  with  a 
standard  formula  for  an  urban  operating  area  which  can  be  modified  for  other  areas  with 
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correction  equations  as  presented  in  [21].  Hata’s  original  model  closely  matches  the 
predictions  of  the  Okumura  model  for  distances  greater  than  1  km.  In  [22],  the  Hata 
model  was  extended  to  cover  distances  on  the  order  of  1  km  and  higher  frequency  ranges. 
We  use  the  extended  model  since  it  better  represents  the  size  and  frequency  ranges  of 
future  cellular  systems. 

The  extended  Hata  model,  or  COST-Hata  model,  uses  the  carrier  frequency  fc , 


the  heights  of  the  base  station  and  mobile  antennae,  /zbase  and  /imobile ,  the  distance  between 
the  base  station  and  mobile,  d,  and  an  operating  area  correction  factor  CM  to  predict 
median  path  loss  in  dB  as  defined  in  [22]  by 


^  =  46.3  +  33.9!og-^-  -  13.82log^  -a(/imobile ) 
MHz  ~ 


44.9-6. 55log 


K 


m 

— ^iog-r~+cM> 


m 


km 


(3.7) 


where 


^mobile 


m 


1.56log— — 0.8 
Mhz 


(dB), 


and 


[0  dB,  for  medium  sized  city  and  suburban  areas 
M  [3  dB,  for  metropolitan  centers. 

As  indicated  in  (3.7),  frequency  is  measured  in  MHz,  antennae  heights  are 
measured  in  meters  (m),  and  separation  distance  is  measured  in  kilometers  (km).  The 
extended  Hata  model  is  restricted  to  the  following  range  of  parameters: 
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fc  :  1500  Mhz  to  2000  Mhz, 
h base  :  30  m  to  200  m, 

^1=  :  1  m  to  10  m, 
d  :  1  km  to  20  km. 


The  Hata  median  loss,  ^ ,  predicted  by  (3.7)  is  comparable  to  the  average  loss, 
£n ,  predicted  by  (3.6)  as  shown  in  Figure  3.10. 


Figure  3.10  Comparison  of  Hata  Loss  and  Average  Loss. 

We  computed  the  losses  depicted  in  Figure  3.10  for  a  range  of  distances, 
1  <  d  <  5  at  a  frequency,  fc  =  2000  MHz  .  We  computed  the  Hata  loss,  ,  for  a  large 
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city  with  the  base  station  antenna  height  at  30m  and  the  mobile  antenna  at  lm.  For  the 
average  path  loss  L„ ,  we  used  a  range  of  path-loss  exponents,  which  correspond  to  a 
shadowed  large  city  environment  from  Table  3.4.  We  used  a  power  reference  distance, 
d0 ,  of  10m. 

It  is  clear  that  at  any  specific  distance,  there  exists  a  reasonable  path-loss 
exponent,  n,  such  that  =  £„  .  That  is,  from  Figure  3.10  given  any  specific  distance  d, 

we  could  adjust  n  so  that  the  curves  representing  ,£H  and  ln  cross  at  d.  Accordingly, 
we  use  the  median  path  loss  predicted  by  the  Hata  model  in  a  similar  way  that  the 
average  path  loss  calculated  using  a  path-loss  exponent  is  used. 

3.  Lognormal  Shadowing 

The  path  loss  predicted  by  the  above  methods  represents  an  average  or  median 
value  for  a  given  distance.  However,  the  actual  path  loss  for  mobile  users  in  different 
locations  but  at  the  same  distance  from  the  base  station  can  vary  due  to  terrain 
differences,  clutter,  or  foliage  along  their  specific  path.  We  call  these  variations  in  large- 
scale  propagation  shadow  fading,  or  simply  shadowing.  Shadowing  causes  the  actual 
path  loss  at  a  given  distance  to  differ  significantly  from  the  predicted  average  or  median 
path  loss.  These  differences,  measured  in  decibels,  can  be  represented  by  a  Gaussian 
random  variable,  since  the  measured  path  loss  varies  around  a  distance  dependent 
average  or  median  value.  Our  path  loss  with  shadowing  £%  in  terms  of  decibels  is 

modeled  by 

£,(<*)  =  £(<*)  +  £  (dB),  (3.8) 
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where  9C  is  a  zero-mean  Gaussian  random  variable  %  ~  Af  (0,crdB)  with  standard 
deviation  adB  and  where  £(d)  is  a  distance-dependent  path  loss.  Since  JS(cf)  is  a 
deterministic  value  and  9C  is  a  zero-mean  Gaussian  random  variable,  our  shadowed  path 
loss  is  a  Gaussian  random  variable  £%  ~  JV (£,crdB)  where  £  is  the  mean  and  the  median 

value  of  the  random  variable.  The  median  value,  or  50th  percentile,  of  a  random  variable 
is  that  smallest  value  of  the  random  variable  where  its  cumulative  distribution  function 
(cdf)  is  1/2.  In  terms  of  the  probability  distribution  function  (pdf),  the  median  is  the  point 
at  which  the  area  under  the  pdf  to  the  right  of  the  point  is  equal  to  the  area  under  the  pdf 
to  the  left  of  the  point.  For  the  Gaussian  random  variable,  the  mean  and  median  values 
are  coincident.  Accordingly,  as  model  for  £(d)  in  (3.8),  we  can  use  the  predicted 

average  path  loss  £„  defined  by  (3.6)  as  in  [19],  or  a  predicted  median  path  loss  as  in 
[15].  In  our  model,  we  will  use  the  median  path  loss  .£H  as  predicted  by  the  extended 
Hata  model,  which  is  defined  by  (3.7).  Accordingly,  we  revise  (3.8)  to  incorporate  the 
Hata  model  as  follows: 

^=J£h+9C  (dB).  (3.9) 

Our  shadowed  path  loss  £%  is  a  Gaussian  distributed  random  variable 
(^H  ’  ^dB  )  with  mean  fidB  =  ,£H  and  standard  deviation  oiB ,  since  it  is  the  sum  of  a 

Gaussian  random  variable  9C  and  a  deterministic  value  .  We  can  convert  (3.9)  from 
decibel  to  ratio  format  using 

=10logLx, 

«£h  =10logLff,  and 
9C  =  10logX. 
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The  Gaussian  random  variable  9C  ~  JV(0,adB)  (in  decibels)  is  converted  to  a  lognormal 
random  variable  X  ~  A(0,AcrdB),  with  parameters  fix  =  XndB  =  0  and  ox  =  XadB , 


where  X  =  In  10/10  as  derived  and  defined  in  Appendix  DI-A.l.  Furthermore,  our 
shadowed  path  loss  in  ratio  form,  Lx  ,  has  been  converted  to  a  lognormal  random 
variable  Lx  ~  A(2£H,/l<Tdfi)  as  well,  which  yields  the  lognormal  random  variable  Lx 
described  by 

Lx=LhX.  (3.10) 


The  parameters  nx  and  ox  of  the  lognormal  random  variable  X  ~A (jux,crx) 

are  not  the  distribution’s  mean  and  standard  deviation.  We  define  the  pdf  and  the 
moments  for  the  lognormal  random  variable  in  Appendix  HI-A.2. 

We  have  created  a  model  for  path  loss  that  incorporates  the  Hata  model  and  the 
effects  of  lognormal  shadowing.  In  order  to  fold  this  path-loss  model  into  our  forward 
signal,  we  must  look  at  the  link  analysis.  In  general,  we  know  that  the  power  received, 
PT ,  will  be  related  to  transmitted  signal  power,  gains  of  the  antennae,  and  the  system  and 
path  losses  as  follows: 


R  = 


mg, 

LLS 


(3.11) 


where  we  have  normalized  the  antennae  gains  and  system  loss  to  one.  However,  in  our 
forward  signal  model,  we  need  the  flexibility  to  adjust  the  signal  power  allocated  to 
individual  channels.  For  example,  we  may  want  to  increase  the  power  in  the  pilot 
channel  to  enhance  synchronization  between  the  base  station  and  the  mobile  users,  and 
consequently  aid  the  demodulation  of  the  data  signals.  If  we  use  Pt  as  our  baseline 
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signal  power,  meaning  most  channels  will  be  transmitted  with  this  signal  power,  then  we 
can  relate  the  signal  power  JPtk  in  each  channel  k  to  this  baseline  using  a  power  factor 

as  follows: 

(3-12) 

In  general,  our  power  factor  will  be  =  1  for  all  channels,  unless  specifically  stated 
otherwise.  We  can  now  modify  (3.11)  to  be  channel-specific  using  (3.12)  as  follows: 


L  L 


(3.13) 


Finally,  if  we  apply  our  Hata-lognormal  model  L  =  LHX  for  path  loss,  we 
transform  the  power  received  in  the  £-th  channel  into  a  random  variable  Pk  as  follows: 


P  _ 

k~LH{d)X ’ 


(3.14) 


where 


=  the  power  factor  used  to  adjust 
the  power  in  the  £-th  channel, 

JPt  =  the  baseline  signal  power, 

Lh  =  the  median  path  loss  using  the  Hata  model,  and 
X  =  the  Lognormal  random  variable  A  ( 0,AcrdB ) . 

As  shown  in  Appendix  HI-C,  the  inverse  of  a  lognormal  random  variable  is 
another  lognormal  random  variable.  Accordingly,  the  power  received  in  the  k th  channel 

is  a  lognormal  random  variable Pk  ~  A^juPk ,AadB ) ,  where  juPt  =ln [<^FkJPjLH). 
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In  this  section,  we  have  developed  a  model  for  the  power  received  by  a  mobile 
user  that  incorporates  large-scale  propagation  losses  and  lognormal  shadowing.  In  the 
next  section  we  will  model  the  effects  of  small-scale  fading. 


4.  Small-Scale  Fading  Due  to  Multipath 

Small-scale  fading  of  a  received  signal  is  characterized  by  rapid  fluctuations  in 
the  strength  (or  amplitude)  of  the  signal.  The  fluctuations  are  a  result  of  multiple  copies 
of  the  transmitted  signal  bouncing  off  various  obstacles  and  arriving  at  the  receiver  via 
different  paths  as  depicted  in  Figure  3.11.  The  copies  of  the  signal  travel  along  paths  of 


Figure  3.11  Multipath  Propagation. 
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different  lengths  and  therefore  sustain  losses  of  different  amounts.  When  they  arrive  at 
the  receiver,  some  copies  are  time-delayed  as  a  result  of  the  longer  paths  that  they  have 
traveled,  and  some  of  the  copies  have  been  reflected  off  moving  objects,  which  induce  a 
doppler-shift  in  the  frequency  of  the  signal.  The  receiver  must  process  the  sum  of  the 
multipath  components  in  order  to  extract  the  information  bit  from  the  signal. 

The  success  or  failure  of  the  receiver  to  correctly  process  each  information  bit  (or 
symbol)  depends  upon  the  small-scale  fading  characteristics  of  the  channel  as  compared 
with  the  signal  characteristics  such  as  bandwidth  and  bit  (or  symbol)  period.  Channel 
fading  characteristics  are  classified  according  to  two  independent  fading  phenomena  as 
described  in  [19]: 

❖  Time  Dispersion  -  fading  due  to  the  collective  time  delay  of  the  multipath 
components  of  the  received  signal.  (Flat/Frequency  Selective) 

❖  Frequency  Dispersion  -  fading  due  to  the  collective  doppler  spread  in  the 
multipath  components  of  the  signal.  (Fast/Slow) 

The  changes  in  signal  amplitude  due  to  doppler  spread  over  one  bit  period  can  be 
characterized  as  fast  or  slow  fading.  Fast  fading  is  caused  by  frequency  dispersion  in  the 
signal,  which  results  in  distortion  of  the  signal.  If  the  frequency  dispersion  is  great 
enough  to  dramatically  affect  the  amplitude  of  the  received  multipath  signal  over  the  bit 
period,  then  it  is  called  fast  fading.  In  reality,  fast  fading  environments  have  only  been 
observed  [19]  when  dealing  with  communication  systems  with  very  low  data  rates, 
(which  result  in  very  long  bit  periods  as  compared  to  the  rate  of  change  of  the  channel.) 
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In  a  slow-fading  channel,  the  doppler  spread  of  the  multipath  components  of  the 
signal  does  not  significantly  affect  the  amplitude  of  the  received  signal  over  the  bit 
period.  Accordingly,  we  consider  the  channel  to  be  static  (for  the  bit  duration)  with 
regard  to  frequency  dispersion/distortion  effects.  Since  our  model  is  concerned  with  high 
data  rates,  our  channel  is  modeled  as  slow-fading. 

The  effects  on  the  received  signal  amplitude  due  to  time  delay  or  time  dispersion 
of  the  multipath  components  of  the  signal  result  in  either  flat  or  frequency-selective 
fading.  Frequency-selective  fading  occurs  when  the  collective  time  delay  of  the 
multipath  components  exceeds  the  bit  duration  and  consequently  (negatively)  affects  the 
processing  of  the  next  information  bit;  this  is  called  Inter-Symbol  Interference  (ISI).  If 
we  think  of  converting  the  multipath  delay  time  of  the  channel  into  the  frequency 
domain,  we  can  describe  the  channel  by  a  corresponding  characteristic  bandwidth.  The 
channel  bandwidths  are  inversely  proportional  to  the  time  delay  spread,  which  means 
long  time  delays  correspond  with  narrow  channel  bandwidths.  When  the  channel 
bandwidth  is  less  than  the  bandwidth  of  the  transmitted  signal,  we  have  frequency- 
selective  fading  and  ISI  in  the  channel. 

On  the  other  hand,  in  a  flat-fading  environment,  the  time  dispersion  of  the 
multipath  signal  is  less  than  the  bit  period,  and  the  receiver  is  able  to  resolve  the 
multipath  signal.  Rat-fading  channels  cause  the  amplitude  of  the  received  signal  to  vary 
with  time  over  the  bit  due  to  the  effects  of  the  multipath  propagation,  but  don’t  interfere 
with  the  adjacent  information  bits.  The  channel  is  time-varying  because  it  changes  as  the 
receiver  moves,  and  as  reflectors  or  scatterers  (of  the  multipath  components)  along  the 
paths  move  or  change.  This  characterizes  multipath  propagation.  The  most  common 
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method  for  modeling  the  amplitude  variations  in  a  flat  fading  channel  is  by  assuming  that 
the  amplitudes  are  distributed  as  a  Rayleigh  random  variable  R  as  defined  in  Appendix 

B. 

Accordingly,  we  will  use  the  Rayleigh  flat-fading  channel  model  to  represent  the 
small-scale  fading  effects  due  to  multipath  propagation  in  our  channel  model.  In  the  next 
section  we  will  combine  the  effects  of  large-scale  path  loss  and  small  scale  fading  into  a 
Rayleigh-lognormal  channel  model. 

C.  RAYLEIGH-LOGNORMAL  CHANNEL  MODEL 

In  Section  IH.B,  we  described  the  large-scale  and  small-scale  propagation  effects 
that  are  possible  when  a  signal  passes  though  a  mobile  channel.  In  Section  m.C  we  will 
combine  these  effects  to  develop  the  Rayleigh-lognormal  channel  model.  The  Rayleigh- 
lognormal  channel  model  can  be  described  as  a  slow-flat-Rayleigh  fading  channel  with 
lognormal  shadowing  and  with  path  loss  as  predicted  by  the  extended  Hata  model.  We 
apply  the  Rayleigh-lognormal  channel  model  to  the  signals  in  the  DS-CDMA  cellular 
system,  and  develop  a  received  signal  by  a  mobile  user  in  center  cell. 


1.  The  Forward  Signal  s0  (t) 

In  order  to  develop  the  effects  of  our  channel  on  the  transmitted  signal  S0  (r) ,  as 
defined  by  (3.2),  we  will  consider  the  signal  in  its  complex  envelope  form  as  follows: 

^(r)  =  Re{^(t)^'2^},  (3.1-5) 
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where  ^  (?)  is  the  complex  envelope  of  the  transmitted  signal  defined  by 

<3'16> 

k=  0 

As  described  in  Section  RB,  we  know  that  the  complex  envelope  S0(t)  is 
affected  by  small-scale  fading  as  modeled  by  the  Rayleigh  random  variable  as  R  -50  (?) . 
The  large-scale  path  loss  endured  by  the  signal  together  with  the  lognormal  shadowing  is 
modeled  by  (t)/JZ^X  ,  using  the  path  loss  as  predicted  by  the  extended  Hata  model. 
Furthermore,  there  is  a  relative  time  delay  rd  introduced  by  the  channel  and  a  phase 
discrepancy  6d  between  transmitting  stations  and  the  mobile  user.  Finally,  the  channel 
introduces  additive  noise  n(t )  and  interference  g(t)  to  the  total  received  signal  r(r), 

which  we  investigate  in  later  sections.  The  complete  baseband  channel  model  is 
graphically  depicted  in  Figure  3. 12. 


Figure  3.12  Baseband  Rayleigh-Lognormal  Channel  Model. 
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We  have  used  the  symbols  s0  (t},  g(t),  n(?),  and  r(?)  to  represent  the  lowpass 
equivalent  signals  for  s0 (?),  £■(?),  n(?),  and  r(?) . 

The  forward  signal  s0  (?)  that  will  be  received  by  a  mobile  user  in  the  center  cell 
as  a  result  of  the  transmitted  signal  S0(?)  passing  through  the  Rayleigh-lognormal 
channel  model  has  a  complex  envelope  described  by  the  following: 

K- 1  p 

k= o  v  (3.17) 

k~0 

k= 0 

where  Pk  is  defined  by  (3.14).  Further,  since  the  time  and  phase  delays  are  relative 

amongst  the  base  stations,  we  take  rd  =  6d  =  0  for  the  forward  signal  from  the  center 

cell’s  base  station.  We  can  describe  the  forward  (bandpass)  signal  using  its  complex 
envelope  as 

,0(?)  =  Re{5~(?)^} 

4-1  , _  (3.18) 

=  ZRpPkbk  (?)  w*  (?)c(?)cos(2^-/c?). 

0 

A  mobile  user  extracts  his  traffic  from  the  composite  signal  s0  (?)  as  described  in 
Section  HI. A.  However,  as  shown  in  Figure  3.12,  the  signal  received  by  a  mobile  user 
includes  Additive  White  Gaussian  Noise  (AWGN)  n(?)  and  co-channel  (or  intercell) 
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interference  g(t)  from  the  adjacent  cells  in  addition  to  the  forward  signal  s0(t) .  We 
develop  the  co-channel  interference  term  in  the  next  section. 

2.  Co-Channel  Interference  ?(?) 

Just  as  the  mobile  user  receives  the  composite  signal  s0  (?)  from  the  base  station 
in  his/her  (the  center)  cell,  he/she  also  receives  a  similar  composite  signal  si  (?)  from  the 

base  station  in  each  adjacent  cell.  The  sum  of  these  constitutes  co-channel  interference. 
The  base  station  in  each  adjacent  cell  generates  a  forward  signal  for  users  in  their  cell  that 
passes  through  the  channel  and  is  received  by  a  user  in  the  center  cell  as  defined  by 

S,(t)  =  ^Riyl2^bij(t+Ti)\Nv(t  +  ti)ci{t  +  Ti)COs(2^fct  +  (pi),  (3.19) 

> 0 

with 

i  =  the  adj  acent  cells  i  =  1,2,.. .,6, 

ij  =  mobile  user  or  channel,/  in  adjacent  cell  i, 

Kt  =  the  number  of  active  channels  in  adjacent  cell  i, 

R.  =  Rayleigh  fading  random  variable  for  signals  from  adjacent  cell  i, 

Py  =  Lognormal  Random  Variable  representing  the  average  power 

received  from  the yth  channel  in  adjacent  cell  i  as  defined  below, 
btj  (?)  =  the  information  signal  for  the  yth  user  channel  in  adjacent  cell  i, 

VJy  (?)  =  Walsh  function  for  the  yth  user  channel  in  adjacent  cell  i, 
c.  (?)  =  PN  spreading  signal  for  the  adjacent  cell  i,  and 

fc  =  the  carrier  frequency  of  the  signal. 

z.  =  the  time  delay  from  adjacent  cell  i,  relative  to  the  time  delay  from 
the  center  cell  base  station. 

(p.  -  the  phase  delay  from  adjacent  cell  i,  relative  to  the  phase  delay 

from  the  center  cell  base  station. 
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We  adapted  Pi;  from  (3.14)  to  represent  the  average  power  received  by  a  mobile  user  in 
the  center  cell  from  forward  signal  st  (?)  in  adjacent  cell  i  channel  j  as  follows: 


MAR’ 


(3.20) 


where 

=  the  factor  used  to  adjust  the  power  in  the y-th 
channel  in  adjacent  cell  i, 

Lh(D  )=  the  median  path  loss  using  the  Hata  model  at  a 
separation  distance  Dt, 

Di  =  the  distance  separating  the  receiver  from  the 
base  station  in  adjacent  cell  i,  and 
X i  =  Lognormal  random  variable  A  (0 ,erdB ) . 

The  co-channel  interference  £  (?)  is  the  sum  of  the  composite  signals  from  each  of 


the  adjacent  cells  as  follows: 

?(0  =  E  fJRi^bij(t  +  Ti)yNij(t  +  ri)ci(t  +  Ti)cos(2^fct  +  (pi).  (3.21) 

i=l  j= 0 

With  the  co-channel  interference  term  defined,  we  will  now  proceed  to  define  the 
complete  signal  received  by  a  mobile  user,  r  (?) . 


3.  The  Received  Signal  r(?) 

As  described  in  Section  m.C.1,  the  signal  received  by  a  mobile  user  in  the  center 
cell  is  comprised  of  the  forward  signal  s0  (?)  sent  from  the  user’s  base  station. 
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interference  g  (t)  from  the  six  adjacent  base  stations,  and  Additive  White  Gaussian  Noise 


(AWGN)  V] .  Accordingly,  our  received  signal  r(t)  is  defined  by 


r(r)  =  s0(i)  +  $-(t)  +  n(r) 


=  ^R^2Pkbk  (t)  wt  (t)c{t)cos(27Tfct)  (3.22) 

k~0 

+S  ^Ri^h(t  +  ti)^ijit  +  Ti)ci(t  +  Ti)cM2^fct  +  (Pi)  +  n{t). 

/= 1  ;'=0 

We  will  use  the  received  signal  to  analyze  the  performance  of  the  DS-CDMA  system  in  a 
Rayleigh-fading  channel  with  lognormal  shadowing. 


D.  SUMMARY 


In  Chapter  IE,  we  have  developed  a  mobile  radio  channel  model,  the  Rayleigh- 
lognormal  channel  model,  which  combines  both  large-scale  and  small-scale  propagation 
effects  into  a  single  model.  The  Rayleigh-lognormal  channel  model  incorporates  the 
large-scale  median  path  loss  as  predicted  by  the  extended  Hata  model  and  it  includes  the 
effects  of  lognormal  shadowing.  Small-scale  propagation  effects  are  built  into  our 
channel  model  using  a  slow-flat  Rayleigh  fading  model.  Combining  large-scale  and 
small-scale  propagation  models  in  this  manner  allows  us  to  more  accurately  characterize 
the  mobile  radio  channel  in  a  single  unified  model.  Accordingly,  our  model  will  allow 
for  more  complete  analysis  of  the  mobile  radio  channel  than  past  models,  which  typically 
incorporated  small-  or  large-scale  propagation  effects  but  not  both. 

We  have  also  developed  the  forward  channel  of  a  typical  DS-CDMA  cellular 

system  consisting  of  traffic  from  a  cell’s  base  station  to  mobile  users  in  the  cell.  We 
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applied  the  Rayleigh-lognormal  channel  model  to  the  signal  transmitted  from  the  base 
station  in  order  to  characterize  the  input  signal  to  a  mobile  user’s  receiver.  The  received 
signal  includes  the  traffic  intended  for  the  mobile  user,  traffic  for  other  users  in  the  cell 
(intracell  interference,)  traffic  for  users  in  adjacent  cells  (intercell  interference,)  and 
additive  noise.  We  will  use  this  received  signal  to  develop  a  Signal-to-Noise  plus 
Interference  Ratio  (SNIR)  and  for  Bit  Error  Rate  (BER)  analysis  in  Chapter  IV. 
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APPENDIX  III-A  LOGNORMAL  RANDOM  VARIABLE 


1.  Lognormal  Random  Variable  Derived  from  a  Gaussian  Random  Variable 


Let  be  a  Gaussian  random  variable  with  mean  and  standard  deviation 
<7^,  denoted  as  ~  JV{f*dB,crdB).  The  pdf  of  the  random  variable  is  given  by 


Pjr(*)  =  ' 


1 


dB 


-°°  <  X  <  °°. 


(3.23) 


V2 nadB 

We  can  transform  the  random  variable  from  units  of  decibels  to  ratio  form  using 


=  lOlogX  ->  X  =  1010  =  g(jC), 
where  g  (±¥)  is  the  transformation  function  and  its  derivative  is  given  by 

=  =— X  =  2X, 

K  dJC  10  10 


where 

.  InlO 

A  — - . 

10 

We  know  that  since  g  (JC)  is  a  monotonically  increasing  function,  the  pdf  of  the 
transformed  variable  can  be  developed  as  follows: 
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1 


Pxw=rap-*(x) 


Pjf(lOlogx) 


-(lOlogx-//^) 


Ax-JIncr 


2 adB 


dB 


1 


-A2(l01og  x-/^g)2 

e  U2(TdB 


-(Al0logx-^g)2 

e  2{XadB] 


1 


-((inio)logx-A/^g)2 

e  4^®) 


1 


-(hx-A^)2 
^  2(^<b) 


-(Inx-//X)2 
„  2<Jv 


(3-24) 


x  >  0, 


where 


Mx  —  A,fidB  and  <7X  — 


(3.25) 
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We  say  that  the  random  variable  X  has  a  lognormal  distribution  [23]  with  parameters 
fly.  and  <rx ,  and  is  represented  as  X~A(/ox,<Jx)  or,  alternatively,  in  terms  of  the 

original  parameters,  X  -  A(AjudB,AcrdB) . 


The  nth  moment  of  X  about  the  origin  is  given  by 


£{X',}  =  /'Ux+  2  =  exp 
Accordingly,  the  mean  of  X  is 


nMx 


ncri 


£{X}  =  /x+C?, 


the  second  moment  about  the  origin  of  X  is 

£{X2}  =  e2Mx+2a*, 


and  the  variance  of  X  is 


C yar{X}  =  e2Mx+axlecrx  -1 


Furthermore,  the  median  of  X  is 


Med(X)  =  eMx. 


2.  Transformation  Y  =  a /X  of  the  Lognormal  Random  Variable 


If  X  -  A(//x,crx) ,  we  can  define  a  transformation  from  X  to  a  new 
variable  Y  as  follows: 


a  >0 . 


(3.26) 


(3.27) 


(3.28) 


(3.29) 


(3.30) 


random 
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We  can  find  the  distribution  of  Y  as  in  Appendix  IHA.l,  where  our  transform 
function  is  defined  by 

£'(x)  =  gr¬ 

ille  pdf  of  Y  is  developed  as 


-(lna-!ny-/zx)2 

1 

e  2crx 

yV27rcrx 

■  o 

-(lny-lna+//x)2 

1 

■e  2cr* 

Jyf27T(Jx 

-flny-^y)2 

1 

p  2o$ 

yV2 7T<7y 

which  is  in  the  form  of  a  lognormal  random  variable  as  defined  by  (3.24).  Accordingly, 


our  transformed  random  variable  Y  ~  A  (jUy,aY)  is  also  lognormal  with 
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//Y=lna-//X  and  <7Y=<7x- 


(3.31) 


Furthermore,  if  we  want  to  represent  the  parameters  in  terms  of  decibels,  we  can 
convert  as 

lna  =  lnlOloga  =  A^, 

with 

cu  =  10  log  a. 

Finally, 

fjy  —  In  a  -  =  hex/ — —  /l  {cu  —  pidB ) . 
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IV.  DS-CDMA  PERFORMANCE  ANALYSIS 


In  Chapter  m,  we  developed  a  Rayleigh-lognormal  channel  model  and  described 
the  signal  received  by  a  mobile  user  on  the  forward  channel  of  a  typical  DS-CDMA 
cellular  system,  which  has  passed  through  the  Rayleigh-lognormal  channel  model.  In  this 
chapter  we  assess  the  performance  of  such  a  system  under  various  conditions,  which  can 
be  adjusted  in  our  model.  We  apply  Forward  Error  Correction  (FEC)  Coding  to  the 
forward  signal  in  the  form  of  convolutional  codes  to  obtain  the  coded  probability  of  bit 
error.  Due  to  the  complexity  of  our  channel  model,  we  use  statistical  modeling  of 
solutions  and  probability  distributions  in  cases  where  no  analytical  solutions  can  be 
(reasonably)  developed. 

In  order  to  assess  the  performance  on  the  forward  channel,  we  take  mobile  user 
one  in  the  center  cell  as  the  receiving  mobile  user  without  any  loss  of  generality  in  our 
final  solution.  We  assume  that  mobile  user  one  is  located  in  the  least  optimum  position 
within  the  cell,  which  is  in  any  comer  of  the  hexagon  as  shown  in  Figure  4.1. 


Figure  4.1  Mobile  User  One  in  the  Seven  Cell  Cluster. 
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The  distance  from  the  base  station  in  the  center  cell  to  mobile  user  one  is  d , 
while  the  distance  from  mobile  user  one  to  the  base  station  in  adjacent  cell  i  is/>. .  The 

length  of  each  side  of  the  hexagon  is  length  d  as  well,  which  implies  that  the  hexagon 
actually  comprises  six  isosceles  triangles.  Using  the  geometry  of  the  triangles,  as  shown 
in  Figure  4.2,  we  can  determine  the  distances  from  each  base  station  as  a  function  of  d. 

d 


Figure  4.2  Geometry  of  the  Cell  Hexagon. 


Accordingly,  the  distances  from  adjacent  cell  base  stations,  Z) ,  to  mobile  user 
one  for  the  worst-case  distance,  d,  from  the  base  station  in  the  center  cell  are  as  follows, 

d,  i  =  4,5 

B  =  2d,  i  =  3,6  (4.1) 

yjid,  i  =  1,2 
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These  distances,  Dl ,  are  used  in  computing  the  Hata  median  path  loss  for  signals 

transmitted  from  the  respective  base  stations,  which  contribute  to  the  co-channel 
interference  in  the  received  signal.  The  median  path  loss  of  the  information  carrying 
signal  from  the  center  cell’s  base  station  is  computed  using  the  distance,  d. 

We  are  now  prepared  to  proceed  with  the  performance  analysis  of  this  system 
model.  We  develop  the  Signal-to-Noise  plus  Interference  (SNIR)  along  with  the 

corresponding  uncoded  probability  of  error  3*e  in  Section  IV. A.  We  apply  convolution 
encoding  to  our  system  and  derive  a  (coded)  probability  of  error  ^  in  Section  IV.B.  In 
Section  IV.C,  we  develop  a  statistical  model  that  allows  us  to  analyze  Q  more 
completely,  which  we  accomplish  in  Section  IV.D. 

A.  PERFORMANCE  OF  THE  BASIC  SYSTEM 

In  order  for  mobile  user  one  to  extract  the  traffic  intended  for  him  from  his  base 
station,  he  must  process  the  received  signal  r[t)  as  defined  in  (3.22)  one  bit  at  a  time. 
Mobile  user  one  processes  each  bit  of  the  signal  as  shown  in  Figure  4.3. 

The  upper  branch  in  Figure  4.3  is  the  data  branch,  in  which  the  data  intended  for 
mobile  user  one  is  despread,  demodulated,  and  finally  integrated  over  the  bit  to  form  a 
decision  statistic,  Y.  The  decision  statistic  Y  is  used  to  determine  whether  the  processed 
bit  was  a  logical  1  or  0.  The  lower  branch  in  Figure  4.3  is  the  pilot  recovery  branch.  In 
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the  pilot  recovery  branch,  the  pilot  tone,  which  is  used  to  coherently  demodulate  the 
received  signal,  is  extracted  from  the  received  signal.  As  shown  in  Figure  4.3,  we  use  a 


Figure  4.3  Processing  the  Received  Signal  r(t). 


very  narrow  filter  around  the  carrier  frequency  to  extract  the  pilot  tone  from  the  signals 
that  remain  spread  over  the  band. 

We  examine  the  components  of  Figure  4.3  in  detail  in  Sections  IV.A.l-rV.A.2.  In 
Section  IV.A.3  and  IV.A.4,  we  develop  the  SNIR  and  Q  ,  respectively. 

1.  The  Demodulated  Signal  y2  (f ) 

In  this  section,  we  develop  the  demodulated  signal  y2(r).  Let  us  start  with  the  data 

branch  of  Figure  4.3.  We  apply  the  cell’s  PN  sequence,  which  is  synchronized  with  the 
base  station,  followed  by  the  Walsh  sequence  (in  the  form  of  time  signals)  in  order  to 
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despread  and  remove  the  orthogonal  Walsh  cover  from  the  desired  information  bit  as 


follows: 


Vi  (t)  =  r(t)c{t)w1(t) 

=W0+?M+"(#WwiM 

=  J0  (t)c(t)  Wj  (?)  +  S‘(0C(0  W1  (0  +  n  Wc  (O' W1 W  (4-2) 

' - y - /  ' - v - t  ' - - - / 

h+Yi  Si  "i 

=  7,  (?)  +  ft  (?)  +  (?)  +  (?) , 

where  we  have  broken  ft  (?)  into  four  components  as  defined  below. 

The  first  two  components  result  from  the  processing  of  the  forward  signal  from 
the  base  station  in  the  center  cell  as  follows: 

7i(0  +  ^W  =  JoWcWwiW 

=  A(0W*  (r)c(r)C0S(2?r/c?)jc(?)w1  (?) 

=  (R-J2PA  (?)  w,  (?)c(?)cos  (2^/c?))  c(?)  w,  (?) 


X Rj2Pkbk  (?)  W*  (?) c(?) COS (2^r/c?)  c (?) W,  (?) . 


The  desired  information  bit,  bx  (?) ,  is  contained  in  the  signal  Ix  (?) : 

I 1  (?)  =  (Ryl2PA  (?)  Wj  (?)c(?)cos(2^/c?))c(?)  wt  (?) 
=  (?)cos(2^/c?). 


Intracell  interference  is  contained  in  the  term  ft : 
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c(Ow-.(0 


(4.4) 


Yi  (0 


A 


jr-i  _ 

E  ^V2^  W  w*  (0  c(t)cos[2n  fct) 

k= 0 


=  E^V2^*  Ww*  (0  w,(r)cos(2ff/cf). 


*=0 

**1 


The  intercell  interference  is  contained  in  the  term  gx  (?) : 


^  6  K,  -l  _  "N 

?i(0=  E  E^V^y(r  +  T0Wy(r  +  r«)c-(r  +  r/)COS(27r/c?  +  ^)  c(r)wi(f) 

I  i=l  j= o  j 


«  ir,-i 


(4.5) 


E  W  (r + T<- )  (* + T' ) ' wi  (0  c«  (r + T«  )c  (0cos(2?r//+ % )  ■ 

1=1  ;=0 


Finally,  the  thermal  noise  contribution  is  contained  in  the  term  nx  (?)  as  follows: 

n1(t)  =  n(t)c(t)w1(t).  (4.6) 

The  pilot  tone  for  each  base  station  is  transmitted  in  Walsh  channel  zero,  using 
w0  (?) .  Accordingly,  to  recover  the  pilot  tone  and  demodulate  the  information  signal  we 
despread  the  received  signal  and  split  off  a  copy  to  the  pilot  recovery  branch,  as  shown  in 
Figure  4.3.  The  unfiltered  pilot  signal,  p  (?) ,  can  be  described  by 

p(t)  =  r(t)c(t)\N0(t) 

=  {s0(t)  +  g(t)  +  n(t))c(t)\N0(t) 

=  s0(?)c(?)w0(?)  +g(t)c(t)w0  (?)  +  n(?)c(?)w0  (?)  (4-7) 

h+Yo  So  "o 

=  h  (0  +  Yo  (0  +  So  (r)  +  no  (0  » 

where  we  have  (once  again)  broken  the  unfiltered  pilot  into  four  components. 

The  first  two  components  are  a  result  of  the  forward  signal  in  the  center  cell  and 
are  defined  by 
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/oW+roW  =  5o(OcWwoW 

WcWcos(2tf/cf)jc(r)w0(f) 
=  ( Ryj2P0b0  (t )  w0  (r )  c(t) cos ( 2ft  fct ))  c{t)  w0  (t ) 


k=0 
fK-1 


j^R^2Pkbk  (t)\Nk  (t)c{t)  cos  (2  ft  fct)  c(t)w0  (r), 


*:=  1 


where  the  desired  pilot  tone,  cos  (2ft fct ) ,  is  contained  in  the  pilot  recovery  signal  I0(t) 


as  follows: 

70  (r)  =  (R^2P0b0  (t)  w0  (r)c(t)cos(2;r/ct))c(t)  w0  (t) 

=  R^2^  cos  (2  ft  fct), 

since  b0  (t)  =  1  for  all  t. 

The  y0  ( t )  pilot  recovery  term  contains  intracell  traffic  and  is  defined  by 

[K—l  "N 

E  (0  w,  (f )  <# ) cos  (2 ftfct)  c(t)  w0  (t) 

k=i  J 

=  E  (0  w*  (0  cos  (2?r/cr )  • 


(4.8) 


(4.9) 


jfc=l 


The  pilot  recovery  term  g0  (t)  contains  the  intercell  traffic  given  by 


&«=("£  Kf^b„  (t  +  T^W^t+T;)  Ci  ( t  +  T; )  COS  (2ft  fct  +  (pt ) 
i= 1  j= o 

c(r)w0(r)  (4.10) 

=  £  E Ri^ibn (*  +  ) wc/ (*  +  W*  +  ^)c(r) cos(2ft fct  +  <pi). 


6  Kr\ 

S  X 

1=1  ;=0 


The  pilot  recover  term  n0  (t)  contains  the  noise  as  follows: 


«oW  =  nWc(r)wo(0  =  nWcW- 


(4.11) 


97 


We  see  that  the  terms  Y0(t),  (?),  and  nx  (?)  inp(r)  are  still  spread  spectrum 

signals  as  compared  to  the  narrowband  component  70  (?) .  Accordingly,  when  we  filter 
p(t)  with  a  narrow  bandpass  filter  centered  at  the  carrier  frequency  fc  only  the  70  (?) 
term,  which  contains  the  desired  pilot  tone,  will  pass  through. 

We  demodulate  the  received  signal  by  applying  the  filtered  pilot  70  (?)  to  yi(0 


which  generates  the  demodulated  signal  y2  (?)  as  follows: 


(0  =  ft  (0  h  W  =  {h  (0  +  ft  (0  +  ft  (0  +  ft  (0)7o  (0 
=  Vo  rJo  ft^o  "*■  Vo- 


(4.12) 


The  desired  information  bit,  £>,(?)  is  now  contained  in  the  demodulated  signal  7,70 


defined  by 


Vo  -  (R-fiPA  (t)cos(2xfct))(Rj2ii C0S(2 */,<)) 

=  2  «!7w>i  (<)cos2  (2»v) 

=  2 S’Vwf-.W  i(l  +  COs(te/tl)) 

=  R2J%PA(>)  (l  +  COs(4ffV)). 


The  intracell  interference  is  contained  in  the  term  ^70 : 


ft70  = 


E7 lj2Pkbk  (?)  w*  (? )  w,  (?)  cos  ( 2;r/c?) 


fc=0 

**1 


V 

tf-1 


(t?a/2^COs(2^/c?)) 


=  E  (0  w*  (O' W,  (?)  cos2  (2nfct) 

*= 0 

=  E/?2V^^Ww*Wwi(f)(1+cos(4Vc0)- 


£-0 


(4.13) 


(4.14) 
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The  intercell  interference  terms  contained  in  the  demodulated  signal  are  contained  in  the 


term  $J0 : 


(  6  K;~  1 


i=l  j-0 


Sih  =  IS  2 Pvby  ( t  +  rf ) \Ntj  ( t  +  zi )  Wj  (r) cjt  +  zt)c (t)cos (2nfct  +  (pt ) 

R^2P0(cos(2nfct)) 

=  rJp0'L  ’S2/?1.^(r  +  r1>,(r  +  TI>1(t)ci(r  +  ri)c(r> 
cos  ( 2  7i  fct  +  <pt )  cos  ( In  fct ) 

=  *V^E  w&(f+^)w1(0ci(?  +  ri)c(f). 

i  [cos  ($?,. )  +  cos  (4^/cr  +  (pt )] 

=  RyfPo  E  E  ^ M  (r  +  )  w.y  +  ^ )  wi  (0  c«  (?  +  )c (0* 

[cos  (p. )  +  cos  (4;r/cr  +  (pi)\ 


6  a:,— i 

s  s 

i—1  j-0 

6  K- 1 

s  s 

/=1  j= 0 

6  AT,— 1 

s  s 

i=l  7=0 


The  noise  contained  in  the  demodulated  signal  is  found  in  the  nxl0 : 

Vo=(«WcWw1W)(^cos(2^)) 

=  R\f^n(t)c{t)\Nl  (t)cos(2nfct). 


(4.15) 


(4.16) 


Accordingly,  the  demodulated  signal  y2  (?)  as  defined  by  (4.12)  through  (4.16),  is 


sent  into  the  integrator  to  determine  the  decision  statistic,  Y. 


2.  The  Decision  Statistic  Y 

In  this  section  we  will  develop  the  decision  statistic  Y,  which  will  be  used  to 
determine  whether  the  bit  bx(t),  for  re  (0,T),  is  a  logical  zero  or  a  logical  one.  We 
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calculate  Y  at  time  T  by  sending  our  demodulated  signal  y2  (?)  through  an  integrator  as 
shown  in  Figure  4.3.  This  results  in 

T  ,r  rT  ,r  <4-17> 

J0  IJodt  +  l  YJodt  +  jo  ^I0dt  +  f\I0dt, 

where  we  integrate  y2  (?)  term  by  term.  Furthermore,  we  will  fix  the  Rayleigh  fading 
random  variable  at  R  =  r  and  the  lognormal  random  variable  representing  the  power 
received  at  Pk  =  pk  in  order  to  develop  our  decision  statistic.  This  means  that  our 

decision  statistic  will  be  conditioned  on  R  =  r  and  Pk  =  pk ,  denoted  by  Y |  .  We  will 

remove  the  conditioning  from  our  solution  later,  by  integrating  out  the  dependence  on  the 
random  variables.  Accordingly,  our  decision  statistic  conditioned  on  the  random 
variables  R  and  Pk ,  is  defined  by 


rL 

—  J0  hhdt  +J0  Yj0dt  +  dt  +Jq  n} IQdt 


(4.18) 


r>Pk 


We  will  develop  each  of  these  four  components  separately. 

We  begin  with  the  first  term,  Y ,  which  is  the  output  of  the  integrator  containing 


our  desired  information  bit.  We  define  Yby 
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(4.19) 


Y  =  />.<* 

=  J0T  r2JpMb i  (?)  [l  +  cos  (4^  /cr )]  * 

=  ^b,  }0T^  +  J0rcos(4^/cr)^r 

where  b*  e  {±1}  corresponds  with  the  time  function  bk  (t) ,  which  is  constant  over  the 
period  (0,r).  We  assume  fc  is  an  integer  multiple  of  1 JT .  Accordingly,  is  the 
desired  information  bit  in  the  decision  statistic  Y I 

1  r’Pk 

The  intracell  interference  is  completely  eliminated  from  our  decision  statistic  due 
to  the  orthogonality  property  of  the  Walsh  functions  as  follows: 

K—  1 

JoT Yxhdt  =  \l X  (0  wt  (O' wx (0 [l +  cos (4;r/c0] dt 

r*Pjt  fc=o 
k*l 


=  0, 


where  the  first  integral  is  zero  due  to  the  orthogonality  of  the  Walsh  functions,  and  the 
second  integral  is  zero  due  to  the  fact  we  have  further  assumed  that  fe  is  an  integer 
multiple  of  l/Tc ,  or  N/T.  Since  the  product  of  Walsh  functions  is  another  Walsh 
function,  we  can  break  the  interval  (0,r)  into  N  pieces  integrating  over  each  chip 
interval. 

The  intercell  interference  contribution  to  the  decision  statistic  is  contained  in  the 
term  C  ■ 
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(4.21) 


T  6  Ki- 1 

=lr4^^TRi4^h(t+T‘)vjAt+Ti)yNi{t)ci(t+ri)c{ty 

i'=l  ;=0 

[cos  ( <p. )  +  cos  (■ 4/r/cr  +  (p. )  ]  Jr 

_  r  6  JC,-1 

=  W ^0  J0  S  2  (r  +  r' )  WV  (?  +  ^ ) W1  (f )  c«  +  ^  )c (^ )  c°s (^.  >?r. 

;=i  ;= o 

The  additive  noise  contribution  to  the  decision  statistic  is  contained  in  the  term  r] : 

V  =  j\l0dt 

r.P  * 

=  Jo  r^/2p0n(r)c(r)w1(r)cos(2^’/ct)jt  (4.22) 

=  r^Jo  n(r)c(r)w1(t)cos(2^:/cr)<ir. 

We  can  combine  the  noise  and  intercell  interference  terms  in  our  decision  statistic 
into  a  single  term  £  : 


7L=y+^ 

=  Y  +  £ 


(4.23) 


Accordingly,  our  decision  statistic,  conditioned  on  the  random  variables 
R  =  r  and  Pk  =  pk ,  is  defined  in  (4.18)  through  (4.23).  As  defined,  Y\  is  not 

\r,pk 


particularly  helpful  in  further  developing  the  performance  analysis  for  our  system  model. 
In  order  to  move  forward  with  our  analysis,  we  will  use  the  Gaussian  Approximation 
from  [25]  and  [26]  to  model  our  decision  statistic.  As  the  name  implies,  we  model 
Y|  as  a  Gaussian  random  variable  The  mean  value  is  deterministic  (for  R  and  Pk 

fixed)  and  is  given  by 


£{^}  =  Y, 


(4.24) 


where  Y  is  defined  by  (4.19).  The  variance  of  our  decision  statistic  is  the  sum  of  the 
variances  of  our  co-channel  interference  term  and  noise  term,  which  are  assumed  to  be 
independent: 

Q?ar{2/}  =  Q7ar{C}  +  Q&rfy}  =  Wartf]  =  o\.  (4.25) 

Accordingly,  we  have  modeled  our  decision  statistic  as  a  Gaussian  random 
variable,  ~  JV{y,(T^  .  This  approximation  of  the  decision  statistic  is  based  on  the 

assumption  that  the  co-channel  interference  term  £  and  the  additive  noise  term  r/  can  be 
modeled  as  independent  zero-mean  Gaussian  random  variables,  as  in  [19].  We  will  use 
this  approximation  of  our  decision  statistic  to  develop  the  SISTER,  and  probability  of  error 

?• 


3.  Signal-to-Noise  Plus  Interference  Ratio 

In  this  section  we  will  develop  a  conditional  SNIR  for  the  DS-CDMA  forward 
channel  in  the  Rayleigh-lognormal  fading  channel.  We  will  not  remove  the  conditioning 

on  the  random  variables  R  =  r  and  Pk—pk  until  we  develop  the  probability  of  error  T?  in 

Section  IV.A.4.  We  define  the  SNIR  at  the  output  of  the  receiver  as  a  ratio  of  the  average 
power  in  the  message  component  of  the  signal,  to  the  average  power  of  the  noise  and 
interference  components,  both  taken  at  the  output  to  the  receiver  after  [27].  As  developed 
in  Section  IV.A.2,  our  decision  statistic  ^  is  made  up  of  the  message  component  of  the 

signal  Y  and  the  noise  plus  interference  term  £ .  Accordingly  the  conditional  SNIR  is 
defined  by 
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Since  our  co-channel  interference  contribution  £  is  modeled  as  a  zero-mean 
random  variable,  we  define  its  variance  by 

C%r{f}  =  2{r}  =  -%pt|'E{«,2}s{Pf},  (4.27) 

1=1  ;=0 

where  the  complete  derivation  of  this  term  can  be  found  in  Appendix  IV-A. 

Similarly,  the  additive  noise  contribution  to  the  decision  statistic  is  zero-mean  and 
its  variance  is  developed  as 

Q?ar{rj\  =  £{V}  =  £  Wx  {t)  COS  (2n  fct)dt^  j 

2r2Po[\ln(t)c{t)w,  (t)cos(2nfct)dt^ 

(X,)cos(2^/cX)^xj 

=  2 r2p0E | £  Jq  n  (t)  n  (X)  c (t)  c (X)  Wj  (t )  Wj  (X) COS  (2^r/cr )  COS  {in  fcX) dtdX J 
=  2r2p0jo  Jq  £{n(t)n(X)}c(r)c(X)w1(r)w1  (%)cos(27t  fct)cos{2n  fcX)dtdX, 

where  we  have  taken  the  expectation  operator  into  the  integral.  We  assume  n(t)  is  a 
wide-sense  stationary  white  noise  process.  Accordingly,  the  autocorrelation  of  the 
process  is  defined  after  [24]  by 

£{n(r)n(X)}  =  ^£(r-X).  (4.28) 

We  substitute  the  autocorrelation  function  of  the  white  noise  process  and  continue 
developing  the  variance  of  rj  as 
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07ar{n}  =  2r2p0fjTo^S(t-%)c(t)c  (X) W,  (?)w,  (X> 

cos  ( In  fct )  cos  (2 7tfc%)  dtdX 
=  lr2p0  ^  c1  (r)w=  (t)oos2  (2nfrt)  dt 

=  -'2PoN0JoTi[l  +  COs(4ar/tt)]<i( 


=  ^PqNq 
2 


/•T  *7 

Jo  dt  +  Jo  COS  (4;? /c?)  dt 


r2p0N0T 
"  2  ’ 

where  c2  (?)  =  1,  and  wf  (?)  =  1. 

Using  (4.27)  and  (4.29),  we  can  update  the  variance  of  the  decision 
defined  by  (4.25)  as  follows: 

(T?  = 

_  r2p0T2  ^ 


EE2K}2k}- 


=r2POT 


3iV  l=1  ;=0 

j  JL  fc1 

3N  ^ 


rzA,Nor 

2 


+ 

J--H  '  s£2K}2M+i^ 


Accordingly,  we  define  the  conditional  SNIR  by 


SNIR  = 


_  (7)2  - 


'r,Pk  a* 


a ; 


r2p0T 


6  if,— 1 


^EE2KMnb- 


r  Pi 
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No. 
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(4.29) 


statistic 


(4.30) 


(4.31) 
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By  way  of  comparison,  we  can  remove  the  fading  and  shadowing  from  our  model 
by  letting  pj  =  Prl  be  the  average  received  power  in  our  information  channel,  and 


Py  =  Pr4j  be  the  average  received  power  from  the  yth  channel  in  adjacent  cell  i.  The 
revised  SNIR  without  fading,  shadowing  or  the  Hata-model  incorporated  would  be 


SNR  = 


1  6  K~l  ]SJ 

— y  y  p .. + — - 

3  r,y  IT 


1 


(4-32) 


Nft 


1  6  K- 1  n 

_Lv  V  _ 

3w^^/»rJT2r/>„ 

This  result  corresponds  with  a  similar  analysis  by  Rappaport  on  the  reverse  channel  of  a 
CDMA  Spread  Spectrum  system  without  fading  or  shadowing  in  [19]. 


4.  Probability  of  Bit  Error  Q 

In  order  to  develop  the  unconditional  probability  of  bit  error,  we  return  to  the 
decision  statistic  ,  which  by  definition  is  conditioned  upon  random  variables 

R  =  r  and  Pk  =  pk.  When  the  decision  statistic  is  a  positive  value  the  threshold  detector 

will  decide  that  the  bit  is  positive  (or  a  logical  zero).  Similarly,  when  the  decision 
statistic  is  a  negative  value  the  threshold  detector  will  choose  the  negative  bit  (or  a  logical 
one).  Accordingly,  we  will  have  a  decision  error,  or  bit  error,  if  the  decision  statistic 
<  0,  when  b:  =  +1  was  actually  transmitted.  We  will  also  have  a  bit  error  if  the 
decision  statistic  >  0,  when  =  -1  was  the  transmitted  bit.  Since  the  chances  of  a 
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logical  zero  or  logical  one  being  transmitted  are  equally  likely,  our  probability  of  bit  error 
can  be  described  as 


=-z 

1  ~ 

z 

+-z 

e 

r<Pk  2  e 

2 

r.Pt.  bj=l 

(4.33) 


Accordingly,  we  can  determine  the  probability  of  error  given  that  b,  =  +1  was  the  bit 
transmitted  by 


r.Pk.b  1=1 


=  Pr 


b,-.<0] 

■pr[5U+«<0] 

=  Pr[r2s[M>T  +  4<0], 

Since  our  decision  variable  is  modeled  as  a  Gaussian  random  variable,  -  JV  ) , 
the  probability  that  ^|b  =1  <  0  is  the  area  under  the  left  tail  from  to  0  given  by 


9. 
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(4.34) 


—a? 
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e  2  das 


V  ’  *  J 


=Q(i/SMRU)- 


where  Q(;t)  is  the  Q -function  defined  by 

“  1  zf. 


(4.35) 
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We  find  the  conditional  probability  for  bx  =  -1  and  y\h  __i  >  0  in  a  similar 
manner.  In  this  case  the  integral  is  evaluated  from  0  to  °° ,  and  since  Y  =-Y 

b,=-l  b,=l 

the  result  for  Q  is  the  same  as  (4.34).  Accordingly,  the  total  probability  of  bit 
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We  see  that  the  probability  of  bit  error,  conditioned  only  on  Pl  =  px ,  is  a  function  of 

R2  =  r2 ,  which  is  a  chi-square  random  variable  with  two  degrees  of  freedom. 
Accordingly,  we  transform  the  conditional  SNIR  after  [24],  by  letting 


Wi  = 


R2P i 


3^2  £*{*’}*{'.}■ 

ML  j= 0 


No. 

IT 


R2 


(4.38) 


where 


Pi 


i=l  j- 0 


No 
2  T 


(4.39) 


Accordingly,  our  SNIR  y/2  (conditioned  on  Pl  =  px)  is  a  chi-square  random  variable 


with  two  degrees  of  freedom,  whose  mean  can  be  described  by 


(4.40) 


where  we  have  normalized  the  underlying  Rayleigh  random  variable,  such  that 
£{i?2]  =  1 .  We  revise  (4.37)  using  the  transform  variable,  yr2,  and  corresponding  chi- 

square  pdf  from  [24]  as  follows: 
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2  +  ^2  y 


(4.41) 


The  details  of  the  integration  in  (4.41)  can  be  found  in  Appendix  IV-B.l.  Thus,  we  have 
removed  the  conditioning  on  one  of  our  random  variables,  initially  conditioned  as  R-r , 
from  our  probability  of  bit  error. 
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To  remove  the  dependency  of  fPe  on  the  lognormal  random  variable  P1  =  px ,  we 


must  integrate  the  solution  over  the  range  of  the  pdf  of  pl  as  follows: 


where 


oo 

%  =  m  Pn(Pi)4Pi 

**  'Pl 

•oo 

jhL 

-I _ ! _ f _ ! _ 

2  2^271  Ac dB  o  XyjOP'x  + 1 


Lh  (g)  £ 
Of  M^) 


y  LH(d)  |  N0 
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(4.42) 


(4.43) 


and  where  we  have  normalized 


Z{i?;2}  =  f{£2}=l. 


We  also  have  taken 


£'{l/Zi}  =  ZjZ,}  =  exp|:^^j  as  developed  in  Appendix  HI-A  where  =  ~k\idB  =0. 
Furthermore,  we  have  introduced  a  new  variable  Eb  =  <^PtT /LH(d) ,  which  represents 


a  baseline  received  bit  energy  without  the  effects  of  fading  or  shadowing.  The  details  of 
the  conversion  of  (4.42)  and  (4.43)  into  the  original  model  terms  can  be  found  in 
Appendix  IV-B.2. 

We  have  developed  the  probability  of  bit  error  for  the  uncoded  system  given  in 
(4.42)  and  (4.43).  The  performance  of  the  uncoded  cellular  system  under  normal 

operating  conditions  is  quite  poor  (Q- 1/2).  Accordingly,  for  our  performance  analysis 

to  be  meaningful,  we  must  add  forward  eiror  correction  (FEC)  coding  to  the  system, 
which  we  will  accomplish  in  the  next  section. 
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B.  CONVOLUTIONAL  ENCODING  OF  DS-CDMA 


In  order  to  improve  the  performance  of  the  DS-CDMA  cellular  system  operating 
in  the  Rayleigh-lognormal  channel,  we  institute  forward  error  correction  in  the  form  of 
convolutional  codes.  We  use  an  (a,/S)  encoder,  which  implies  that  for  every 

/£  information  bits,  n  coded  bits  are  produced  by  the  encoder  and  that  the  code  rate  Rcc 
is  given  by  Rcc  =/£/«.  The  coded  bits  are  spread,  orthogonally  covered,  BPSK 

modulated,  and  transmitted  by  the  base  stations,  in  the  same  manner  as  the  uncoded  bits. 
In  order  to  preserve  the  bit  rate  of  the  system,  the  bit  duration  for  the  coded  system  is 
reduced  as  Tcc  =  7RCC  =  T  £j/i . 

The  convolutional  encoder  has  a  constraint  length  of  v,  which  means  that  a  single 
information  bit  can  affect  the  output  of  the  encoder  for  a  maximum  of  v  shifts,  as 
defined  in  [27].  Accordingly,  at  least  one  of 'the  /S  shift  registers  that  make  up  the 
encoder  has  a  length  of  v  - 1 . 

1.  Viterbi  Branch  and  Path  Metrics 

To  decode  the  information  bit  stream,  we  use  the  Viterbi  Algorithm  with  soft- 
decision  decoding.  The  output  of  our  demodulator,  Y  from  Figure  4.3,  will  be  the  input 
to  the  decoder.  The  decoder  examines  sequences  of  coded-demodulated  bits,  which  we 
denote  to  determine  the  most  likely  information  bit  stream  for  mobile  user  1.  The 

decoding  decisions  are  based  upon  metrics  developed  from  the  input  sequences.  We  will 
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analyze  the  performance  of  the  coded  system  in  a  manner  similar  to  that  of  Proakis  [24], 
in  which  the  Viterbi  algorithm  was  applied  to  an  additive  white  Gaussian  noise  channel 
with  soft-decision  decoding. 

In  order  to  simplify  the  performance  analysis,  we  can  assume  that  the  all  zero 
sequence  (or  \  (?)  =  1  for  all  ?)  was  transmitted  [24].  Accordingly,  we  will  redefine  our 


demodulator  output  from  (4.23)  to  accommodate  the  analysis  of  coded  sequences  as 
follows: 


>Pk,/n 


(4.44) 


where  the  subscript  /  specifies  the  branch  (or  information  bit)  in  the  trellis  of  the 
decoder,  while  the  m  specifies  the  position  of  the  coded  bit  with  in  the  /th  branch  such 
that  m  -  \,2,...n ,  and  where 


(4.45) 


and 


£C  _  rc  c 

b  fa  b  jm  1 1  jtrt  ’ 


6  Kr  1 


=  0SE^^(^  +  ^)w,(?  +  r>1(?)c,.(?  +  T,Oc(?)cos(^>/?+  (4.46) 
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with 


/(m-l)  = 


/(/rz-l), 

(/  —  l)  a,  m—  1  and  /  ^  1 . 
0,  /w  =  1  and  /  =  1 
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Statistically,  the  decoder  output  which  is  conditioned  upon 

R  -  r?n  and  Rk  =  Pk,pn » is  approximated  by  a  Gaussian  random  variable,  just  like  Y\rpt 


in  the  uncoded  system.  Its  moments  are  adapted  from  the  uncoded  case  and  are  defined 
by 


(4.47) 


and 


_  PO  ,j/7i  ~  cc 

~  3N 

=  r^Po^Tcc 
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xX^KHnl 

r=l  j-0 
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2 

N« 


(4.48) 


The  Viterbi  soft-decision  decoder  branch  and  path  metrics  are  developed  using 
branches  per  path  through  the  trellis,  and  n  coded  bits  per  branch.  For  simplicity,  we 
view  each  coded  bit,  as  estimated  by  the  eth  path  through  the  trellis,  as  a  logically  coded 

bit,  denoted  e  {0,1} ,  rather  than  a  voltage,  bcJ‘!  e  {±l},  such  that,  bcJ‘]  =  1  -  2 . 

The  Viterbi  algorithm  branch  metrics  in  each  path  /  for  branch  /  are  calculated 
after  [24]  as 


/'/)=X<vM‘£’)-  <4-49> 

nt-1 


These  are  summed  over  PS  branches  to  form  path  metrics 
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(4.50) 


=EE*>(1_2C2)- 

/=1  «=1 

For  example,  if  we  let  the  all-zero  path  through  the  trellis  (which  is  the  correct  path)  be 
denoted  as  /  =  0 ,  then  =  0,  for  all  jm ,  and 


/“I 

c4-5 1) 

/=1  /*=1 
yy 

/ =1  /W=l 

Furthermore,  if  we  let  any  competing  path  through  the  trellis,  which  decodes  as  a 
combination  of  zeros  and  ones,  be  denoted  as  /  =  1 ,  then  its  path  metric  can  be  described 
by 


CM'1' 

/"I 

=ee**.(i-2c2)- 

/=1  /w=l 


(4.52) 


For  any  incorrect  path  £  =  1 ,  we  know  that  C^2  =  1  in  a  finite  number  of  coded 

bits,  which  we  will  denote  as  d  bits.  Accordingly,  when  we  compare  the  coded  bits  in 
the  correct  path  with  the  coded  bits  in  the  incorrect  path,  they  are  different  in  only  the 

d  positions  in  which  =  1 .  We  can  use  this  fact  to  develop  the  probability  of  error  for 

any  path  through  the  trellis  that  is  a  distance  d  from  the  all-zero  path. 


2.  First-Event  Error  Probability 

The  first-event  error  probability  is  defined  as  the  probability  that  at  a  node  £B  in 
the  trellis,  one  of  the  competing  paths  merges  with  the  correct  path  for  the  first  time  and 
the  competing  path  has  a  path  metric  that  exceeds  the  all -zero  path  metric  [24].  When 
this  happens,  the  decoder  discards  the  correct  path  in  favor  of  the  incorrect  path,  which 
results  in  decoding  error,  and  consequently  bit  errors.  Using  the  path  metrics  defined  by 
(4.51)  and  (4.52),  we  can  compute  a  first-event  error  probability  for  a  path  /  =  1  that  is  a 
distance  d  from  the  correct  path  as  follows: 

<P2(d) [  p  =  Pr{CM(1)  > CM(0)} 

=  Pr jcM(1)  - CM(0)  >  oj 

=  Pr{-2i£^>0 

[  /=1  m=l 

=  Prj|>;<o}, 

where  the  new  index  £  runs  over  the  set  of  d  bits  in  which  the  two  paths  differ,  or 
f/i  ~  f°r  c%  =  1-  Accordingly,  the  set  {^}  are  independent  and  identically 
distributed  Gaussian  random  variables,  whose  sum 

£  =  (4-54) 

/=1 

is  also  a  Gaussian  random  variable  with  a  mean 


(4.53) 
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/=i  \LH(d)xe 

_  y  2  1 

M«0  6 

_  Vc^gf^cc  y  r2; 

4M 


(4.55) 


We  converted  the  (fixed)  terms  p0/  and  pu ,  back  to  the  original  random  variables,  using 


(3.14)  adapted  as  follows: 


_ 


(4.56) 


where  X,  =  x/  is  our  lognormal  random  variable,  X,  ~  A(0,/icrds).  As  shown  in 
Appendix  1II-A.2,  the  transformation  of  X(  =  1 /X/  results  in  another  lognormal  random 
variable,  where  Xe  ~  A(0,/l<7^) . 

The  variance  of  the  Gaussian  random  variable  ^  is  defined  by 


Q 7ar{%]  =  W}  =  Ea#; 


(4.57) 


=  or«- 


Accordingly,  we  can  modify  (4.53)  using  (4.54)  through  (4.57)  as  follows. 
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(4.58) 


=  Q 


v  J 


where  we  have  introduced  a  new  random  variable  ^  ■  which  is  the  sum  of  d 

multiplicative  chi-square(with  2  degrees  of  freedom)-lognormal  random  variables  given 
by 


(4.59) 


and  where 
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(4.60) 


Here  we  have  reverted  to  the  original  lognormal  random  variable  using  (4.48),  and  taken 
£{-K,2 )  =  1 ,  and  £  j-^-j  =  expp-f^j ;  Ec  =  {/£/ a)  Eb  is  the  coded  bit  energy. 


If  we  assume  that  we  can  find  the  pdf  pz^  (z^)  for  the  new  random  variable,  then 

we  can  remove  the  conditioning  of  the  first-event  error  probability  by  integrating  across 
its  pdf  to  obtain 


J^ML  PzM&m 


=  1°  PzM&j’ 


\aj 


(4.61) 


where  9^  (d)  is  the  probability  of  first-event  error  for  a  path  that  is  distance  d  from  the 

all-zero  path  and  that  merges  with  the  all-zero  path  at  node  SS .  For  a  given 
convolutional  code,  there  are  many  such  paths  of  differing  distances  that  could  merge 
with  the  all-zero  path  at  node  dt .  Accordingly,  we  can  calculate  an  upper  bound  on  our 
first-event  error  probability,  92 ,  by  summing  all  the  possible  paths  as  follows  [24], 


9 

J2 


s  2 

d=d^ 


A/^2  {d). 


(4.62) 
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where  Ad  is  the  total  number  of  paths  that  is  a  distance  d  from  the  all-zero  path  and  that 

mergers  with  the  all  zero  path  for  the  first  time  at  node  We  can  also  calculate  an 
upper  bound  on  the  probability  of  bit  error  9>e  using  the  total  number  of  information  bit 
errors  that  is  associated  with  selecting  a  path  of  distance  d  from  the  all  zero  path: 

t.At,  PAW-  <4-ra> 

For  any  particular  convolutional  code,  the  values  of  Ad  and  j3d  can  be 
calculated.  Accordingly,  we  have  developed  a  probability  of  error,  ^ ,  for  the  coded  DS- 
CDMA  cellular  system  in  the  Rayleigh-lognormal  channel,  contingent  upon  finding  the 
distribution  for  our  random  variable  Zd,  which  is  a  non-trivial  undertaking.  We  address 
the  difficulties  and  propose  a  solution  to  this  problem  in  the  next  section. 

C.  APPROXIMATING  THE  SUM  OF  MULTIPLICATIVE  CHI-SQUARE2- 
LOGNORMAL  RANDOM  VARIABLES 

In  Section  IV .B,  we  developed  an  upper  bound  on  the  probability  of  bit  error  t?e 
for  the  coded  system.  This  bound  depends  upon  all  possible  path  errors  92{, d ),  as 
defined  by  (4.61),  that  are  at  distance  d  from  the  correct  path.  In  order  to  calculate 
92(d)  for  a  specific  distance  d ,  we  must  first  find  the  probability  density  function  for 

the  sum  of  d  chi-square2-lognormal  random  variables  as  defined  by  (4.59),  where  the 
chi-square2-lognormal  notation  represents  a  multiplicative  chi-square  (with  2  degrees  of 
freedom)-lognormal  random  variable.  In  Section  IV.A,  we  developed  the  probability  of 
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bit  error,  9e ,  for  the  uncoded  case  by  integrating  out  its  dependency  on  the  chi-square 
and  the  lognormal  random  variables  separately.  The  fact  that  $,(d)  is  dependent  on 
the  sum  of  d  chi-square2-lognormal  complicates  the  situation  for  the  coded  system.  If 
we  let  M  =  R2X  represent  our  single  chi-square2-lognormal  random  variable,  then  we 
can  reasonably  find  the  pdf  ioiJt,  which  we  represent  as  pM  (/») .  The  random  variable 

Zj,  which  is  a  sum  of  d  chi-square2-lognormal  random  variables,  has  a  probability 
distribution  function  that  can  only  be  determined  by  convolving  pM  (m.)  with  itself  d 
times  as  follows: 

Vzd  M  =  PA  («i )  *  Vm,  («*2  )  *  -  *  Vm,  («* )  • 

v _  v  _ / 

Furthermore,  for  any  particular  convolutional  encoder  we  require  a  series  of  Qid)  for 
d  =  d^,  d^  + 1,  d^  +2, .. .  in  order  to  calculate  the  upper  bound  on  the  probability  of 
bit  error,  as  defined  in  (4.63).  If  we  truncate  the  series  by  using  only  the  first  five 
terms,  e.g., 

I  ^/W+4 
«  d=dfut 

we  will  need  to  perform  5 d^  +10  convolutions,  or  integrations,  to  determine  Q  for  a 

particular  operating  scenario.  In  our  analysis,  we  will  consider  convolutional  encoders 
with  a  code  rate  of  1/2  and  constraint  lengths  ranging  from  7  to  9.  Typically,  d^  for 

such  encoders  ranges  from  10  to  12.  Accordingly,  for  each  operating  scenario  that  we 
wish  to  consider,  we  would  need  to  perform  between  60  and  70  convolutions  in  order  to 
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calculate  the  upper  bound  on  the  probability  of  bit  error,  ^  for  that  set  of  data.  The 

computational  cost  in  performing  these  integrations  is  too  high  a  price  to  pay,  which 
explains  why  analysis  of  coded  DS-CDMA  systems  that  include  both  fading  and 
shadowing  models  has  not  previously  been  accomplished. 


Our  situation  would  be  greatly  improved  if  our  channel  model  only  considered  the 
case  of  Rayleigh  fading.  For  example,  we  can  modify  (4.58)  to  eliminate  the  lognormal 
shadowing,  which  results  in 


(4.64) 


where  \f/2d  is  the  sum  of  d  chi-square  random  variables,  \ffe ,  with  two  degrees  of 
freedom.  The  distribution  of  y/2d  is  chi-square  with  2d  degrees  of  freedom  [24]. 


Accordingly,  we  remove  the  dependence  on  the  random  variable  y/2d  by  integrating 


across  its  pdf ,  obtaining 


(d)=j?2(d)\ Lp~  (w2y)did2s 

-CO  Wld 

=  JQ 


121 


which  has  a  closed  form  solution  as  defined  in  [24]  for  L-fold  diversity  and  adapted  for 
L  =  d  by 


k= 0 


d-\  +  k\ 


('-4. 


(4.65) 


where 


and  where 


is  defined  by 
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i-l  j-0 


+  - 


a(L„(I>,)  2E, 
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For  example,  Figure  4.4  depicts  the  upper  bound  on  the  probability  of  bit  error  for 
a  DS-CDMA  in  a  Rayleigh  fading  channel  using  (4.65)  and  (4.63)  using  a  rate  Vi 
convolutional  encoder  with  constraint  length  v  =  8  as  compared  with  the  uncoded 
system.  The  extended  Hata  model  for  a  large  city  is  incorporated  with  parameters 
d  =  1  km;  fc  =  2  GHz;  =  30  m;  hmohile  =  1  m.  By  way  of  comparison,  Figure  4.4 
also  includes  the  probability  of  bit  error  for  the  coded  and  uncoded  system  in  an  AWGN 
channel  with  an  average  received  Signal-to-Noise  Ratio  (SNR)  per  bit  of  yb  =  Ej N0 . 
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The  average  received  SNR  per  bit  for  the  Rayleigh  fading  channel  is  developed  as 


follows: 


n=£ 


R^PJ 


N  „I»  (d) 


since  we  have  normalized  E {l?2}  =  1. 

As  shown  in  Figure  4.4,  the  DS-CDMA  system  can  achieve  a  probability  of  bit 
error  in  the  Rayleigh  fading  channel  of  10-6  using  EEC  (Rcc  =  l/2;v=8)  with  an  average 
SNR  of  approximately  9  dB  with  10  users  per  cell.  In  the  uncoded  case,  we  see  the 


Figure  4.4  Performance  of  a  DS-CDMA  System  in  a  Rayleigh  Fading  Channel. 


123 


interference  floor  for  the  Rayleigh  fading  channel  and  the  AWGN  channel  developing  at 
an  average  SNR  per  bit  of  20  dB.  Accordingly,  in  the  case  of  Rayleigh  fading  only,  we 
can  calculate  the  upper  bound  on  the  probability  of  bit  error  defined  by  (4.63)  using  the 
closed  form  solution  in  (4.65)  for  (V)  without  performing  any  convolutions. 


If  we  consider  the  case  of  lognormal  shadowing  only  (no  Rayleigh  fading),  we 
modify  (4.61)  as  follows: 


9*(J) 


f 

!>/ 

=  Q 

\ 

/=i 

\ 

'  a 

K 

J 

(4.66) 


=  Q 


a 

v  J 


where  a  is  defined  by  (4.60).  We  have  also  introduced  a  new  random  variable 

d 

~  &d  ~  'f2uxe  ,  which  is  a  sum  of  lognormal  random  variables  Xe  ~  A(0,A(JdB)  fixed 


at  1^=^.  In  order  to  remove  the  conditioning  on  <zd ,  we  must  detennine  its 
probability  distribution  function. 


Since  the  sum  of  lognormal  random  variables  appears  in  many  important 

communications,  radar,  and  antennae  problems,  considerable  effort  has  been  given  in 

finding  its  distribution  function  [30)-[33],  While  no  closed  form  solution  for  the 

distribution  function  has  been  found,  there  is  agreement  that  the  sum  of  lognormal 

random  variables  can  be  approximated  as  another  lognormal  random  variable. 

Accordingly,  in  order  to  model  the  sum  of  independent  lognormal  random  variables  as 

another  lognormal  random  variable,  we  must  determine  suitable  parameters  for  the 
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model.  There  have  been  many  different  approaches  to  finding  the  lognormal  model 
parameters,  such  as  moment  matching  or  cumulant  matching  [33].  For  example,  the 
cumulant  matching  approach  seeks  model  parameters,  for  which  the  Cumulative 
Distribution  Function  (CDF)  of  the  model  lognormal  random  variable  matches  the  CDF 
of  the  sum  of  lognormal  random  variables  (as  simulated  using  Monte  Carlo  methods.) 
Once  the  lognormal  model  parameters  are  determined  using  one  of  the  methods  outlined 
in  [33],  the  approximated  distribution  is  completely  specified.  Accordingly,  for  the  case 

of  strictly  lognormal  shadowing,  we  can  remove  the  conditioning  of  9*2  (d)  on  Z2d-a:d 

by  approximating  it  as  another  lognormal  random  variable.  Consequently,  we  can  solve 
for  the  upper  bound  on  the  probability  of  bit  error  using  (4.63)  for  lognormal  shadowing 
only. 

In  summary,  the  distribution  of  the  sum  of  d  chi-square  random  variables  with  2 
degrees  of  freedom  is  chi-square  with  2d  degrees  of  freedom.  Furthermore,  the 
distribution  of  the  sum  of  lognormal  random  variables  can  be  approximated  as  lognormal. 
In  light  of  these  two  facts,  we  will  approximate  our  sum  of  d  chi-square2-lognormal 
random  variables,  Zd,  as  a  multiplicative  chi-square(with  Id  degrees  of  freedom)- 
lognoimal  random  variable. 

1.  Defining  the  Model  Zd 

In  order  to  approximate  Zd  as  a  chi-square  2d  -lognormal  random  variable,  we 
first  represent  our  model  Zd  in  terms  of  its  two  factors  defined  by 
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(4.67) 


where  ^  is  a  chi-square  random  variable  with  2d  degrees  of  freedom,  whose 
distribution  is  completely  specified  based  on  the  underlying  Rayleigh  random  variable, 
i.e.,  Vu  =  ^R2  •  Accordingly,  the  moments  of  are  defined  by 

£{iy = ^ 

£{]&}  =  (4.68) 

and  C yar{yi,}  =  d, 

where  we  have  normalized  £{i?2}  =  1.  Furthermore,  Xz  ~  A(juz,<jz)  is  the  lognormal 

component  of  Zd ,  where  the  model  parameters  pLz  and  az  are  yet  to  be  developed.  The 
moments  of  Xz  are  given  in  Appendix  III-A  by  (3.27)  through  (3.29). 

We  develop  the  moments  for  our  model  by  assuming  that  and  Xz ,  from 
which  we  obtain 


44  -  = e{3ue{x,} 


=  d  exp 


p2  > 


Az  +■ 


(4.69) 


The  second  moment  of  Zd  can  be  developed  from  the  second  moments  of  y2d  and  Xz 
as  follows: 


441  =  £{&*}}  =  £{]&}£{*J} 

1  J  "*  (4.70) 

=  [d  +  d1 )  exp{lpz  +  2crJ. ) . 

Finally,  the  variance  of  Zd  can  be  developed  from  its  first  and  second  moments  defined 
by 
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Q?ar{Zd]  =  de2Mz+^z  (1  +  d)e°z  -  d 


(4.71) 


The  remaining  step  in  approximating  Zd  as  chi-square  2d  -lognormal  is  in 

determining  the  model  parameters  pLz  and  a\  for  Xz.  We  develop  these  two  parameters 

by  relating  the  mean  (4.69)  and  variance  (4.71)  of  our  model  Zd  with  a  scaled  version 

of  mean  and  variance  of  the  original  chi-square2-lognormal  random  variable,  R2X . 

To  establish  a  value  for  juz,  which  is  related  to  the  original  chi-square2- 

lognormal  channel  conditions,  we  observe  in  simulations  that  the  mean  of  Zd  is 
proportional  to  d  times  the  mean  of  the  original  chi-square2-lognormal.  Accordingly, 
we  equate  the  mean  of  our  model  with  d  times  the  mean  of  R2X  scaled  by  a  factor  of 
as  follows: 

£{^}  =  ^E{R!x}  =  g fy^exp 

where  <y\  =  A2crjB  is  the  parameter  from  our  original  lognormal  shadowing  random 
variable,  X  ~  A(0,  crx ) .  By  applying  (4.69),  we  solve  for 

(cr2-<xj.)  .  . 

2  -+ln(g-.)-  <4-72) 

In  order  to  define  the  parameter  az,  we  observe  that  the  variance  of  Zd  is 
proportional  to  d  times  the  variance  of  the  original  chi-square2-lognormal. 
Accordingly,  we  equate  the  variance  of  our  model  with  d  times  the  variance  of  R2X 
scaled  by  a  factor  of  gp2  as  follows: 
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Qfcr{^}  =  %'1dQ?ar{R2x)=$'2decr* 


2e°x  -1 


By  applying  (4.71),  we  solve  for 


—  In(l  + 


(4.73) 


The  scaling  factors  ^  and  gp2  of  the  mean  and  variance,  respectively,  vary  with 
oiB  and  are  chosen  such  that  the  CDF  of  our  model  Zd  best  matches  the  CDF  of  Zd  for 
d  =  8, 9,  ...,16.  Table  4.1  depicts  the  values  of  and  for  each  adB  as  determined 


by  simulation. 


°dB 

^2 

2 

l 

.9875 

3 

l 

.9625 

4 

l 

.9125 

5 

l 

.8125 

6 

.9375 

.4875 

7 

.9125 

.3375 

8 

.8625 

.2 

9 

.8 

.1 

Table  4.1  Values  of  gf  and  gp2  for  Zd . 


With  the  model  parameters  completely  specified,  we  can  now  develop  the  pdf  of 
Zd  using  the  marginal  densities  of  its  components;  namely  yd  and  Xz .  We  return  to 
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the  transformation  defined  by  (4.67)  and  fix  Xz  =  x ,  which  implies  that  Zd  -  x\f2d. 
The  conditional  pdf  of  Zd  is  given  by 


We  remove  the  conditioning  of  Zd  on  Xz  =  x  as  follows: 


where  jiz  and  o2z  are  defined  by  (4.72)  and  (4.73),  respectively. 

We  have  developed  a  distribution  for  the  chi-square  'Id  -lognormal  random 
variable  ~Zd ,  which  approximates  the  distribution  of  a  sum  of  d  chi-square2-lognormal 

random  variables.  In  the  next  section  we  will  look  at  specific  examples  of  the  model  and 
compare  performance  results  using  our  model  with  the  results  as  predicted  by  Monte 


Carlo  simulations. 


2. 


Testing  the  Model  Zd 


In  this  section  we  will  explore  the  usefulness  of  our  model  as  it  applies  to  bit  error 
analysis  of  DS-CDMA  with  FEC.  As  an  example,  we  will  consider  the  case  of  oAB  =  7 

with  rate  Vz  convolution  encoding  with  v  =  8  and  approximate  the  probability  of  bit  error 
using  the  first  five  terms  of  the  union  bound  from  (4.63).  Accordingly,  we  will  examine 

our  model’s  ability  to  approximate  #  (d)  for  d  =  10  through  14,  since  these  are  the  key 

parameters  when  employing  rate  Vz  convolutional  encoding  with  constraint  lengths  7  and 

8. 

Figures  4.5  through  4.9  show  histograms  for  Zd  with  a dB  - 1  and  £{i?2]  =  l, 

which  were  formed  by  generating  d  sets  of  chi-square2-lognormal  distributed  data 
consisting  of  400,000  independent  samples  per  set  and  summing  over  the  d  sets.  We 
have  over-plotted  the  histogram  with  the  pdf  of  our  model,  Zd . 


Figure  4.5  Histogram  of  ZlQ  and  the  PDF  for  Zl<1  for  =7. 
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Figure  49  Histogram  of  4  and  the  PDF  for  2^  for  adB  =7. 


We  know  that  the  distribution  of  Zd  is  not  chi-square  2d  -lognormal  as 

approximated  by  Zd .  The  real  question  is  whether  Zd  can  well  approximate  Zd  in 

calculating  the  first-event  error  probability  as  defined  (4.61).  In  order  to  determine  this 
answer,  we  evaluate  the  integral  in  (4.58)  by  simulation  and  compare  with  the  result 
predicted  by  our  model.  We  simulate  the  integral  by  generating  d  independent  samples 
from  the  chi-square2-lognormal  distribution  and  summing  them  to  form  one  realization 
for  zd  and  consequently  one  realization,  pl ,  for  3^  (d) .  We  repeat  this  process  100,000 

times  and  form  our  point  estimate,  p ,  for  {d')  as  follows: 


_  j  I05 

p"io?§A' 


(4.75) 


Figures  4.10  through  4.14  show  the  point  estimates  for  the  simulated  %(d)  as  compared 
with  its  approximation  using  our  model  for  the  case  of  a ^  =  7  with  10  users  per 
adjacent  cell  and  d  =  10  through  14.  In  this  case,  our  model  well  approximates  the 
simulated  result  for  3i(d). 
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Figure  4.13  First  Event  Error  Probability  ?2  (13)  foradB  =  7 
with  10  users  per  cell. 


Figure  4.14  First  Event  Error  Probability  92  (14)  for  adB  =  7 
with  10  users  per  cell. 


We  incorporate  these  results  into  (4.63)  to  determine  the  probability  of  bit  error 
using  the  first  five  terms  in  the  union  bound.  For  a  rate  Vi  convolutional  encoder  with 
constraint  length  of  v  =  8 ,  we  obtain  a  modeled  and  simulated  probability  of  bit  error  as 
shown  in  Figure  4.15.  The  average  received  SNR  per  bit  for  the  Rayleigh-lognormal 
channel  in  Figure  4.15  is 


n  =  £ 


El 

Nn 


(4.76) 
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By  selecting  a  constraint  length  of  v  =  8  as  an  example,  we  incorporate  the 
modeled  and  simulated  results  of  ?2(d)  for  d  =  10  through  14  since  d^  =10  and 
/3l0  =  2,  Pn  =  22,  =  60,  =  148,  fil4  =  340 .  In  the  case  of  10  users  per  cell,  small 

differences  between  the  modeled  and  simulated  ?2(d),  as  shown  in  Figures  4.10  through 
4.15,  grow  larger  when  their  weighted  sums  are  applied  in  the  union  bound  (4.63), 
although  the  differences  are  still  quite  small  (on  the  order  of  10~4).  Nevertheless,  these 
differences  do  not  outweigh  the  usefulness  of  our  model  in  analyzing  the  performance  of 


Figure  4.15  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  [odB  =  7)  using  a  Rate  Vi  Convolution  Encoder  with  v  =  8 . 
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DS-CDMA  cellular  systems.  In  particular,  the  model  closely  matches  the  simulated 
result  in  the  range  from  10  to  15  dB,  which  is  a  practical  range  in  which  DS-CDMA 
cellular  systems  operate.  Furthermore,  the  differences  in  the  modeled  and  simulated 
results  are  most  noteworthy  when  the  interference  floor  is  in  a  range  between  10”3  and 
10-7,  where  even  the  slightest  differences  are  magnified  due  to  the  logarithmic 
presentation  of  results. 

We  should  note  that  both  the  modeled  and  simulated  results  predict  that  the  coded 
system  has  inferior  performance  to  the  uncoded  system  for  low  SNR  per  bit  (less  than  5 
dB).  This  misleading  prediction  is  a  result  of  using  the  union  bound  to  upper  bound  the 
probability  of  bit  error. 

Appendix  IV-C  provides  graphical  results  similar  to  Figure  4.15  for 
cdB  =  2,  3, ...  9.  In  the  next  section  we  will  use  these  results  to  explore  the  performance 

analysis  of  a  DS-CDMA  cellular  system  with  FEC  operating  in  a  Rayleigh  fading  and 
lognormal  shadowing  environment. 

D.  BIT  ERROR  ANALYSIS  OF  DS-CDMA  WITH  FEC 

In  Section  IV.C,  we  developed  a  model  which  approximates  the  distribution  for  a 
sum  of  multiplicative  chi-square(with  2  degrees  of  freedom)-lognormal  random 
variables.  In  this  section  we  use  that  approximation  to  analyze  the  performance  of  DS- 
CDMA  with  FEC.  For  example,  using  the  same  parameters  from  the  previous  section, 
we  can  compare  the  probability  of  bit  error  for  AWGN,  Rayleigh  fading,  and  Rayleigh- 
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lognormal  channels.  Figure  4.16  shows  these  results  for  10  users  per  adjacent  cell,  in 
which  we  vary  the  lognormal  shadowing  parameter  from  2  to  9. 


SNR  per  bit  (dB) 

Figure  4.16  Probability  of  Bit  Error  for  DS-CDMA  with  FEC  in  Various  Channel 
Conditions  with  10  users  per  cell  (Rcc  =  1/2  and  v  =  8  ) . 

As  depicted  in  Figure  4.16,  when  lognormal  shadowing  is  added  to  a  pure  Rayleigh 
fading  Channel,  the  probability  of  bit  error  drifts  away  from  that  of  the  pure  Rayleigh 
channel  as  we  increase  the  lognormal  shadowing  parameter  from  2  to  9. 

We  can  also  consider  performance  of  the  DS-CDMA  system  with  FEC  in  terms  of 
the  number  of  users  per  cell  that  can  be  supported  at  a  given  SNR  per  bit.  If  we  fix  the 
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average  SNR  per  bit  at  15  dB,  we  can  predict  the  probability  of  bit  error  for  a  range  of 
users  and  for  a dB  =  2, 3,...,  9  as  shown  in  Figure  4.17.  The  number  of  users  per  cell 


Figure  4.17  Probability  of  Bit  Error  for  DS-CDMA  with  FEC  in  Various  Channel 
Conditions  with  the  Average  SNR  =  15  dB  (Rcc  =  1/2  and  v  =  8  ) . 


determines  the  amount  of  intercell  interference  received  by  a  mobile  user.  As  the 
intercell  interference  increases,  so  does  the  probability  of  bit  error.  In  order  to  achieve  a 
desired  probability  of  bit  error,  system  designers  could  limit  the  number  of  users  based  on 
the  channel  conditions  as  indicated  in  Figure  4.17.  For  example,  a  DS-CDMA  system 
with  FEC  operating  in  a  channel  with  Rayleigh  fading  and  lognormal  shadowing  with  a 
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fading  parameter  of  <7^  =  7  would  be  limited  to  1 1  users  (active  channels)  in  order  to 
maintain  a  probability  of  bit  error  of  10"3  or  less. 

We  can  further  limit  the  amount  of  intercell  interference  by  sectoring  the  cells 
into  3  or  6  sectors  of  120°  or  60°  sectors,  respectively.  This  is  accomplished  by  replacing 
the  omni-directional  antenna  at  the  base  station  with  several  directional  antennas,  which 
radiate  within  a  desired  sector.  Accordingly,  if  we  assume  that  the  users  are  evenly 
distributed  in  the  adjacent  cells,  the  amount  of  intercell  interference  is  reduced  to  1/3  or 


0  10  20  30  40  50  60  70  80  90  100 

Users  per  cel! 


Figure  4.18  Probability  of  Bit  Error  for  DS-CDMA  using  Sectoring  for 
adB  =  7  with  an  SNR  per  Bit  of  15  dB  (Rcc  =  1/2  and  v  =  8  ) . 
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1/6  of  its  original  value.  By  reducing  the  intercell  interference  with  sectoring,  the 
performance  of  our  DS-CDMA  cellular  system  operating  in  a  channel  with  Rayleigh 
fading  and  lognormal  shadowing  is  greatly  improved.  Figure  4.18  shows  the  dramatic 
increase  in  performance  for  the  case  of  odB  - 1 . 


E.  SUMMARY 

In  this  chapter  we  developed  a  signal-to-noise  plus  interference  ratio  and  the 
probability  of  bit  error  for  a  DS-CDMA  cellular  system  operating  in  a  Rayleigh  fading 
channel  with  lognormal  shadowing  which  incorporates  the  extended  Hata  model.  We 
added  convolutional  encoding  for  forward  error  correction  and  developed  an  expression 
for  the  upper  bound  on  the  probability  of  bit  error,  which  depends  on  the  sum  of 
multiplicative  chi-square(with  two  degrees  of  freedom)-lognormal  random  variables.  In 
order  to  evaluate  the  expression  for  specific  operating  conditions,  we  approximated  the 
sum  of  d  chi-square-lognormal  random  variables  as  a  multiplicative  chi-square(with 
2 d  degrees  of  freedom)-lognormal  random  variable.  This  approximation  allowed  us  to 
analyze  the  performance  of  the  DS-CDMA  cellular  system  operating  in  a  Rayleigh- 
lognormal  channel  under  a  range  of  operating  conditions. 

In  analyzing  the  performance  of  our  DS-CDMA  system  in  this  chapter,  we  placed 
our  intended  mobile  user  in  the  worst-case  position:  at  the  edge  of  the  cell.  Under  this 
scenario,  two  of  the  interfering  base  stations  are  as  close  to  the  mobile  as  is  his/her  own 
base  station.  Each  active  channel  from  those  two  interfering  base  stations  is  received  as 


140 


interference  at  roughly  the  same  power  level  as  that  of  intended  traffic.  Accordingly,  our 
results  are  very  pessimistic.  In  the  next  Chapter,  we  provide  more  realistic  results  by 
assuming  that  the  mobile  user’s  position  is  randomly  distributed  in  the  cell. 
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APPENDIX  IV- A  VARIANCE  OF  THE  COCHANNEL  INTEREFERENCE 


In  developing  the  variance  for  the  cochannel  interference  term,  £ ,  we  separate 
the  term  according  to  the  individual  channel  contributions  as  follows: 


i= I  y=0 


where  Itj  is  the  contribution  to  the  co-channel  interference  term  from  the  individual  co¬ 
channel,  y  in  the  adjacent  cell,  i.  We  can  manipulate  ItJ  into  a  simpler  form  as  follows: 

la  =  J0T  (r  +  K('  +  )  wi  {* )  < ct  (t  +  rt  )c(t)  cos  (p, )  dt 

=  Rt  ^  cos  )  JQr  (t  +  x. )  wiy  (r  +  r, )  c,  (t  +  x, )  w,  (t)c(t)  dt  (4.78) 

ayit+Zi)  d^t) 

=  cos (n  )fo  atj  (t  +  t; )d1(t)  dt, 

where  (r)  =  fy,.  (r)wi;/.  (r)c,  (r)  and  <ij  (r)  =  w,  (f)c(r)  are  simply  other  PN  signals, 
whose  precise  form  is  of  little  significance  in  our  analysis. 

We  develop  the  interference  power  contribution  from  an  individual  intercell 
channel  from  the  second  moment  since  the  co-channel  interference  is  modeled  as  zero- 
mean.  Accordingly, 
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2 W}  =  2 1(^ cwfaOjS  (*  +  r‘M  M<*J  } 

=  2{^  cos2  (^.)(joTa&.  (f +  *i)4  (*)*)(£"#  (X  +  7.M  (*)<**)} 

=  ~  ~  ~  2  {IT  aij  {*  +  l  H(X  +  Ti )  di  (0  di  (X) ^X} 

=  -Lil_Lil  Jj  £■  l{aij(t  +  Ti)aij(X  +  Tl)}  dtdX 


P{t-X) 


P[t-K) 


£K}*W 


Jo  Jq  p2{t-X)dtdX, 


where  fi(u)  is  defined  by  (2.14)  to  be  the  autocorrelation  function  for  a  random  binary 
signal  and  used  here  to  represent  the  autocorrelation  function  of  a  PN  signal. 
Additionally,  we  determined  S' {cos  ^  }  =  1/2 ,  by  assuming  the  phase  discrepancy  (pt  to 

be  uniformly  distributed  between  (0,2;r) .  In  order  to  continue  our  analysis,  we  perform 
a  change  of  variables,  from  (r,A.)  to  ( u,v ) ,  defined  by 

u  =  t-X,  and  u  +  v  =  2t, 

v  =  t  +  X  ,  u-v  =  -2X, 


thus, 


r  =  i(M  +  v),  and  X  =  ^(v  -u). 


We  calculate  the  Jacobian  determinant  of  the  transformation  as  follows, 


lx  =  det 
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The  new  limits  of  integration  are  determined  to  be 


-T  <u<T  and  |n|<v<2T-ju 


Figure  4.19  Transformation  of  the  limits  of  integration,  (t,  X)  ->  (u,  v ) . 


Applying  the  change  of  variables,  from  (t,X)  to  (u,v) ,  we  find  that 
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where  we  took  advantage  of  the  fact  that  the  region  of  integration  is  symmetric  about  the 
v  axis  as  shown  in  Figure  4.19. 

We  finish  the  development  of  £{/?}  by  substituting  for  0(u)  as  defined  by 
(2.14)  as  follows: 


{2T_2u)du 
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(4.79) 


Accordingly,  we  can  define  and  develop  the  cochannel  interference  term  using  the 
contributions  of  the  individual  channels,  assuming  those  contributions  are  independent  of 
one  another: 


Q*r{f}  =  £{f2}  =  E 
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APPENDIX  IV-B  DEVELOPING  THE  UNCODED  PROBABILITY  OF  BIT 
ERROR 


1.  Removing  Dependency  on  the  Chi-square  R. V. 


In  order  to  remove  the  dependency  of  %  on  the  chi-square  random  variable,  , 

we  convert  the  Q-function  to  the  error  function  and  integrate  across  the  chi-square  pdf, 
obtaining 

=j;|  l-erf  j 

=  f  J7  Pr,  (VijdVi -|j0"erff ^  V„ 

V  ) 

=  erf  f 

V.  '  J 

where 

V  J- 


Since  y/2  is  a  chi-square  random  variable  with  mean  value  of  y/2  as  defined  in  (4.40),  we 
replace  its  pdf : 
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We  have  a  change  of  variables,  from  ys2X.o  x ,  where 
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which  results  in 
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(4.82) 


where  the  integral,  7  (x) ,  has  an  analytical  solution  as  defined  in  a  general  form  in  [28] 


by  equation  4.3.8  as 


+  .f  +  1*#  j.  («3) 

provided  that  Re{f2]  >  0,  and  Re{j?}>-2,  and  where  2Fl  {ava2,b{,z)  is  a 


hypergeometric  series,  also  called  the  Gaussian  hypergeometric  function. 
The  hypergeometric  series  is  defined  after  [29]  equation  9.100  as 
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+  (fy  +  l)(/bj  +  2)»1«2»3 

which  converges  inside  the  unit  circle,  |z[  <  1 ,  provided  ava1,bl  >  0 . 


(4.84) 
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In  our  case,  a=  1,  6  =  /=,  and»  =  l.  Accordingly,  our  integral,  l(x) 
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where  the  hypergeometric  series  converges,  since  we  know  our  that  our  SNIR  is  much 


less  than  1,  and  consequently, 


-Vi 


<  1 .  Furthermore,  the  form  of  the  hypergeometric 


series  in  (4.85),  constitutes  a  special  form  of  the  series,  which  can  be  represented  as  an 
elementary  function  as  defined  by  equation  9.121.1  in  [29]  as 

2Fl(-n,faf3\-z)  =  (  l  +  zf,  (4.86) 

where  in  our  case,  z  =  “~>  n  =  — “•  Applying  (4.86),  to  our  result  in  (4.85),  we  find 


that 


J(x)  =  - 


.2, 


V  2, 


(4.87) 


=  Vt. 

4  \2  +  ^: 


Substituting  into  (4.82),  with  I(x)  found  in  (4.87),  we  obtain  a  probability  of  bit 


error 
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2.  Removing  Dependency  on  the  Lognormal  R.V. 

In  order  to  remove  the  dependency  of  ^  on  the  lognormal  random  variable,  we 
a  =  a[p1)  into  its  original  form,  defined  by  (4.39),  as  follows: 
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The  lognormal  random  variables  p,  and  Ptj  revert  back  to  their  original  form 


using  (3.14)  and  (3.20)  respectively,  as  follows: 
1  1 
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where  Z  ~  A(0,/1<t^)  and  Eb  =  q^1*T / LH  (d) .  Accordingly,  the  unconditional 


probability  of  bit  error  is  defined  by 
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APPENDIX  IV-C 


MODELED  AND  SIMULATED  PROBABILITY  OF  BIT 
ERROR  FOR  THE  RAYLEIGH-LOGNORMAL  CHANNEL 
USING  120°  ANTENNA  SECTORING 


Figure  4.20  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( aiB  =  2)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 
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Figure  4.21  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (adB  =  3)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 
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Figure  4.22  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( crdB  =  4)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 
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SNR  per  bit  (dB) 

Figure  4.23  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( crdB  -  5)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 


Figure  4.24  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (crdB  =  6)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 
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Figure  4.25  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  =  7)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 


Figure  4.26  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (adB  =  8)  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 


157 


V.  INCORPORATING  USER  DISTRIBUTION  AND  POWER 
CONTROL  INTO  THE  FORWARD  CHANNEL  MODEL 

In  Chapter  IV,  we  analyzed  the  performance  of  the  forward  channel  in  a  DS- 
CDMA  cellular  system  by  assuming  that  our  intended  mobile  user  was  situated  in  the 
outer  comer  of  the  hexagonal  cell.  We  can  obtain  a  more  realistic  view  of  typical 
performance  of  the  forward  channel  by  assuming  that  the  users  are  randomly  distributed 
within  the  cell  according  to  a  specified  probability  distribution.  The  position  of  the 
intended  mobile  in  the  cell  is  represented  by  a  random  variable  from  the  specified 
distribution,  which  we  incorporate  into  a  revised  probability  of  bit  error  in  Section  V .A. 

In  Section  V.B  we  build  on  the  user  distribution  analysis  to  incorporate  power 
control  into  the  forward  channel  model.  That  is,  instead  of  transmitting  at  a  constant  rate 
of  power  in  each  user  channel  as  in  Section  V.A,  the  base  station  adjusts  the  transmit 
power  in  each  channel  to  ensure  the  power  received  by  all  mobile  users  is  equal  and 
adequate.  By  reducing  the  power  in  the  user  channels  to  only  that  which  is  required,  we 
aim  to  reduce  co-channel  interference  and  enhance  overall  system  performance. 

To  simplify  our  analysis,  we  will  replace  our  hexagonal  cells  with  overlapping 
circular  cells  as  depicted  in  Figure  5.1.  The  overlap  accounts  for  soft-handoff  of  mobile 
users  between  cells.  We  take  the  center  base  station  to  be  the  origin  of  a  polar  coordinate 
system.  Accordingly,  the  position  of  our  intended  mobile  user  in  the  center  cell  is 
represented  by  (*,0),  where  0<*<1  and  -7t  <  6  <  7t .  Based  on  the  underlying 
hexagonal  geometry  depicted  in  Figure  4.2,  we  know  that  the  distance  from  the  base 
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Figure  5.1  Circular  Seven-Cell  Cluster. 


station  in  the  center  cell,  to  each  of  the  adjacent  base  stations  is  V3  ,  since  we  normalized 
d  —  1.  The  adjacent  base  stations  are  separated  by  an  angle  of  7t/3  with  respect  to  the 
origin.  Using  the  geometry  of  the  cell  cluster,  we  determine  the  distance,  fZ) ,  from  each 
adjacent  base  station  to  the  mobile  user  in  the  center  cell  as  a  function  of  %  and  0,  as 
follows: 


The  development  of  5)  can  be  found  in  detail  in  Appendix  V-A. 


Since  we  aim  to  randomize  the  intended  mobile  user’s  position  within  the  cell,  we 
are  also  randomizing  the  distances  S)r  These  distances  are  used  in  calculating  the  path 

loss  for  the  intended  information  signal  and  the  co-channel  interference.  The 
randomization  of  the  distances  complicates  the  use  of  the  Hata  model  in  predicting  path 
loss.  Accordingly,  the  performance  analysis  in  this  chapter  will  use  the  average  path  loss 

Ln[d)  as  predicted  by  (3.5),  with  a  path  loss  exponent  of  n  =  4. 


We  modify  the  analysis  from  Chapter  IV  to  reflect  this  change.  Specifically,  if 
we  assume  that  the  fixed  position  of  the  intended  mobile  user  is  known  precisely  ,0f ) , 


we  revise  (4.60)  as  follows: 

3 N  £„(*!)  2 

exp(^)f  V  ,  N„ 

3  N  tr£c 2  Ee’ 

where  5)  is  defined  by  (5.1)  and  where  Ec  is  given  by 


(5.2) 


(5.3) 


Using  a  we  can  determine  the  first-event  error  probability  and  the  upper  bound  on  the 
probability  of  bit  error  using  (4.61)  and  (4.63),  respectively.  Figure  5.2  shows  the 
probability  of  bit  error  using  a  path  loss  exponent  of  n  =  4 ,  as  it  compares  with  the 
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Figure  5.2  Comparison  of  Probability  of  Bit  Error  for  DS-CDMA  in  a  Rayleigh- 
Lognormal  (adB  =  7)  Channel  using  120°  Sectoring  (Rcc  =  1/2  and  v=8) . 


probability  of  bit  error  using  the  extended  Hata  model  as  from  Chapter  IV  where 


\  =  d  =  1  km  and  where  the  worst-case  SNR  is  given  by 


R'-qJP,T  ' 
N  „£»WX. 


=  £{i!J}E{Vx}i  =  E{X}^. 


(5.4) 
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A. 


USER  DENSITY  IN  THE  CELL 


In  this  section  we  assume  that  our  users  are  physically  distributed  in  the  cell 
according  to  a  given  distribution,  rather  than  placing  our  intended  user  in  the  worst-case 
position  in  the  cell.  This  allows  us  to  assess  the  typical  performance  of  the  forward 
channel  in  our  DS-CDMA  cellular  system.  In  order  to  incorporate  user  distribution  into 
our  performance  analysis  of  the  forward  channel,  we  require  a  joint  probability  density 
function  in  terms  of  %  and  0  that  represents  a  user’s  position  within  the  cell. 

We  define  a  user  density,  (?(*•),  as  the  normalized  number  of  users  per  unit  of 

area  at  a  distance  %  from  the  base  station  [34].  Since  the  user  density  is  a  function  of 
(only)  distance,  k  ,  from  the  base  station  and  not  a  function  of  0 ,  we  have  assumed  that 
the  user  density  is  circularly  symmetric  about  the  origin.  Accordingly,  the  number  of 
normalized  users  in  a  differential  ring  of  width  cU  at  a  radius  %  is  given  by 

G(*)2TVuk. 

The  total  normalized  number  of  users  within  a  distance  %  from  the  base  station 
can  then  be  viewed  as  the  cumulative  distribution  function  for  an  individual  user’s 
distance,  k  ,  from  the  base  station  given  by 

K 

FK(%)  =  fG(s)2nsds,  s>0. 

0 

Accordingly,  the  probability  density  function  pK  (*)  for  a  user’s  distance  from 
the  base  station  is  given  by 
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(5.5) 


Pk  (*)  =  ^-^1  =  2mG(t),  *>  0. 

d*c 

Since  we  assume  circular  symmetry,  the  probability  density  function  pe  (0)  for 
the  angular  position  of  a  user  is  uniformly  distributed  as  follows: 

Pe(Q)=^~>  -7t<0<7t.  (5.6) 

zn 

By  assuming  independence  of  k,  and  0 ,  we  develop  a  joint  pdf  for  a  user’s 
position  within  the  cell: 

P*,e  (*>e)  =  P*  (*)Pe  (6)  =  *G(i).  (5.7) 

Consequently,  for  any  valid  pdf  pR  (<t) ,  we  can  determine  the  user  density  function 
G(*)  using  (5.5)  and  the  joint  pdf  for  a  user’s  position  within  the  cell  using  (5.7).  For 
example,  if  we  assume  the  user  density  in  the  cell  is  uniform  as  defined  by 

G(*)  =  -,  0<*<1,  (5.8) 

71 

we  can  determine  the  distribution  for  the  user’s  distance  from  the  base  station  using  (5.5). 
In  this  case,  the  distribution  is  linear: 

p«(*)  =  2*,  0<  *  <1.  (5.9) 

The  joint  pdf  for  a  user’s  position  in  the  cell  is  obtained  by  applying  (5.7)  as 
follows: 

=  0<*<1,  -n<e<n  (5.10) 

71 
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To  incorporate  this  result  into  our  performance  analysis,  we  modify  our 
expression  for  the  first-event  enror  probability  defined  by  (4.61)  and  (5.2)  (conditioned  on 
H  =  k  and  0  =  0)  as  follows: 


-n  0 


(z<rKe  (*$)dzdd>tdQ, 


where 


a  = 


(5.11) 


(5.12) 


In  view  of  the  fact  that  the  baseline  bit  energy  received  energy,  Ec ,  is  a  function  of  *,  it 
is  also  random  variable  defined  by 


(5.13) 


where  Ec  is  the  baseline  bit  energy  for  the  fixed  worst-case  position  given  by  (5.3). 


In  the  case  of  linear  user  distribution,  the  first-event  error  probability  is  developed 
using  equations  (5.10)  through  (5.13)  as  follows: 
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The  upper  bound  on  the  probability  of  bit  error  is  given  by  (4.63)  using  the  first-event 
error  probability  defined  in  (5.14).  Furthermore,  we  can  implement  antenna  sectoring  as 
in  Chapter  IV,  which  simply  reduces  the  number  of  sectored  users  in  cell  i  to  KjS, 

where  S  is  the  number  of  sectors. 

We  simulated  the  integral  in  (5.14)  for  the  sectored  system  using  100,000  Monte 
Carlo  trials.  Figure  5.3  depicts  performance  results  using  the  linear  user  distribution  and 
120°  antenna  sectoring  as  it  compares  with  our  worst-case  user  position  results  from 


Figure  5.3  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  F ading  and 
Lognormal  Shadowing  (o'dB  =7)  using  a  Linear  User  Distribution  with  120° 

Sectoring  and  FEC  (Rcc  =  1/2  and  v=8) . 
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Chapter  IV  using  the  SNR  per  bit  as  defined  by  (5.4).  As  shown  in  the  figure, 
accommodating  the  linear  user  distribution  into  the  analysis  indicates  better  performance 
(~  3  to  4  dB)  than  that  predicted  by  our  worst-case  scenario.  The  linear  distribution, 
however,  still  presumes  that  the  user’s  position  is  more  likely  to  be  toward  the  edge  of  the 
cell,  since  £{»}  =  2/3,  and  the  edge  of  the  cell  is  at  %  - 1.  Accordingly,  using  the  linear 

user  distribution  in  the  performance  analysis  of  the  forward  channel  is  more  realistic  than 
the  worst-case  scenario,  however,  it  is  not  an  overly  optimistic  view. 


B.  RANDOMIZING  TRANSMIT  POWER  ON  THE  FORWARD  TRAFFIC 

CHANNELS 

In  the  previous  section,  we  added  user  distribution  to  our  analysis  in  order  to 
better  assess  the  performance  of  the  forward  channel.  In  this  section,  we  will  take 
advantage  of  the  user  distribution  to  adjust  the  transmit  power  in  each  user  channel  in  a 
way  that  ensures  the  power  received  by  each  user  is  equal.  By  reducing  the  transmit 
power  in  the  user  channels  to  only  that  which  is  necessary,  we  intend  to  reduce  the  co¬ 
channel  interference  experienced  by  users  in  other  cells. 

The  goal  in  our  power-controlled  system  is  to  ensure  that  the  power  received  by 
all  mobile  users  is  at  the  target  power  level  P  by  adjusting  the  transmit  power  in  each 
channel,  k,  with  the  power  control  factor  qs£ .  We  can  relate  the  target  power  received  to 

the  baseline  transmit  power,  Pt ,  by  an  attenuation  factor  as  follows 

P=£,  (5.15) 

a 
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where  JP,PV  and  a  are  constant  for  all  users  in  all  cells.  When  the  mobile  user  receives 

the  information  signal,  the  actual  power  received  is  measured  and  reported  back  to  the 
base  station  for  adjustment.  At  the  point  the  mobile  user  measures  his  received  power 
and  reports  it  to  the  base  station,  it  is  no  longer  a  random  variable  as  described  by  (3.14), 
rather  it  is  a  realization  of  the  random  variable,  which  is  known  precisely  and  represented 
in  our  analysis  as  follows: 


p  _ 

‘  4K)*»  4 (dt)xk~  d-xt' 


(5.16) 


where  dk  is  the  user’s  actual  distance  to  the  base  station,  and  xk  is  the  shadowing 
experienced  by  the  user.  Accordingly,  in  order  for  the  power  received  by  all  users  to  be 
at  the  target  level/*,  the  base  station  must  adjust  the  power  factor  such  that 


dn  x 

a 


(5.17) 


With  perfect  power  control,  the  power  received  by  intended  mobile  is  simply 


(5.18) 


By  adjusting  the  power  factor  to  achieve  the  target  power  received  by  all  users,  we  have 
defeated  the  lognormal  shadowing  effect  on  the  information  signal. 

We  develop  the  power  received  by  our  mobile  user  from  the  forward  signal  in 
adjacent  cell  i  channel  j  from  (3.20)  as  follows: 


P  c £XJP 

“  L,(3>,)X,  L,(S>,)X,' 


(5.19) 
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where 


Xt  =  Lognormal  random  variable  A(0,A<JdB )  representing  shadowing  on 
the  forward  signal  from  base  station  i,  as  received  by  the  mobile  in  the 
center  cell. 

=  Random  varible  representing  the  distance  between  mobile  user  j  in 
adjacent  cell  i  and  the  base  station  in  cell  i.  The  probability  distribution 
function  is  also  pK  (*). 

Xij  =  Lognormal  random  variable  A{0,AadB )  representing  shadowing  on 
the  forward  signal  from  base  station  i,  as  received  by  mobile  user;  in 
cell  i. 

Accordingly,  the  power  received  from  the  adjacent  cells  is  a  random  variable  based  upon 
user  distribution  and  lognormal  shadowing.  When  we  incorporate  (5.18)  and  (5.19)  into 
our  performance  analysis  we  find  that  the  first-event  error  probability  is  defined  by 
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(5.20) 


where  £ jx,yj  =  £{X,}  =  £{l/X;}  =  exp|/.2a^/2j  and  where  we  fix  PTCC=EC  as 

defined  by  (5.3)  for  comparative  purposes.  Accordingly,  we  remove  the  dependency  on 
re  as  in  (4.65)  as  follows: 


dtA(d-l  +  k^ 

» SI  k 


k= 0 


(5.21) 


where 


since  £ 


is  defined  by 


Finally,  we  remove  the  dependency  on  the  random  variables  *  and  6,  which  are 
contained  in  the  terms  5),  as  follows: 


n  1  , 


$(</)  =  J J  9i(d)  bKe(>i,Q)d’(dQ 


-*  o  v  "1’e  y 

Vr  ^S/V-l  +  JlA 


(5.22) 
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Using  the  linear  user  distribution  from  Section  V.A,  we  define  the  first-event  error 
probability  by 


-nO  ^=ol  K  J 


(5.23) 


We  can  calculate  an  upper  bound  on  the  probability  of  bit  error  using  (4.63)  and 
(5.23).  Furthermore,  we  implement  antenna  sectoring  as  in  Chapter  IV.  Figure  5.4 
compares  the  probability  of  bit  error  for  a  cellular  system  using  forward  power  control 
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Figure  5.4  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (crdB  =  7)  using  Forward  Power  Control  and  Linear  User 

Distribution  with  120°  Sectoring  and  FEC  (Rcc  =  1/2  and  v=8) . 


with  the  fixed-power  results  from  Section  V. A.  The  average  received  SNR  per  bit  for  the 
Rayleigh-lognormal  channel  with  perfect  power  control  is  developed  as  follows: 


In  the  case  of  10  users  per  cell,  the  power-controlled  system  clearly  outperforms 
the  system  with  fixed  power(~5dB).  For  40  users  per  cell,  however,  the  power- 

controlled  system  outperforms  the  fixed-power  system  for  yb  <  12  dB,  at  which  point  the 
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interference  floor  begins  to  develop.  We  see  that  the  interference  floor  for  the  power- 
controlled  system  is  actually  higher  than  the  interference  floor  for  the  fixed-power 
system.  In  order  to  explain  this,  let  us  examine  the  isolated  co-channel  interference 
contribution  from  (5.20)  as  it  compares  with  its  equivalent  expression  in  (5.14)  for  the 
fixed-power  system.  In  order  to  ensure  that  the  received  power  for  all  users  is  equal  in 
the  power-controlled  system,  we  adjusted  the  power  in  the  forward  channel  to  overcome 
lognormal  shadowing  effects.  By  doing  so,  we  have  in  effect  increased  the  expected 

value  of  the  co-channel  interference  by  a  factor  of  f{X}  =  exp(A2a^/2)  over  the  fixed- 

power  system.  This  extra  factor  of  £{X}  is  a  result  of  randomizing  the  transmit  power 

in  the  forward  channel. 

In  spite  of  the  increased  co-channel  interference,  the  performance  of  the  power- 
controlled  system  is  far  superior  to  that  of  the  fixed-power  system  for  lognormal 
shadowing  environments  of  odB  <  6  with  120°  antenna  sectoring.  Under  these 

circumstances,  the  benefits  eliminating  the  lognormal  shadowing  from  the  intended 
signal  outweigh  the  increased  co-channel  interference.  Consequently,  we  can  achieve  a 
probability  of  bit  error  of  less  than  10~3  without  entering  the  interference-limited  region 
of  operation.  As  shown  in  Appendix  V-B,  the  power-controlled  system  can 
accommodate  as  many  as  100  users  per  cell  using  120°  sectoring  for  adB  <6.  For 
cdg  =  6 ,  the  power-controlled  system  can  accommodate  up  to  80  users  using  120 
sectoring  and  still  achieve  the  desired  probability  of  bit  error  of  10  3 . 
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In  designing  a  system  with  power  control  that  operates  in  more  intense 
lognormal  shadowing  environments  (odB  >  6) ,  we  must  ensure  that  our  system  is  not 
working  in  the  interference-limited  region.  For  example,  if  we  expect  our  system  to  work 
in  the  y  b  =  15  dB  range,  the  system  could  accommodate  up  to  roughly  30  users  per 
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Figure  5.5  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 

Lognormal  Shadowing  ( oiB  =7)  vs.  Users  per  Cell  using  Linear  User  Distribution 

with  120°  Sectoring  and  FEC  (Rtt  =  1/2  and  v=8) . 

cell  in  a  o ^  =  7  shadowing  environment  using  120°  sectoring.  Adding  more  users 
would  degrade  performance  due  to  operating  in  the  interference-limited  region  as 


indicated  in  Figure  5.5. 


If  a  greater  user  capacity  is  required,  we  could  increase  the  number  of  sectors  per 
cell  to  six.  The  power-controlled  system  operating  in  a  =  7  shadowing  environment 

could  support  up  to  70  users  per  cell  using  60°  antenna  sectoring  without  operating  in  the 
interference-limited  region  as  shown  in  Figure  5.6. 


9 

e 


10 


10 


10 


10 


icf4 


10 


10 
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Figure  5.6  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (crdB  =  7)  vs.  Users  per  Cell  using  Linear  User  Distribution 

with  60°  Sectoring  with  and  FEC  (Rcc  =  1/2  and  v=8) . 


The  use  of  sectoring  to  alleviate  the  co-channel  interference  in  support  of  power 
control  only  works  up  to  a  point.  For  a  shadowing  environment  of  =  8 ,  for  example, 
the  6-sector  system  can  accommodate  up  to  roughly  27  users  prior  to  becoming 
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interference-limited  at  an  SNR  per  bit  of  yb  =  15  dB,  which  could  be  an  acceptable  user 
capacity.  In  an  environment  where  adB  =  9 ,  however,  even  with  60°  sectoring  we  are 

interference  limited  with  as  few  as  12  operators  working  at  yb  =15.  Accordingly, 

implementing  power  control  in  such  environments  demands  some  additional  method  of 
interference  mitigation  for  the  system  to  be  viable.  Appendix  V-B  contains  complete 
performance  results  for  the  power-controlled  system  operating  in  a  lognormal 
environment  of  >  6  using  60°  sectoring. 

Accordingly,  we  have  shown  that  by  carefully  adding  power  control  to  the 
forward  traffic  channel,  we  can  greatly  improve  the  performance  of  our  DS-CDMA 
cellular  system  operating  in  a  Rayleigh-lognormal  channel. 

C.  SUMMARY 

In  this  chapter,  we  improved  upon  the  performance  analysis  of  Chapter  IV  for  a 
DS-CDMA  cellular  system  operating  in  a  Rayleigh  fading  and  lognormal  shadowing 
environment,  by  randomizing  user  positions  in  the  cell,  rather  than  assuming  the  worst- 
case  scenario.  We  presented  a  revised  upper  bound  for  the  probability  of  bit  error,  which 
incorporated  the  user  distribution  within  the  cell.  We  analyzed  the  performance  of  the 
system  using  a  linear  user  distribution  based  upon  the  density  of  users  being  uniform 
throughout  the  cell.  We  found  that  the  performance  of  the  system  improved  by  3  to  4  dB 
for  our  test  case  of  odB  =  7. 
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In  Section  V.B,  we  implemented  a  power-control  technique  that  ensures  the 
power  received  by  all  users  is  equal.  In  doing  so,  we  randomized  the  transmit  power  in 
each  user  channel  to  overcome  the  effects  of  lognormal-shadowing  and  distance 
dependent  path  loss.  We  found  that,  with  proper  implementation,  the  power-controlled 
system  can  greatly  increase  performance  and  user  capacity  of  the  forward  channel.  A 
system  designer  must,  however,  steer  clear  of  the  interference-limited  region  in  heavily 
shadowed  environments. 
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APPENDIX  V-A  GEOMETRY  OF  PATH  LOSS  DISTANCES 


In  this  appendix  we  develop  the  path  loss  distances  from  the  mobile  user  to  each 
of  the  adjacent  base  stations  in  terms  of  %  and  0.  In  all  cases,  we  use  the  geometry  of  the 
triangle  depicted  in  Figure  5.7  in  order  to  determine  S)r 


Base  Station  i 


Center  Cell 
Base  Station 

Figure  5.7  Key  Triangle  of  Relationships. 

We  define  0,  as  the  angle  between  the  mobile  user  and  base  station  i.  Using  side- 
angle  relations  in  plane  triangles,  we  determine  that 

fZ)2  =  *2  +  3  -  2^3*  COS  0,.  (5.24) 

Accordingly,  we  seek  to  determine  0;  as  a  function  of  0 ,  which  ensures  that  each  5)  is 
only  a  function  of  k  and  0. 
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1. 


The  Development  of  S\ 


Figure  5.8  Geometry  for  . 


In  determining  0j ,  we  must  consider  two  cases: 


Case I:  O<0<7t  (Figure  5.8),  0!=©. 

Case  II:  -rt<0<O  01=-0. 

Since  COS0  =  cos(-0) ,  we  take  0j  =  0  to  determine  from  (5.24)  as  follows: 

2\=yft^~  2V3^cos0!+3  (52 

=  V*2  -2V3*cos0  +  3. 
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2.  The  Development  of  S)2 


Figure  5.9  Geometry  for  S>2 . 
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3.  The  Development  of  S)3 


Figure  5.10  Geometry  for  2)3 . 
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In  determining  03 ,  we  must  consider  four  cases: 


Ott  Ott 

Case  I:  0  <  0  <  —  (Figure  5.10)  03  =  —  -  0  =  - 
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Case  II: 


2n 
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e,  =  8-— . 
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,  we  can  take  0,  =0-—  to  determine 


5X  from  (5.24)  as  follows: 

S)3  =yff- 2V3ycos0^T3 


=  .L2  2V3*cos  f  0 

V  l  3 


+  3 
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2n 


.  2n 


cos  0  cos —  +  sin  0sin — 
3  3 


+  3 


=  V*2+>/3*cos0-3*sin0  +  3. 


(5.27) 
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4.  The  Development  of  2)4 


Figure  5.11  Geometry  for  f£>4 . 

In  determining  04 ,  we  must  consider  two  cases: 

Case  I:  O<0<7t  (Figure  5.11)  04=7t-0. 

Case  II:  -7t<0<O  04=rt  +  0. 

Since  cos(7t  +  0)  =  COs(ti-0)  =  -cos0,  we  take  COS04  =-cos0  to  determine 
3)a  from  (5.24)  as  follows: 

®4  = 


ft2  -2^*.cos04  +3 
U2  +  2>/3*.cos0  +  3. 


(5.28) 


183 


5.  The  Development  of  S>s 


Figure  5.12  Geometry  for  5). 


In  determining  05 ,  we  must  consider  four  cases: 


Case  I: 


0<  0<  — 
3 


05  =  0  + 


27t 


71  I  27 1 

Case  II:  —  <  0  <  n  (Figure  5.12)  05  =  -  0  +  — 


Case  III: 


<eso 
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e,=e+— . 
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Case  IV: 
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Since  cos 


(  2n^ 

0  +  — — 

=  cos 

f  9+— Y) 

3  J 

\ 

l  3  JJ 

,  we  can  take  05  =0  +  —  to  determine 


%  from  (5.24)  as  follows: 


fZ>5  =^2  -2>/3»cos05  +3 


cos  0  cos - sin0sin — 

V  3  3  J 


=  V^+V^cos  0  +  3*  sin  0  +  3. 


6.  The  Development  of  S)6 
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In  determining  06 ,  we  must  consider  four  cases: 


Case  I: 

Case  H: 

Case  III: 

Case  IV: 

Since  cos 

(5.24)  as  follows: 


271  7T 

0  <  0  <  —  (Figure  5.13)  06=0  +  -^. 
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,  we  take  06  =  0+—  to  determine  S>6  from 


5)6  =  -  2^3*  cos  06  +  3 


=  Jx2  -2S% cos 


+  3 


=  U2-2-j3% 


COS0COS—  —  sin0sin— 


+  3 


(5.29) 


=  V*2  -\/3*cos0  +  3*.sin0  +  3. 
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APPENDIX  V-B  PROBABILITY  OF  BIT  ERROR  FOR  THE  RAYLEIGH- 
LOGNORMAL  CHANNEL  USING  FORWARD  CHANNEL 
POWER  CONTROL 


Figure  5.14  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( <rdB  =2)  using  Power  Control  and  120°  Sectoring 

(Rcc  =  1/2  and  v=8) . 
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SNR  per  bit  (dB) 


Figure  5.15  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( <rdB  =3)  using  Power  Control  and  120°  Sectoring 

(Rcc  =l/2andv=8). 


SNR  per  bit  (dB) 


Figure  5.16  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( crdB  =  4)  using  Power  Control  and  120°  Sectoring 

(Rcc  =  1/2  and  v=8) . 


Figure  5.17  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (<rdB  =5)  using  Power  Control  and  120°  Sectoring 

(Rcc  =  l/2andv=8). 


Figure  5.18  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( <rdB  =  6)  using  Power  Control  and  120°  Sectoring 

(Rcc  =  1/2  and  v=8) . 


SNR  per  bit  (dB) 


Figure  5.19  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( adB  =7)  using  Power  Control  and  60°  Sectoring 

(Rcc  =  1/2  and  v=8) . 


SNR  per  bit  (dB) 

Figure  5.20  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  (cr^  =  8)  using  Power  Control  and  60°  Sectoring 

(Rcc  =  1/2  and  v=8) . 
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Figure  5.21  Probability  of  Bit  Error  for  DS-CDMA  with  Rayleigh  Fading  and 
Lognormal  Shadowing  ( crdB  =9)  using  Power  Control  and  60°  Sectoring 

(Rcc  =  1/2  and  v=8) . 
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VI.  CONCLUSIONS  AND  FUTURE  WORK 


In  this  thesis,  we  have  explored  and  developed  many  tools  for  analyzing  the 
performance  of  the  forward  channel  of  a  sectored  DS-CDMA  cellular  system  with  EEC 
operating  in  a  Rayleigh  fading  and  lognormal  shadowing  environment.  We  will  now 
summarize  the  work  and  conclusions  and  outline  areas  of  future  work. 

j 

A.  CONCLUSIONS 

In  Chapter  n,  we  examined  the  algebraic  properties  and  generation  of  a  set  of 
Walsh  functions,  and  proved  the  closure  property  for  a  set  of  Walsh  functions.  We 
developed  the  notion  of  extending  the  orthogonality  property  for  a  reduced  set  Walsh 
functions  so  that  product  of  any  three  Walsh  functions  from  the  reduced  set  is  identically 
zero  when  integrated  over  the  bit  period.  We  proposed  a  method  to  generate  this  reduced 
set  and  determined  the  maximum  number  of  Walsh  functions  in  the  reduced  set  that 
ensures  extended  orthogonality.  We  concluded  the  chapter  by  developing  the 
autocorrelation  functions  for  a  set  of  Walsh  functions  and  by  describing  the  average 
autocorrelation  function  for  the  entire  set. 

In  Chapter  HI,  we  developed  a  forward  channel  model  for  the  DS-CDMA  cellular 
system  that  combined  the  effects  of  large-scale  and  small-scale  propagation  into  a  single 
unified  model.  We  used  the  extended  Hata  model  to  predict  large-scale  path  loss  and 
incorporated  lognormal  shadowing.  We  also  included  small-scale  propagation  effects 
into  the  model  using  slow-flat  Rayleigh  fading.  We  developed  the  composite  signal 
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received  by  a  mobile  user  in  the  worst-case  position  consisting  of  the  signals  from  the 
center  and  adjacent  base  stations,  which  have  traveled  through  the  Rayleigh-lognormal 
channel. 

In  Chapter  IV,  we  developed  a  Signal-to-Noise  plus  Interference  Ratio  (SNIR) 
and  the  bit-error  probability  for  the  DS-CDMA  cellular  system  in  a  Rayleigh-lognormal 
channel  by  using  Gaussian  approximation.  We  incorporated  Forward  Error  Correction 
(EEC)  in  the  form  of  convolutional  encoding  with  soft-decision  decoding  into  our 
forward  channel  model  and  develop  a  revised  upper  bound  on  the  probability  of  bit  error. 
We  approximated  a  sum  of  d  multiplicative  chi-square(two  degrees  of  freedom)- 
lognormal  random  variables  as  a  multiplicative  chi-square(with  2d  degrees  of  freedom)- 
lognormal  random  variable,  in  order  to  evaluate  the  probability  of  bit  error  for  the  DS- 
CDMA  system  with  FEC  operating  in  the  Rayleigh-lognormal  channel.  Our  modeled 
performance  results  compared  favorably  with  the  Monte  Carlo  simulated  results. 
Specifically,  we  found  that  the  model  was  quite  accurate  in  the  SNR  per  bit  range  of 

yb  = 10  to  15  dB.  In  predicting  the  interference  floor,  the  model  effectively  upper 

bounded  the  Monte  Carlo  simulated  results.  Using  our  model  to  approximate  the  sum  of 
d  multiplicative  chi-square(two  degrees  of  freedom)-lognormal  random  variables  in  the 
performance  analysis  of  DS-CDMA  cellular  systems  operating  in  a  Rayleigh-lognormal 
channel  is  a  fast,  effective  and,  accordingly,  a  very  valuable  tool  for  system  designers. 

We  concluded  Chapter  IV  by  incorporating  antenna  sectoring  into  our  forward 
channel  model.  Lastly,  we  analyzed  the  performance  of  our  forward  channel  model 
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operating  in  the  Rayleigh-lognormal  channel  under  a  wide  range  of  shadowing  conditions 
and  for  various  numbers  of  users  per  cell. 

In  Chapter  V,  we  modified  the  performance  analysis  developed  in  Chapter  IV  to 
incorporate  user  distribution  within  the  cell  by  assuming  that  the  user’s  position  is 
random  according  to  a  specified  distribution.  We  revised  the  upper  bound  on  the 
probability  of  bit  error  for  the  DS-CDMA  cellular  system  with  EEC  in  a  Rayleigh- 
lognormal  channel  to  accommodate  distributed  users.  Accordingly,  we  developed  a  more 
realistic  view  of  typical  performance  for  our  forward  channel  model.  We  analyzed  the 
performance  of  the  system  using  a  linear  user  distribution  and  found  that  system 
performance  was  3  to  4  dB  better  than  predicted  by  our  worst-case  analysis. 

We  concluded  Chapter  V  by  introducing  forward  channel  power  control  as  a 
means  of  enhancing  the  overall  performance  of  the  system.  We  analyzed  the 
performance  of  the  system  with  power  control  under  a  wide  range  of  shadowing 
conditions  and  for  various  numbers  of  users  per  cell.  We  found  that  in  light  to  medium 
shadowing  environments  (odB  <  6)  power  control  can  greatly  enhance  the  performance 

of  the  forward  channel  and  accordingly  increase  user  capacity.  In  heavily  shadowed 
environments,  performance  is  enhanced  provided  that  the  system  is  not  operating  in  the 
interference-limited  region,  which  is  a  function  of  required  user  capacity  and  the 
shadowing  parameter  G^. 
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B.  FUTURE  WORK 


There  is  much  opportunity  for  further  research  in  the  area  of  DS-CDMA  cellular 
systems  operating  in  a  Rayleigh  fading  and  lognormal  shadowing  channels.  For  example, 
implementation  of  FEC  using  Serially  Concatenated  Convolutional  Codes  (SCCC)  would 
enhance  our  analysis.  The  use  of  concatenated  codes  with  interleavers  can  improve  the 
performance  of  the  cellular  system  by  overcoming  any  burst  errors  that  may  arise  due  to 
multipath  propagation  effects. 

Our  analysis  in  Chapter  V  incorporated  a  user  distribution  based  on  uniform  user 
density  within  the  cell.  There  are  many  other  user  densities  and  distributions  that  could 
be  explored  and  that  may  better  reflect  user  distribution  in  practice. 

In  implementing  fast  power  control  in  Chapter  V,  we  assumed  that  the  base 
station  was  not  limited  in  the  total  amount  of  power  transmitted.  That  is,  in  order  to 
achieve  the  target  receive  power  level  for  every  user  and  overcome  shadowing  regardless 
of  severity,  the  base  station  was  free  to  boost  power  to  whatever  level  was  required.  Of 
course  in  practice,  the  base  station  would  be  limited  in  total  transmit  power. 
Accordingly,  an  analysis  that  caps  the  base  station’s  total  transmit  power  and  considers 
power  allocation  to  the  service  channels  would  be  insightful. 

Finally,  the  analysis  of  the  forward  channel  developed  in  this  thesis  should  be 
extended  to  the  reverse  channel  of  a  DS-CDMA  cellular  system  in  a  Rayleigh  fading  and 
lognormal  shadowing  environment. 
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